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Part  of  the  913-detector  array  of  the  new  Disk  Chopper  Spectrometer 
(DCS)  that  will  be  available  to  users  in  2000.  To  view  the  detectors  as 
found  in  the  instrument,  the  picture  should  be  rotated  90°.  Information 
about  the  DCS  and  other  instruments  at  NCNR  is  available  at 
http://www.ncnr.nist.gov/. 
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FOREWORD 


Once  again,  it  is  a  pleasure  to  be  able  to  reflect  on  the  accom- 
plishments of  the  NIST  Center  for  Neutron  Research  over 
the  past  year.  In  the  reactor  operations  area,  1999  was  another 
outstanding  year.  In  spite  of  an  unsheduled  maintenance  shutdown, 
the  reactor  operated  250  days,  with  a  reliability  factor  of  better  than 
90  %.  The  cold  source  availability  for  the  period  was  98  %;  i.e.,  the 
cold  source  held  the  reactor  from  operation  4  days  during  the  year. 
The  remaining  spent  fuel  in  the  storage  pool  was  shipped,  providing 
space  for  at  least  five  years  of  operation.  Also,  an  order  has  been 
placed  for  a  new  cooling  tower  which  will  not  only  provide  needed 
capability  for  the  next  25  years,  but  will  also  reduce  the  plume 
visible  during  cold  weather.  Last,  steady  progress  has  been  made  in 
preparing  for  a  license  renewal  application  to  the  Nuclear  Regulatory 
Commission,  in  order  to  extend  the  period  of  operation  beyond  2004. 

We  have  also  made  great  progress  in  instrumentation,  with 
the  back  scattering  spectrometer  operational;  with  the  spin  echo 
spectrometer  now  being  used  for  real  measurements;  and  with  the 
disk  chopper  spectrometer  being  commissioned.  All  three  of  these 
instruments  will  be  available  to  "friendly"  users  in  the  next  proposal 
cycle.  Other  work  is  also  advancing  well-the  perfect  crystal  small 
angle  scattering  spectrometer  (part  of  the  NSF/NIST  CHRNS)  is 
being  installed  at  the  reactor;  the  first  phase  of  the  high  intensity 
filter  analyzer  spectrometer  is  ready  to  begin  installation;  and  the 
design  and  manufacture  of  new  thermal  neutron  spectrometers  is 
underway.  This  simultaneous  development  program  has  put  severe 


strains  on  our  resources,  but  we  can  now  look  forward  to  many  years 
of  benefit  from  the  results. 

Finally,  as  always,  the  results  are  seen  in  the  output  of  the 
researchers  who  use  the  facility.  As  was  done  last  year,  we  are 
presenting  highlights  of  this  work  in  the  following  chapters  of  this 
report.  I  think  all  can  agree  that  the  results  truly  speak  for  them- 
selves. 
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INTRODUCTION  TO 
THE  NIST  CENTER  FOR 
NEUTRON  RESEARCH  (NCNR) 


Modern  technological  society  is  dependent  upon  increasingly 
sophisticated  use  of  materials,  many  of  whose  attributes  are 
dictated  by  their  sub-microscopic  structural  and  dynamical  proper- 
ties. Our  knowledge  of  these  properties  is  provided  by  a  wide  range 
of  scientific  techniques  of  which  the  many  types  of  scattering  (for 
example.  X-rays,  light,  electrons,  neutrons)  are  arguably  the  most 
important.  Of  these  probes,  neutrons  are  perhaps  least  familiar,  but 
they  provide  important  advantages  for  many  types  of  measurements. 

Neutrons,  as  prepared  for  use  at  modern  sources,  are  moving 
at  speeds  comparable  to  those  of  atoms  moving  at  room  temperature, 
thus  providing  the  ability  to  probe  dynamical  behavior.  At  the  same 
time,  neutrons  are  well  matched  to  measurements  at  length  scales 
ranging  from  the  distances  between  atoms  to  the  size  of  biological 
or  polymer  macromolecules.  Neutrons  are  sensitive  to  the  magnetic 
properties  of  atoms  and  molecules,  allowing  study  of  the  underlying 
magnetic  properties  of  materials.  They  also  scatter  quite  differently 
from  normal  hydrogen  atoms  than  they  do  from  heavy  hydrogen 
(deuterium),  allowing  selective  study  of  individual  regions  of  molec- 
ular systems.  Finally,  neutrons  interact  only  weakly  with  materials, 
providing  the  opportunity  to  study  samples  in  different  environments 
more  easily  (at  high  pressures,  in  shear,  in  reaction  vessels,  etc.),  and 
making  them  a  non-destructive  probe.  These  favorable  properties  are 
offset  by  the  relative  weakness  of  the  best  neutron  sources  compared 
to  X-ray  or  electron  sources,  and  by  the  relatively  large  facilities 
required  to  produce  neutrons.  As  a  result,  major  neutron  sources 
are  operated  as  national  user  facilities  to  which  researchers  come 
from  all  over  the  United  States  (and  abroad)  to  perform  small  scale 
science  using  the  special  measurement  capabilities  provided. 

In  addition  to  scattering  measurements,  neutrons  can  be  used 
to  probe  the  atomic  composition  of  materials  by  means  of  capture 
and  resultant  radioactive  decay.  The  characteristics  of  the  decay  act 
as  "fingerprints"  for  particular  atomic  nuclei,  allowing  studies  of 
environmental  samples  for  pollutants  (e.g.,  heavy  metals),  character- 
ization of  Standard  Reference  Materials,  and  many  other  essential 
measurements.  While  the  scattering  and  capture  users  of  neutrons 
are  little  concerned  with  understanding  the  inherent  properties  of  the 
neutron,  there  are  important  areas  in  physics  that  can  be  explored 
by  carefully  measuring  fundamental  neutron  behavior.  Examples 
include  the  lifetime  of  the  free  neutron,  an  important  quantity  in  the 
theory  of  astrophysics;  the  beta  decay  process  of  the  neutron,  the 


details  of  which  are  stringent  tests  of  nuclear  theory;  and  the  effects 
of  various  external  influences  such  as  gravity  or  magnetic  fields  on 
neutrons. 

The  NCNR  utilizes  neutrons  produced  by  the  20  MW  NIST 
Research  Reactor  to  provide  facilities,  including  the  nation's  only 
internationally  competitive  cold  neutron  facility,  for  all  of  the  above 
types  of  measurements  to  a  national  user  community.  There  are 
approximately  35  stations  in  the  reactor  and  its  associated  beams  that 
can  provide  neutrons  for  experiments.  At  the  present  time  26  of  these 
are  in  active  use,  of  which  6  provide  high  neutron  flux  positions  in 
the  reactor  for  irradiation,  and  20  are  beam  facilities.  A  schematic 
layout  of  the  beam  facilities  and  brief  descriptions  of  available 
instrumentation  are  given  below.  More  complete  descriptions  can  be 
found  at  http://www.ncnr.nist.gov/. 

These  facilities  are  operated  both  to  serve  NIST  mission  needs 
and  as  a  national  facility,  with  many  different  modes  of  access.  Some 
instrumentation  was  built  years  ago,  and  is  not  suited  to  general 
user  access;  however,  time  is  available  for  collaborative  research. 
NIST  has  recently  built  new  instrumentation,  and  reserves  33  %  of 
available  time  for  mission  needs  with  the  balance  available  to  general 
users.  In  other  cases,  instrumentation  was  built  and  is  operated 
by  Participating  Research  Teams  (PRT).  PRT  members  have  access 
to  75  %  of  available  time,  with  the  balance  available  to  general 
users.  In  a  special  case,  NIST  and  the  National  Science  Foundation 
established  the  Center  for  High  Resolution  Neutron  Scattering  at  the 
NCNR,  with  a  30  m  Small  Angle  Scattering  (SANS)  instrument,  a 
cold  neutron  triple  axis  spectrometer,  and  a  perfect  crystal  SANS 
under  construction.  For  these  facilities,  most  time  is  available  for 
general  users.  While  most  access  is  for  research,  whose  results  are 
freely  available  to  the  general  public,  proprietary  research  can  be 
performed  under  full  cost  recovery.  Each  year,  approximately  1600 
researchers  (persons  who  participated  in  experiments  at  the  facility, 
but  did  not  necessarily  come  here)  from  all  areas  of  the  country,  from 
industry,  academe,  and  government  use  the  facility  for  measurements 
not  otherwise  possible.  The  research  covers  a  broad  spectrum  of 
disciplines,  including  chemistry,  physics,  biology,  materials  science, 
and  engineering. 
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NIST  CENTER  FOR 
NEUTRON  RESEARCH 
LAYOUT 


1  Horizontal  Sample 
Reflectometer 

Horizontal  surface  of  sample 
allows  reflectivity  measure- 
ments of  free  surfaces,  liquid 
vapor  interfaces,  as  well  as 
polymer  coatings. 

2  Neutron  Interferometry  and 
Optics  Station 

Includes  perfect  silicon  inter- 
ferometer; vibration  isolation 
system  provides  exceptional 
phase  stability  and  fringe 
visibility. 


3  Prompt  Gamma  Activation 
Analysis 

Cold  neutron  fluxes  allow 
detection  limit  for  H  of  1  i^g 
to  10  (jug.  Focused  beams  are 
available  for  profiling. 

4  NG-7  30  m  SANS 

Small  Angle  Neutron  Scattering 
instrument  for  microstructure 
measurement  sponsored  by 
NIST,  the  ExxonMobil  Research 
and  Engineering  Co.,  the 
University  of  Minnesota,  and 
Texaco  R&D. 

5  Neutron  Physics  Station 

A  cold  neutron  beam 
150  mm  x  60  mm,  available 
for  fundamental  neutron  phys- 
ics experiments. 


6  Fermi  Chopper TOF 
Spectrometer 

A  hybrid  time-of-flight  spec- 
trometer for  inelastic  scatter- 
ing, with  wavelengths  between 
0.23  nm  and  0.61  nm.  The 
wavelength  is  chosen  by  a 
focusing  pyrolytic  graphite 
crystal,  while  the  beam  is 
pulsed  by  a  simple  Fermi 
chopper. 

7  Spin  Echo  Spectrometer 

A  neutron  spin  echo  spec- 
trometer offering  neV  energy 
resolution,  based  upon  Julich 
design,  sponsored  by  NIST, 
Jiilich  and  ExxonMobil. 
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10  NG-3  30  m  SANS 

Instrument  for  microstructure 
measurement  sponsored  by  the 
National  Science  Foundation 
and  NIST;  part  of  CHRNS. 

11  Backscattering  Spectrometer 

High  intensity,  very  high 
resolution  backscattering 
spectrometer,  with  many 
innovative  features,  and  energy 
resolution  of  approximately 
1  |xeV. 

12  8  m  SANS 

Instrument  for  polymer 
characterization,  sponsored  by 
Polymers  Division. 


13  Vertical  Sample  Reflectometer 

Instrument  for  measuring 
subsurface  structure  with 
polarization  analysis  capability. 
Reflectivities  down  to  10"8  can 
be  measured. 

14  BT-8  Residual  Stress 
Diffractometer 

Diffractometer  optimized  for 
depth  profiling  of  residual 
stress  in  large  components. 

15  BT-9  Triple  Axis  Spectrometer 

Triple  axis  crystal  spectrometer 
for  measurements  of  excita- 
tions and  structure. 


16  BT-1  Powder  Diffractometer 

Powder  diffractometer  with  32 
detectors;  incident  wavelengths 
of  0.208  nm,  0.154  nm,  and 
0.159  nm,  with  highest  resolu- 
tion of  5d/d  =  8x1 04. 

17  BT-2  Triple  Axis  Spectrometer 

Triple  axis  crystal  spectrometer 
with  polarized  beam  capability 
for  measurement  of  magnetic 
dynamics  and  structure. 

18  BT-4  Filter  Spectrometer 

A  triple  axis  crystal 
spectrometer  with  a  Be  or 
Be/Graphite  filter  analyzer 
option  for  chemical 
spectroscopy. 


19  Cold  Neutron  Depth  Profiling 

A  station  for  quantitative 
profiling  of  subsurface 
impurities  and  coatings,  based 
on  neutron  capture  and 
emission  of  a  charged  particle. 

20  Thermal  Column 

A  very  well-thermalized  beam 
of  neutrons  used  for  radiog- 
raphy, tomography,  dosimetry, 
and  other  experiments. 
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NCNR  IMAGES  1999 
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1.  An  incoming  SANS  image  excites 
the  interest  of  Martin  Vigild  (center) 
as  Frank  Bates  (standing),  Newell 
Washburn  (right)  and  Ken  Hanley, 
all  of  the  University  of  Minnesota, 
look  on. 

2.  NCNR's  Peter  Gehring  describes  the 
analyzer  of  the  spin  polarized  triple 
axis  spectrometer  (SPINS)  to  NCNR 
Summer  School  participants. 


3.  The  Spin  Echo  Spectrometer  (NSE) 
at  NCNR,  commissioned  this  year,  is 
described  in  a  Research  Highlight  in 
this  issue. 

4.  Stephen  FitzGerald  (Oberlin 
College)  loads  a  sample  at  the  HFBS. 

5.  Gudrun  Schmidt  (NIST,  Polymers 
Division)  awaits  the  display  of  an 
updated  SANS  image. 


6.  NCNR  Summer  School  participants 
interacting  at  the  recently 
commissioned  high-flux  back- 
scattering  spectrometer  (HFBS). 

7.  Silke  Rathgeber  and  So  Hyun  Park 
(both  at  NCNR),  ready  to  load  a 
sample  at  the  NSE. 


10.  Dave  Clem  and  Mike  Rinehart 
(both  of  NCNR)  installing  biological 
shielding  at  the  DCS. 
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SUPERLATTICE 
MAGNONS 


A powerful  technique  for  optimizing  material  properties  is  to 
deposit  alternating  layers  of  different  materials  to  form  a  thin 
film  superlattice.  In  particular,  magnetic  and  non-magnetic  materi- 
als grown  as  a  superlattice  exhibit  a  variety  of  tunable  couplings 
between  the  magnetic  layers.  Control  of  these  couplings  by  varying 
the  layer  materials  and  thicknesses  has  lead  to  dramatic  increases  in 
performance  in  certain  device  applications.  A  prime  example  is  the 
recent  use  of  layered  magneto-resistive  films  in  high-performance 
magnetic  recording  heads  and  sensors.  The  inter-layer  magnetic 
coupling  can  survive  across  as  many  as  30  non-magnetic  atomic 
planes  due  to  its  one-dimensional  nature,  but  the  coupling  strength 
is  always  much  weaker  than  bulk  material  magnetic  interactions. 

The  behavior  of  the  magnetic  fluctuations  (magnons)  in  these 
weakly  coupled  layer  systems  is  of  interest  because  it  provides  a 
direct  measure  of  the  magnetic  interac- 
tions responsible  for  the  magnetic  struc- 
ture, and  leads  to  a  better  understanding 
of  the  unique  layer  to  layer  couplings. 
Until  now,  the  only  measurements  that 
have  probed  these  fluctuations  have 
used  Brillouin  light  scattering  [1],  or 
ferromagnetic-resonance  techniques  [2], 
both  of  which  measure  only  the  longest 
wavelength  dynamics.  These  measure- 
ments have  found  interesting  resonances 
associated  with  the  superlattice 
structure. 

It  is  important  to  directly  measure 
the  magnons  at  shorter  wavelengths  in 
order  to  determine  the  dispersion,  which 
directly  relates  to  the  magnetic  interac- 
tions, both  within  and  between  layers, 
and  whether  or  not  the  magnetic  fluc- 
tuation waves  propagate  between  layers. 
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1.70 


Only  inelastic  magnetic  scattering  techniques  can  provide  this  infor- 
mation, but  with  current  neutron  scattering  sources  the  intensities 
from  such  measurements  can  be  prohibitively  low,  because  the 
amount  of  magnetic  material  in  the  films  is  so  small. 

In  order  to  overcome  this  difficulty  we  have  made  neutron 
inelastic  scattering  measurements  on  a  very  large  superlattice  of 
alternating  layers  of  dysprosium  and  yttrium.  Dysprosium  is  the 
magnetic  constituent,  and  it  has  the  strongest  neutron  magnetic- 
scattering  of  all  the  elements.  In  order  to  maximize  the  amount 
of  Dy,  350  bilayers  composed  of  43  A  of  Dy  and  28  A  of  Y 
(designated  [Dy43/Y,g]350)  were  grown  by  MBE  techniques  on  a  2.5 
cm  x  1.3  cm  substrate  resulting  in  3  mg  of  Dy  in  the  sample. 

The  magnetic  structure  of  this  superlattice  has  previously 
been  determined  to  be  the  basal-plane  helical  structure  of  bulk 
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FIGURE  1.  Inelastic  magnetic  scattering  from  a  [Dy43/Y28]350  superlattice, 
obtained  by  subtracting  10  K  data  from  75  K  data,  shown  as  an  intensity  map 
in  Q-energy  space.  The  highest  intensity  is  300  counts/  30  min  and  decreases 
by  30  counts  for  each  level.  A  magnetic  Bragg  peak  is  just  off  the  graph  at 
Q  =  1.97  A1. 
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Dy  with  a  coherence  length  greater  than  4  bilayers.  The  helix 
progresses  remarkably  undisturbed  through  the  conduction  electrons 
of  the  non-magnetic  yttrium,  but  with  a  turn-angle  that  is  different 
than  in  the  dysprosium.  The  coupling  strength  across  the  yttrium 
has  been  measured  by  applying  a  uniform  external  magnetic  field. 
The  coupling  breaks  down  when  the  external  field  provides  0.2  meV 
per  Dy  atom  in  the  basal-plane,  while  the  equivalent  zone-boundary 
magnetic  fluctuation  in  bulk  Dy  is  greater  than  5  meV  per  Dy  atom. 

The  inelastic  neutron  scattering  measurements  were  per- 
formed at  the  NIST  Center  for  Neutron  Research  on  the  cold- 
neutron  spectrometer  SPINS.  A  multi-component  crystal  was  used 
to  analyze  the  energy  of  the  scattered  neutrons  in  order  to  enhance 
the  measured  intensity  while  sacrificing  wave-vector  (Q)  resolution. 
Measurements  were  performed  at  75  K  as  a  compromise  between 
the  size  of  the  ordered  magnetic  moment  and  the  thermal  population 
of  magnons.  The  magnons  of  interest  are  those  propagating  along 
the  superlattice  growth  direction  or  the  c-axis  of  the  hep  rare-earth 
structure.  Lower  magnon  energies  produce  higher  thermal  popula- 
tions, but  become  contaminated  with  elastic  background  scattering, 
so  it  is  necessary  to  subtract  scans  taken  at  low  temperatures  where 
the  magnons  have  become  depopulated  in  order  to  remove  this 


background.  The  resulting  magnetic  signal  is  shown  in  Fig.  1  as 
a  color-coded  intensity  map  in  Q-energy  space.  There  is  a  clear 
ridge  of  intensity  which  moves  to  higher  energies  as  Q  moves  away 
from  the  magnetic  Bragg  peak  at  1.97  A"1.  We  have  concentrated 
on  the  magnon  branches  that  extend  towards  smaller  Q  since  they 
move  away  from  the  intense  magnetic  and  nuclear  Bragg  peaks 
that  produce  a  large  elastic  background.  This  measured  dispersion  is 
compared  with  bulk  Dy  in  Fig.  2,  which  also  shows  the  diffraction 
pattern  along  the  growth  direction  (c-axis)  of  the  superlattice.  The 
bulk  Dy  dispersion  is  shown  as  lines  originating  from  the  magnetic 
Bragg  peaks  at  Q  =  1.97  A1  and  Q  =  2.06  A1.  The  measured 
magnons  are  shown  as  bars  centered  on  the  measured  peak  posi- 
tions and  with  lengths  representing  the  full-width-at-half-maximum- 
intensity.  The  agreement  with  the  bulk  dispersion  is  quite  good. 
There  is  no  evidence  in  these  data  of  the  influence  of  the  yttrium 
layers  other  than  possibly  the  splitting  of  the  dispersion  into  two 
branches  because  of  the  superlattice  structure  as  shown  by  separate 
bars  at  both  Q  =  1.8  A1  and  at  Q  =  1.75  A"1.  The  observed  modes 
are  not  over-damped,  but  we  cannot  measure  the  damping  under 
the  current  experimental  conditions.  Also,  these  data  would  have 
to  be  extended  to  smaller  Q,  in  order  to  approach  the  interface 
thickness.  We  are  currently  designing  additional  focusing  configura- 
tions so  that  the  instrumental  resolution  will  be  better  optimized  for 
measurements  of  this  dispersion  surface. 

REFERENCES 

[1]  B.  Hillebrands  and  G.  Gtintherodt  in  Ultrathin  Magnetic  Structures  I  +  //. 
edited  by  B.  Heinrich  and  J.  A.  C.  Bland  (Springer  \ferlag,  1994). 

[2]  C.  F.  Majkrzak  et  al.,  Adv.  in  Phys.  40,  99  (1991). 


FIGURE  2.  The  dispersion  along  the  growth  direction  (c-axis)  measured  for  the 
[Dy^YJ^  superlattice  is  plotted  against  the  energy  scale  on  the  left  side. 
It  is  compared  to  the  measured  dispersion  in  bulk  dysprosium  (solid  lines) 
expected  to  originate  from  each  of  the  superlattice  magnetic  Bragg  peaks. 
(Only  the  negative  Q  branches  from  the  two  low  Q  Bragg  peaks  are  shown). 
On  the  right  side  the  diffraction  scan  for  this  sample  showing  the  superlattice 
peaks  is  displayed.  This  film  is  so  large  that  the  sample  diffraction  peaks  are 
as  strong  as  the  substrate  peak  shown  in  blue  on  the  right. 
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ANTIFERROMAGNETIC  SPIN 
REORIENTATION  IN  EXCHANGE-BIASED 

Fe304/NiO  SUPERLATTICES 


The  magnetic  hysteresis  loop  of  a  ferromagnetic  material  can  be 
displaced  along  the  field  axis  as  a  result  of  magnetic  exchange 
coupling  to  an  antiferromagnet.  This  "exchange-biasing"  phenom- 
enon was  discovered  over  40  years  ago  in  oxidized  Co  particles  [1] 
and  has  been  observed  in  a  variety  of  thin  films  and  multilayers. 
Artificial  spin-valve  structures  with  exchange-biased  layers  show 
great  promise  for  applications  as  magnetoresistive  sensors  in  read 
heads.  Recent  theories  [2-4]  make  specific  predictions  about  the  ori- 
gin of  exchange  biasing  and  the  response  of  the  antiferromagnetic 
layer. 

High-angle  neutron  diffraction  is  an  ideal  probe  as  the  anti- 
ferromagnet gives  rise  to  distinct  reflections  of  magnetic  origin. 
We  have  performed  neutron  diffraction  studies  of  exchange-biased 
[001]  Fe,0/NiO  superlattices  on  the  SPINS  and  BT-9  triple-axis 
spectrometers.  These  measurements  confirm  that  the  NiO  layers 
retain  their  bulk  antiferromagnetic  structure  in  which  ferromagnetic 
{111}  planes  alternate  direction  along  each  <111>  axis.  Our  data 
show  that  exchange  biasing  is  associated  with  domain  walls  that 
form  and  "freeze"  within  the  antiferromagnetic  NiO  layer.  Upon 
field  cooling  into  the  exchange-biased  state,  magnetic  domains  lock 
within  the  NiO  layers  and  do  not  change  with  subsequent  applica- 
tion of  magnetic  fields.  In  contrast,  the  antiferromagnetic  domain 
sizes  in  unbiased  samples  prepared  by  cooling  in  zero  field  depend 
sensitively  on  the  strength  of  the  applied  field. 

We  focus  here  on  a  Fe,04(6  nm)/NiO(  1 1  nm)  superlattice 
deposited  using  oxygen  plasma-assisted  molecular  beam  epitaxy. 
Measurements  of  the  magnetic  hysteresis  loop  show  little  evidence 
of  exchange  biasing  after  cooling  in  zero  field.  Strong  exchange 
biasing  is  induced  upon  field  cooling  the  superlattices  from  high 
temperatures.  For  our  superlattice,  the  biasing  field  at  30  K  is  0.043 
T  after  field  cooling  from  550  K. 

For  the  diffraction  experiments  the  sample  was  oriented  as 
shown  in  the  inset  of  Fig.  1.  Vertical  magnetic  fields,  H,  were 
applied  in  the  sample  plane.  In  this  configuration,  the  (1 1 1)  and 
(111)  NiO  antiferromagnetic  reflections  lie  in  the  scattering  plane. 
Using  a  horizontal-field  magnet,  the  (1 1 1)  and  (1 1 1)  reflections 


could  be  accessed  by  rotating  the  sample  90°  about  the  growth 
axis.  Figure  1  shows  a  typical  growth-axis  (00/)  scan  through  the 
(111)  reflection  for  the  superlattice.  Structural  stacking  faults  [5]  at 
the  NiO/Fe,0,  interfaces  limit  the  coherence  of  the  Fe,0,  structural 

3    4  3  4 

and  magnetic  order  to  a  single  Fe,04  layer.  As  a  result,  we  can 
easily  separate  the  broad  Fe304  and  narrow  NiO  contributions  to  the 
reflection  and  track  the  latter  as  a  function  of  field. 

The  neutron  scattering  data  reveal  that  the  magnetic  domain 
sizes  in  the  antiferromagnetic  NiO  depend  on  the  presence  or 
absence  of  exchange  biasing.  Figure  2  shows  the  full-width-at-half- 
maximum  (FWHM)  for  the  ( 1 1 1 )  NiO  peak  scanned  along  the  (00/) 
and  (//0)  directions  after  cooling  in  a  6  T  field  (i.e.,  exchange-biased 
state)  and  cooling  in  zero  field  (i.e.,  unbiased  state).  After  field 
cooling,  the  FWHM  of  the  NiO  reflection  scanned  along  the  (00/) 
growth  direction  shows  no  dependence  on  field. 


0.68      0.72        0.76        0.80  0.84 
QZ(A"1) 

FIGURE  1 .  Zero-field  growth-axis  (00/)  scan  through  the  (111)  reflection  for  the 
Fe304(6  nm)/Ni0(11  nm)  superlattice  after  cooling  to  78  K  in  a  6  T  field  parallel 
to  [110].  The  broad  Fe304  peak  is  shown  in  red,  and  the  remaining  scattering 
is  from  the  NiO.  The  green  arrow  marks  the  FWHM  of  the  NiO  peak.  The  inset 
shows  the  scattering  diagram. 
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The  antiferromagnetic  domain  size  along  the  growth  direc- 
tion, determined  from  the  inverse  of  the  FWHM,  remains  constant 
near  750  A  after  field  cooling.  This  contrasts  with  the  pronounced 
field-dependence  of  the  FWHM  observed  after  cooling  in  zero  field 
(Fig.  2).  In  this  state,  the  FWHM  is  smallest  near  zero  field,  but 
approaches  the  constant  field-cooled  value  when  the  magnitude  of 
the  field  is  greater  than  2  T  The  corresponding  domain  size  varies 
from  approximately  1200  A  in  zero  field  to  800  A  in  high  fields. 
Growth-plane  (7/0)  scans  through  the  NiO  reflection  show  a  similar 
difference  between  the  field  dependence  of  the  FWHM  after  cooling 
in  zero  field  and  cooling  in  a  6  T  field  (Fig.  2). 

Upon  field  cooling,  domain  walls  both  parallel  and  perpen- 
dicular to  the  growth  direction  lock  into  the  antiferromagnetic  layer, 
presumably  due  to  the  exchange  coupling  between  the  NiO  antifer- 
romagnetic moments  and  the  Fe,04  moments  that  are  aligned  paral- 
lel to  the  cooling  field.  After  cooling  in  zero  field,  the  domains  are 
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FIGURE  2.  Fu II- width-at-half -maximum  (FWHM)  of  the  (111)  NiO  reflection  as  a 
function  of  field  at  78  K  after  cooling  in  a  6  T  field  and  after  cooling  in  zero 
field.  The  top  plots  show  the  FWHM  of  the  peak  along  the  (00/)  direction.  The 
bottom  plots  show  the  FWHM  from  (110)  scans.  The  black,  blue  and  red  data 
are  from  each  field  cycle. 


larger,  but  reversibly  decrease  in  size  in  high  fields  (Fig.  2)  as  the 
magnetic  frustration  increases.  Consistent  with  several  theoretical 
predictions  [2,3]  we  believe  that  exchange  biasing  may  originate 
from  magnetic  frustration  that  leads  to  "frozen"  domain  walls  in  the 
field-cooled  state. 

We  also  observe  that  the  magnitude  of  the  ordered  NiO 
moments  in  all  four  of  the  {111}  domains  in  the  Fe304/NiO  super- 
lattices  depends  upon  applied  field.  After  cooling  in  a  large  vertical 
field,  the  intensity  of  the  antiferromagnetic  (111)  NiO  reflection 
reversibly  decreases  as  the  field  is  increased.  (The  behavior  is 
qualitatively  similar  after  cooling  in  zero  field.)  Analogous  field- 
cooled  measurements  of  the  (111)  NiO  peak  were  performed  in 
a  horizontal  field  and  surprisingly  show  a  comparable  intensity 
decrease.  Some  of  the  NiO  moments  thus  seem  to  disappear  out 
of  all  four  {111 }  domains.  We  speculate  that  these  NiO  spins 
may  become  disordered  as  a  result  of  the  high  internal  fields. 
Simultaneously,  some  of  the  NiO  moments  reorient  perpendicular 
to  the  Fe,04  magnetization  direction  in  high  fields  [6].  While  this 
"spin-flop"  response  is  favored  by  some  theoretical  models  [4],  it 
does  not  appear  to  be  directly  responsible  for  the  exchange  biasing 
for  these  samples. 

Future  studies  will  focus  on  the  characteristics  and  origin 
of  the  "spin-flop"  response  of  the  NiO  spin  structure  to  large  mag- 
netic fields.  In  addition,  we  will  further  explore  the  differences 
between  the  antiferromagnetic  domains  for  the  exchange-biased  and 
unbiased  conditions. 

REFERENCES: 

[1]  W.  H.  Meiklejohn  and  C.  P.  Bean,  Phys.  Rev.  102,  1413  (1956);  105,  904 
(1957). 

[2]  A.  P.  Malozemoff,  J.  Appl.  Phys.  63,  3874  (1988). 

[3]  T.  C.  Schulthess  and  W.  H.  Butler.  Phys.  Rev.  Lett.  81,  4516  (1998). 

[4]  N.  C.  Koon,  Phys.  Rev.  Lett.  78,  4865  (1997). 

[5]  J.  A.  Borchers,  R.  W.  Erwin,  S.  D.  Berry,  D.  M.  Lind.  J.  F.  Ankner,  E. 
Lochner,  K.  A.  Shaw  and  D.  Hilton,  Phys.  Rev.  B  51,  8276  (1995). 

[6]  J.  A.  Borchers,  Y.  Ijiri,  D.  M.  Lind,  P.  G.  Ivanov,  R.  W.  Erwin,  S.-H.  Lee 
and  C.  F.  Majkrzak,  J.  Appl.  Phys.  85,  5883  (1999). 


NIST  CENTER  FOR  NEUTRON  RESEARCH 


9 


SPIN  DENSITY  WAVE  ORDER  IN  THE 
SUPERCONDUCTING  STATE  OF  La2Cu04+Y 


n  the  high  temperature  (high-T )  superconductors,  the  multiple 
oles  played  by  the  electrons  continue  to  defy  theoretical  under- 
standing. It  appears  that  all  high-T  superconductors  are  based  on 
structures  with  Cu02  planes,  in  which  the  electrons  on  neighboring 
copper  ions  are  strongly  coupled  magnetically.  In  systems  based 
on  La,Cu04,  antiferromagnetism  is  a  dominant  feature  of  the  phase 
diagram  at  low  doping  levels,  and  conventional  itinerant-electron 
behavior  dominates  in  the  high  doping  regime.  Intermediate  doping 
levels  are  described  by  neither,  but  this  is  where  the  superconduct- 
ing properties  are  optimal.  In  our  experiments,  we  focus  on  the  evo- 
lution of  magnetic  properties  from  the  insulating  antiferromagnet  to 
the  superconductor. 

It  is  becoming  increasingly  apparent  that  incommensurate 
spin  structures  are  universal  to  the  high-T  superconductors. 
Especially  noteworthy  is  the  observation  of  static  incommensurate 
magnetic  ordering  coexisting  with  superconductivity  in 
La, ,  Ndn,Sr  CuO,  and  recently  in  La,  Sr  CuO,  [1,  2].  Further 

1.6-x      0.4     x  4  ■>  2-x     x  4  L  J 

experiments  are  needed  to  understand  the  coexistence  of  the  incom- 
mensurate spin  density  waves  (SDW)  and  superconductivity,  the 
details  of  the  spin  structure,  and  the  influence  of  pinning  potentials. 

A  crystal  of  La,Cu04+y  with  a  superconducting  transition 
at  =  42  K  was  produced  by  doping  pure  La2Cu04  electrochemi- 
cally  with  a  large  quantity  of  excess  oxygen  (y  «  0.12).  The 
superconducting  shielding  signal  measured  after  cooling  in  zero 
field  is  shown  in  Fig.  la.  The  transition  is  very  sharp  with  an  onset 
T  =  42  K.  Our  initial  characterization  suggests  that  the  crystal  is 
a  bulk  superconductor  with  a  hole  concentration  similar  in  density 
and  homogeneity  to  that  of  La2  SrCu04  crystals  with  x  0.15. 

At  temperatures  below  the  superconducting  T ,  we  observe 
elastic  magnetic  scattering  at  a  quartet  of  incommensurate  positions 
centered  around  (100),  which  is  the  Bragg  position  for  the  antifer- 
romagnetism in  the  undoped  insulator.  Surprisingly,  we  find  that  the 
incommensurate  wavevectors  are  not  precisely  along  the  Cu-O-Cu 
bond  direction,  but  are  rotated  by  about  3°. 

In  Fig.  2a  we  show  elastic  scans  along  the  in-plane  h 
direction  through  an  incommensurate  position  for  various  tempera- 
tures using  5  meV  neutrons.  Below  42  K,  the  observed  peaks 
are  extremely  sharp  and  are  resolution-limited,  while  above  42  K 
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FIGURE  1.  (a)  Magnetic  susceptibility  measured  after  cooling  in  zero 
field,  (b)  Peak  intensity  of  the  incommensurate  elastic  scattering  as  a 
function  of  temperature.  The  measurement  was  performed  with  two  different 
neutron  energies  of  13.7  meV  and  5  meV.  The  solid  line  denotes  the  BCS 
superconducting  order  parameter  squared,  with  a  T  of     41  K. 

the  peaks  disappear.  The  solid  lines  in  the  figure  are  Gaussians 
convolved  with  the  instrumental  resolution  which  indicate  that  the 
in-plane  static  magnetic  order  is  correlated  over  distances  larger 
than  400  A.  From  this,  we  conclude  that  static  long-range  magnetic 
order  exists  in  the  superconducting  state  of  La,Cu04+y.  Also,  the 
SDW  order  is  not  specific  to  a  tetragonal  crystal  structure  as  previ- 
ously believed  since  this  crystal  is  orthorhombic;  it  is  a  more 
general  phenomenon. 

We  then  investigated  how  the  static  spin  arrangement  is 
correlated  between  CuO,  planes.  The  /-dependence  of  the  incom- 
mensurate scattering  is  shown  in  Fig.  2b.  The  intensity  modulation 
of  both  the  (10/)-  and  (01/)-centered  scattering  is  reminiscent  of  the 
spin  structure  of  the  undoped  parent  compound  La,Cu04.  The  solid 
lines  in  both  panels  represent  fits  to  Gaussian  lineshapes  convolved 
with  the  instrumental  resolution,  assuming  a  model  for  the  stacking 
arrangement  and  spin  direction  identical  to  that  of  pure  La,Cu04. 
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Here,  the  only  free  parameters  are  the  width  and  a  single  overall 
intensity  scale  factor.  The  agreement  is  clearly  satisfactory,  with  the 
fit  indicating  that  the  spins  are  correlated  across  »  3  CuO,  planes. 
We  conclude  that  the  stacking  arrangement  of  the  magnetically 
ordered  planes  in  our  La,Cu04+y  sample  follows  that  of  undoped 
insulating  La,Cu04,  even  though  the  magnetic  order  in  the  CuO, 
planes  is  incommensurate.  This  is  the  first  direct  evidence  that  the 
magnetism  of  the  doped  superconductor  mimics  the  magnetism  in 
the  undoped  insulator  in  such  a  specific  way. 

We  show  in  Fig.  lb  the  peak  intensity  of  the  elastic  signal 
as  a  function  of  temperature  measured  using  both  13.7  meV  and 
5  meV  neutrons.  The  fact  that  one  obtains  identical  results  for  the 
temperature  dependences  of  the  intensities  with  these  two  different 
neutron  energies  and,  concomitantly,  energy  resolutions  indicates 
that  the  scattering  is  truly  elastic.  The  intensity  of  the  elastic  scat- 
tering turns  on  at  approximately  the  same  transition  temperature  as 
superconductivity.  Noting  that  the  intensity  of  the  magnetic  scatter- 
ing is  proportional  to  the  square  of  the  magnetic  order  parameter, 
we  plot  the  square  of  the  BCS  order  parameter  curve  over  the  data 
using  a  T  of    41  K.  The  agreement  indicates  that  the  magnetism 
exhibits  mean  field  behavior  just  like  conventional  superconductiv- 
ity. This  is  very  surprising  given  the  two  dimensionality  of  the 
ordered  magnetism.  The  size  of  the  ordered  moment  is  0.15  [xB, 
which  is  25  %  of  the  ordered  moment  in  pure  La,Cu04.  Our 
results  argue  against  an  itinerant  electron  description  of  the  incom- 
mensurate magnetism  since  it  is  difficult  to  see  how  a  delocalized 
model  can  support  interplanar  spin  correlations  and  choose  the 
same  preferred  spin  direction  as  in  insulating  La2Cu04.  It  appears 
that  the  spins  are  localized  and  ordered  in  this  high  temperature 
superconductor. 
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FIGURE  2.  (a)  Scans  along  the  in-plane  H  direction  over  one  of  the  incom- 
mensurate SDW  peaks  for  various  temperatures,  (b)  The  top  panel  shows  the 
/-dependence  of  the  SDW  scattering  centered  about  the  (10/)  position.  The 
bottom  panel  shows  scattering  centered  about  the  (01/)  position. 
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LOCAL  SPIN  RESONANCE  AND 
SPIN-PEIERLS-LIKE  PHASE 
TRANSITION  IN  ZnCr  O 


I  ost  magnets  order  when  the  thermal  energy  drops  below 
I  a  microscopic  energy  scale  for  magnetic  interactions.  The 
topology  of  certain  lattices  can,  however,  reduce  the  energy  gain 
associated  with  long  range  ordering  thus  favoring  more  entropic 
phases.  Figure  1  shows  one  such  lattice,  the  tetravalent  site  in 
the  pyrochlore  structure,  that  forms  a  network  of  corner-sharing 
tetrahedra.  As  shown  in  the  figure,  antiferromagnetically  (AFM) 
interacting  spins  on  the  lattice  cannot  satisfy  all  their  exchange 
interactions  simultaneously.  This  phenomenon,  called  geometrical 
frustration,  can  lead  to  macroscopic  classical  ground  state  degen- 
eracy and  offers  the  possibility  of  qualitatively  new  states  of  matter. 
Theoretical  studies  have  in  fact  shown  that  spins  with  nearest  neigh- 
bor antiferromagnetic  Heisenberg  interactions  on  the  pyrochlore  lat- 
tice do  not  have  a  long  range  ordered  phase  at  all. 

Pyrochlore  magnets  studied  experimentally  so  far  exhibit  a 
coninuous  phase  transition  at  a  finite  temperature,  Tf,  into  a  glassy 
phase  with  static  short  range  correlations  [1].  Spinel  antiferromag- 
nets,  AB,04,  in  which  the  octahedral  B  site  forms  the  same  mag- 
netic lattice  as  in  the  pyrochlore  structure,  however  behave  quite 
differently.  For  instance,  ZnCr,04  exhibits  a  first  order  phase  transi- 
tion to  a  long-range  ordered  Neel  phase  at  T  =  12.5  K  ,  much  less 
than  the  Curie-Weiss  temperature  10   I  =  393  K.  We  have  explored 
this  ordered  phase  and  the  corresponding  phase  transition  through 
inelastic  neutron  scattering  [2], 

Figure  2  provides  an  overview  of  our  neutron  scattering 
results  as  a  color  image  of  /  (Q,  w)  at  three  temperatures.  For 
T  >  T ,  Figures  2a  and  2b  show  a  constant-Q  ridge  centered  at 
Q  *  1.5  A1  and  extending  beyond  10  meV.  The  ridge  indicates 
quantum  critical  fluctuations  of  small  AFM  clusters,  most  likely 
antiferromagnetically  correlated  tetrahedra,  and  closely  resembles 
those  obtained  in  similar  experiments  on  other  frustrated  AFM. 
For  T  <  T  ,  however,  the  low  energy  spectral  weight  concentrates 
into  a  sharp  constant-energy  mode  centered  at  fid)  =  4.5  meV 
=  IJI  »  kBT..  The  wave  vector  dependence  of  this  resonance 
intensity  reveals  that  it  is  an  excitation  among  antiferromagnetically 
correlated  nearest  neighbor  spins.  Though  they  can  not  be  seen 
in  Fig.  2c,  there  are  in  fact  magnetic  Bragg  peaks  in  the  elastic 
scattering  channel  (see  Fig.  3b),  which  provide  evidence  for  long- 
range  order  for  T  <  T .  It  is  unusual  that  excitations  of  such  local- 
ized character  exist  in  a  long-range  ordered  phase.  The  resonance 


Frustrated 


FIGURE  1 .  A  network  of  corner-sharing  tetrahedra.  When  two  spins  in  a 
tetrahedron  are  aligned  antiparallel  to  each  other  then  the  third  spin 
can  not  satisfy  its  antiferromagnetic  interaction  with  the  other  two  spins 
simultaneously. 

indicates  the  presence  of  weakly  interacting  spin  clusters  within  the 
ordered  phase,  which  is  a  key  feature  of  geometrically  frustrated 
magnets. 

Theoretical  work  indicates  that  magnetic  order  cannot  devel- 
op in  an  isotropic  spin  system  with  nearest  neighbor  antiferromag- 
netic Heisenberg  interactions  on  the  pyrochlore  lattice.  It  is  natural 
to  ask  what  deviation  from  this  model  causes  order  to  develop 
in  ZnCr204?  To  answer  the  question  we  probed  the  temperature 
dependence  of  static  and  dynamic  features  of  this  system  in  the 
vicinity  of  the  phase  transition.  Figure  3  shows  that  long  range 
antiferromagnetic  order  (blue  squares  in  frame  (b))  and  the  local 
spin  resonance  (frame  (a))  appear  simultaneously  in  a  spectacular 
first  order  transition.  It  also  shows  that  magnetic  ordering  is  accom- 
panied by  a  cubic  to  tetragonal  lattice  distortion  (red  circles  in 
frame  (b)).  The  lattice  distortion  plays  a  crucial  role  in  relieving 
frustration  and  allowing  long-range  order  to  develop.  It  is  well 
known  that  exchange  interaction  between  Cr3+  ions  whose  oxygen 
coordination  octahedra  share  an  edge  are  strongly  dependent  on  the 
oxygen  bond  angles  and  hence  the  metal  ion  spacing.  As  a  conse- 
quence the  tetragonal  strain  ea  >  0  and  ec  <  0  yields  weaker  AFM 
interactions  between  spins  occupying  the  same  basal  plane  and 
stronger  AFM  interactions  between  all  other  spin  pairs.  This  reduces 
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the  energy  of  a  particular  long  range  ordered  spin  configuration  with 
respect  to  the  low  energy  degenerate  manifold  thus  allowing  the 
system  to  achieve  long  range  order.  The  overall  picture  that  emerges 
is  that  of  two  distinct  phases  in  competition:  a  cubic  cooperative 
paramagnet  and  a  tetragonal  long-range  ordered  antiferromagnet. 
Though  the  spin  Hamiltonian  has  a  lower  expectation  value  in  the 
latter  phase,  the  lattice  energy  is  greater  and  the  entropy  is  lower 
in  the  tetragonal  phase.  The  phase  transition  occurs  when  the  free 
energy  of  the  tetragonal  low  entropy  phase  drops  below  that  of  the 
disordered  cubic  paramagnet. 

There  are  strong  analogies  between  the  phase  transition  in 
ZnCr,04  and  the  spin-Peierls  (SP)  transition.  In  both  cases  the 
high  T  phase  is  nearly  quantum  critical  and  can  lower  its  energy 
through  a  lattice  distortion.  In  both  cases  the  transition  occurs  from 
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FIGURE  3.  (a)  Image  of  inelastic  neutron  scattering  for  Q  =  1.5  A1,  (b) 

T- dependence  of  magnetic  Bragg  scattering  from  a  powder,  <rm  (blue  squares), 

and  of  lattice  strains  measured  n  a  single  crystal  (red  circles). 

a  strongly  correlated  paramagnet:  T  «  0CW,  and  in  both  cases  low 
energy  spectral  weight  is  moved  into  a  finite  energy  peak. 

There  are  also  important  differences  between  the  two  transi- 
tions. The  low  T  phases  are  qualitatively  different,  the  transition  in 
ZnCr,04  is  a  first  order  one,  while  the  SP  transition  is  second  order, 
and  the  change  in  entropy  at  T  plays  an  important  role  in  ZnCr204, 
but  not  in  a  SP  transition.  The  central  idea  that  finite  lattice  rigidity 
can  preclude  a  spin  liquid  at  T  =  0  however  does  carry  over  and 
should  be  relevant  for  any  frustrated  magnet  when  other  symmetry 
breaking  interactions  are  sufficiently  weak. 
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FIGURE  2.  Contour  maps  of  the  magnetic  neutron  scattering  intensity  at 
temperatures  spanning  the  phase  transition  at  T=  12.5(5)  K.  The  data  were 
taken  by  utilizing  a  flat  analyzer  and  two-dimensional  position-sensitive 
detector  at  the  SPINS  spectrometer. 
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COMPOUND  REFRACTIVE  OPTICS  IMPROVE 
RESOLUTION  OF  30  m  SANS  INSTRUMENT 


■  n  the  most  favorable  cases,  cold  neutrons  can  be  deflected 
I  through  an  angle  of  a  degree  or  two  by  grazing  incidence  reflec- 
tion, but  by  only  an  arc  second  or  two  by  refraction.  Hence  grazing 
incidence  reflection  optics  has  long  been  considered  the  most  prom- 
ising means  for  focusing  neutrons  for  applications  such  as  small- 
angle  neutron  scattering  (SANS).  Numerous  attempts  over  more 
than  30  years  to  produce  highly  reflective  surfaces  for  neutrons 
have  been  vitiated,  however,  by  SANS  from  the  mirror  surfaces 
themselves,  which  blurs  the  focus.  The  best  mirrors  produced  thus 
far  are  only  marginally  better  for  SANS  than  pinhole  collimation, 
i.e.,  circular  apertures  separated  by  long  distances. 

Scientists  at  Bell  Laboratories  recently  took  a  fresh  look  at 
this  problem  and  proposed  that  multiple  refraction  from  relatively 
high  index,  low  absorbing  material  could  be  superior  to  reflection 
optics  or  conventional  pinhole  collimation  for  SANS.  Initial  mea- 
surements [1]  at  Ris0  National  Laboratory,  Denmark,  demonstrated 
the  proposed  focusing  effect,  but  did  not  make  quantitative  com- 
parisons with  reflection  optics  or  pinhole  collimation  for  application 
in  SANS  instruments.  Measurements  made  recently  at  the  NCNR 
in  collaboration  with  the  Bell  Labs  scientists  [2]  have  addressed 
these  issues  and  have  demonstrated  and  quantified  the  significant 


FIGURE  1 .  Now  installed  in  the  pre-sample  flight  path  of  the  NG-7  30  m 
SANS  instrument  are  two  sets  of  MgF2  biconcave  lenses  that  can  be  inserted 
into  the  beam  under  computer  control.  The  28-lens  array  in  the  foreground 
focuses  8.44  A  neutrons  at  a  distance  of  15  m  from  the  lenses,  and  the 
6-lens  set  focuses  18  A  neutrons  at  the  same  distance.  Each  lens  is  25  mm 
in  diameter,  has  a  radius  of  curvature  of  25  mm,  and  is  1  mm  thick  in  the 
center. 

14 


Source  Sample 

aperture  aperture  Detector 


FIGURE  2.  Upper  panel,  conventional  SANS  pinhole  collimation.  The  source  and 
sample  apertures,  A1p  and  A2p,  respectively,  determine  the  shape  and  extent 
of  the  beam  profile,  l(x),  at  the  detector  plane.  Lower  panel,  focusing  lens 
geometry.  Ideally,  the  source  aperture,  Au,  alone  determines  the  beam  profile. 

improvement  in  resolution  that  can  be  achieved  with  compound 
refractive  optics. 

Our  tests  were  made  with  the  same  set  of  cylindrical  bicon- 
cave MgF,  (magnesium  fluoride)  lenses  used  in  the  Ris0  study.  Up 
to  30  lenses  were  placed  end-to-end  near  the  sample  position  of  the 
30  m  SANS  instrument  to  focus  neutrons,  emanating  from  a  circular 
source  aperture  15  m  upstream,  onto  the  plane  of  instrument's  two- 
dimensional  detector.  Figure  1  shows  an  array  of  28  lenses  for 
focusing  8.44  A  neutrons  at  a  distance  of  15  m  from  the  sample, 
next  to  a  set  of  6  lenses  for  focusing  18  A  neutrons  at  the  same 
distance. 

For  this  geometry,  as  depicted  in  Fig.  2,  the  lenses  ideally 
produce  a  1:1  image  of  the  source  aperture  at  the  detector  inde- 
pendent of  the  size  of  the  sample.  Since  the  scattering  signal  is 
proportional  to  sample  size,  the  lens  system  can,  in  principle,  be 
used  to  improve  resolution  more  efficiently,  by  reducing  the  size  of 
the  source  aperture,  than  is  possible  with  pinhole  collimation  where 
both  the  source  aperture  and  sample  size  must  be  reduced  propor- 
tionally to  improve  angular  resolution.  Aberrations  and  small-angle 
scattering  by  the  lenses  could,  however,  blur  the  image  to  such 
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FIGURE  3.  Red  dots  are  measured  points  along  the  vertical  profile  of  the  image 
formed  by  the  28  lens  array  shown  in  Fig.  1  of  a  9.5  mm  diameter  source  of 
neutrons  15  m  upstream  from  the  lenses.  The  blue  dots  are  a  Monte  Carlo 
calculation  of  the  profile  that  includes  the  effects  of  spherical  and  chromatic 
aberrations  as  well  as  the  broadening  caused  by  gravity.  The  shoulder  in  the 
measured  profile  at  Y  =  0  is  due  to  a  residual  fast  neutron  component  in 
the  beam. 

a  degree  that  any  advantage  over  pinhole  collimation,  which  does 
produce  a  sharply  defined  beam  spot  at  the  detector,  would 
be  lost. 

To  accurately  measure  the  intensity  profile  produced  by  the 
lenses,  a  dysprosium  foil  was  positioned  at  the  focal  plane  and 
exposed  to  the  focused  beam  for  approximately  2  h.  The  activated 
foil  was  then  placed  in  contact  with  a  high  resolution  image  plate 
which  stored  the  image  produced  by  the  emitted  gamma  rays  with  a 
spatial  resolution  of  better  than  0. 1  mm.  A  typical  profile  obtained 
from  reading  out  the  image  plate  is  shown  in  Fig.  3.  Also  plotted  in 
the  figure  is  a  Monte  Carlo  calculation  of  the  profile  that  includes 
the  effects  of  spherical  and  chromatic  aberrations  as  well  as  the 
broadening  caused  by  gravity.  The  measured  profile  agrees  with  the 
simulation  down  to  intensity  levels  of  10'3  of  the  peak  intensity  and 
has  an  overall  signal-to-background  ratio  in  the  wings  approaching 
105,  which  is  highly  satisfactory  for  most  SANS  measurements. 

The  practical  benefit  provided  by  the  lenses  is  demonstrated 
by  the  SANS  data  from  voids  in  a  single  crystal  (2.5  cm  in  diameter 
and  1  cm  thick)  of  fast-neutron-irradiated  aluminum  shown  in 
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FIGURE  4.  Small-angle  scattering  from  voids  in  a  single  crystal  (2.5  cm  in 
diameter  and  1  cm  thick)  of  fast-neutron-irradiated  aluminum.  The  measure- 
ments were  made  under  equivalent  resolution  conditions  (i.e.  nearly  identical 
beam  spot  size  at  the  detector)  using  both  simple  pinhole  collimation  and  the 
28  biconcave  lens  array  shown  in  Fig.  1.  The  integrated  gain  in  intensity  due 
to  the  lenses  is  approximately  26. 


Fig.  4.  The  measurements  were  made  under  equivalent  resolution 
conditions  (i.e.,  nearly  identical  beam  spot  size  at  the  detector) 
using  both  simple  pinhole  collimation  and  the  28  biconcave  lens 
array  shown  in  Fig.  1.  The  first  data  point  unaffected  by  the  beam 
stop  in  both  data  sets  occurs  at  Q  «=  0.001  A"1,  but  the  scattered 
intensity  at  the  detector  (counts/cm2/s)  is  more  than  10  times  higher 
by  using  the  lenses  to  illuminate  a  much  larger  area  of  the  sample. 

The  focusing  lenses  shown  in  Fig.  1  are  now  installed  for 
routine  use  in  the  NCNR's  30  m  SANS  instrument  on  guide  NG-7. 
Further  testing  is  planned  to  understand,  and  hopefully  eliminate, 
the  sources  of  parasitic  scattering  that  contribute  to  the  tails  of  the 
beam  profile  seen  in  Fig.  3,  prior  to  installing  a  lens  system  in  the 
NIST/NSF  30  m  SANS  instrument  on  guide  NG-3. 
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POLYMER  BRUSH 

RESPONSE  TO  SOLVENT  FLOW 


The  presence  of  polymer  chains  grafted  or  adsorbed  onto  a 
surface  can  dramatically  alter  the  forces  that  affect  interactions 
between  surfaces.  The  equilibrium  properties  of  such  polymer  brush 
systems  have  been  studied  for  the  past  two  decades,  yielding  gener- 
al agreement  between  theory  and  experiment.  Conversely,  the  non- 
equilibrium  properties  of  polymer  brushes  are  still  under  intense 
theoretical  and  experimental  investigation.  Of  particular  interest 
is  the  response  of  a  brush  to  the  frictional  forces  imposed  by 
solvent  flow.  The  behavior  of  polymer  brushes  subjected  to  flow 
has  important  technological  implications  for  the  rheology  of  col- 
loidal dispersions  stabilized  by  polymer  layers,  for  the  lubrication 
properties  of  polymer  coated  interfaces,  for  biocompatibility  of 
medical  implant  devices,  and  for  permeation  flow  through  polymer- 
containing  porous  media  [1]. 

The  height  of  a  polymer  brush  is  determined  by  the  equi- 
librium conformation  of  the  tethered  chains,  which  depends  on 
both  the  grafting  density  and  quality  of  solvent.  The  basic  physics 
governing  the  static  behavior  of  a  polymer  brush  result  from  a 
competition  between  two  opposing  tendencies:  1)  elastic  contrac- 
tion, as  the  chains  attempts  to  maximize  their  entropy  by  adopting 
random  walk  configurations,  and  2)  monomer-monomer  interac- 
tions, such  as  polymer-polymer  repulsions,  and  polymer-solvent 
wetting  [2,  3]. 

Polymer  chain  stretching  in  densely  grafted  brushes  has  been 
studied  by  many  different  techniques  including  surface  forces  appa- 
ratus [4],  neutron  reflectivity  [5],  and  small  angle  neutron  scattering 
(SANS)  [6].  In  general,  there  is  good  agreement  with  results  from 
experiment,  simulations  and  analytical  calculations  [3]. 

Oscillatory  shear  measurements  performed  with  a  surface 
forces  apparatus  suggest  that  the  normal  forces  between  a  pair  of 
brush  surfaces  are  altered  when  sheared.  However,  these  measure- 
ments do  not  give  the  actual  brush  profile  either  with  or  without 
shear.  Effective  hydrodynamic  thickness  measurements  of  polymer 
brushes  under  shear  indicate  a  thickening  of  the  brush;  whereas  neu- 
tron reflectivity  experiments  on  adsorbed  PS-PEO  block  copolymer 
brushes  on  a  silica  surface  show  no  effect  of  shear  on  the  brush 
density  profile  in  good  solvent,  and  a  slight  increase  in  poor  solvent. 
These  earlier  reflectivity  measurements  were  limited  to  shear  rates 
of     10,000  s"1  [7]  since  the  adsorbed  PS-PEO  block  copolymer 


FIGURE  1.  Schematic  illustrating  the  effect  of  shear  on  grafted  brushes 
predicted  by  Miao  et  al.  [8]  Such  an  effect  would  be  consistent  with  the 
present  data. 
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FIGURE  2.  Shear  cell.  Arrows  denote  solvent  flow  direction.  The  polystyrene 
brush  is  grafted  onto  the  Si  crystal. 

brushes  have  a  tendency  to  come  off  the  surface  at  higher  shear 
rates. 

Predictions  from  theoretical  calculations  of  brush  profiles 
under  shear  span  the  gamut  of  possibilities,  ranging  from  brush 
thickening  to  brush  compression,  including  no  effect  of  shear  flow 
on  the  density  profile  [3].  Miao  et  al.  [8]  predict  that  the  response  of 
a  brush  to  the  solvent  shear  flow  is  displayed  as  chain  tilting  toward 
and  chain  stretching  along  the  direction  of  flow.  However,  the 
overall  conformational  properties  such  as  brush  thickness  remain 
essentially  unaffected  (Fig.  1 ). 

We  have  performed  neutron  reflectivity  measurements  on  a 
chemically  grafted  polymer  in  both  good  and  poor  solvents  at  shear 
rates  over  an  order  of  magnitude  greater  than  previously  reported. 
Our  neutron  reflectivity  experiments  measure  the  segment  density 
profile  of  the  polymer  brushes  under  shear  in  an  experimental 
cell  similar  to  the  one  used  by  Baker  et  al.  [9],  (Fig.  2).  We 
use  deuterated  polystyrene  (d-PS),  83  kg/mol,  with  a  trichloro- 
silane  end  group  to  bind  the  d-PS  brush  chemically  onto  a  single 
crystal  Si  surface  [5].  We  used  a  good  solvent,  toluene,  and  a  poor 
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FIGURE  3.  Representative  shear  data  for  deuterated  polystyrene  brush 
under  shear  in  cyclohexane  (a),  and  toluene  (b).  In  both  cases,  black 
circles  represent  shear  of  0  s\  whereas  red  squares  represent  data  taken 
at  30,000  s1. 

solvent,  cyclohexane,  as  the  solvent  media.  The  dry  brush  height 
was  measured  by  X-ray  reflectivity  to  be  17.5  nm.  Without  shear 
the  brush  extends  to  31  nm  in  cyclohexane  and  75  nm  in  toluene. 
We  measured  the  brush  profile  at  several  shear  rates,  up  to  130,000 
s1,  yet  we  see  no  effect  of  shear  on  the  brush  density  profiles  in 
either  solvent  (Fig.  3).  No  desorption  of  the  polymer  brush  was 
ever  observed.  In  fact,  the  neutron  reflectivity  profiles  at  0  s1  and 
130,000  s"1  look  identical,  indicating  less  than  a  2  %  to  3  %  change 
in  the  brush  density  profile. 

We  have  also  been  able  to  establish  that  the  slight  shear 
induced  swelling  reported  by  Baker  et  al.  [7]  for  a  poor  solvent 
(cyclohexane),  was  probably  due  to  frictional  heating  of  the  solvent. 
In  high  shear  fields,  heat  generated  from  friction  between  the 
solvent  and  interior  surfaces  of  the  apparatus  does  not  readily 
dissipate,  causing  a  =  2  °C  to  3  °C  rise  in  the  temperature  of 
the  shear  cell.  We  were  able  to  demonstrate  that  the  brush  height 
in  cyclohexane  is  unaffected  by  shear  when  the  cell  temperature 
is  carefully  controlled.  An  elevated  cyclohexane  temperature  swells 
the  brush  as  the  solvent  quality  improves.  Naturally,  these  effects 
were  not  observed  in  toluene. 

Thus,  our  neutron  reflectivity  data  represent  the  first  compre- 
hensive measurements  of  shear  effects  on  the  density  profile  of  a 
grafted  polymer  brush  into  regimes  that  are  predicted  by  some  to 
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display  an  effect.  We  cannot  at  this  time  determine  if  the  brush 
responds  as  predicted  by  Miao  et  al.  [8],  or  if  there  is  insufficient 
solvent  penetration  into  the  brush  to  exert  enough  force  on  the 
chains  to  induce  conformational  changes.  Future  measurements  will 
distinguish  between  these  cases. 
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MACROMOLECULAR  CONFORMATION 
IN  ULTRATHIN  POLYMER  FILMS 


Since  thin  polymeric  films  are  ubiquitous  in  technological  appli- 
cations such  as  paints,  lubricants,  and  adhesives,  a  critical 
characterization  of  their  thermophysical  properties  is  essential.  A 
central  premise  in  the  development  of  theories  for  predicting  the 
properties  of  polymer  melts  in  confined  geometries  is  that  chains 
maintain  their  unperturbed  Gaussian  conformation,  which  they 
adopt  in  the  bulk,  in  the  direction  parallel  to  the  surfaces  under  all 
conditions  [1,2].  These  assumptions,  which  form  the  foundations  of 
the  important  field  of  polymer  thin  films,  have  been  questioned  on 
the  basis  of  indirect  experimental  findings  [3,4].  We  have  utilized 
the  power  of  small  angle  neutron  scattering,  especially  the  high 
neutron  flux  at  NCNR,  to  unequivocally  characterize  the  chain 
structure  and  conformation  in  ultrathin  polymer  films,  and  thus  have 
resolved  this  important  fundamental  question. 

The  experimental  measurements  of  chain  conformations  and 
system  thermodynamics  in  thin  films  have  remained  elusive  due 
to  the  small  amounts  of  sample  material  involved.  To  illustrate 
this  point,  a  thin  film  of  10  nm  incurs  a  decrease  in  signal  by 

lxl04from  a  typical  bulk  polymer  sample.  In  this  case,  the  noise 
is  comparable  to  the  signal,  complicating  the  experiments.  Prior  to 
upgrades  of  the  cold  neutron  source  at  the  NCNR,  data  collection 
times  were  prohibitive.  Improvements  in  sample  preparation,  which 
are  discussed  in  detail  elsewhere  [5,6],  have  allowed  us  to  measure 
molecular  size  and  conformation  of  an  isotopically  labeled  blend 
of  polystyrene  (25  wt%  d-PS/75  wt%  h-PS)  in  films  as  thin  as 
12  nm.  Two  blends  of  nominally  matched  molecular  weight,  Mn, 
of  270,000  and  650,000,  respectively,  were  utilized.  These  were 
labeled  270k  and  650k,  respectively.  Solutions  of  the  blends  were 
spin  cast  on  silicon  substrates  (Semiconductor  Processing)  and 
annealed  at  120  °C  (T  «  105  °C). 


The  scattered  intensity  for  a  15  nm  thick  film  of  the  270k 
blend  is  compared  to  an  analogous  bulk  sample  in  Fig.  1.  It  is  clear 
from  the  figure  that,  on  a  unit  volume  basis,  the  thin  film  scattering 
is  higher  than  that  of  the  bulk.  We  postulated  that  this  difference  is 
attributable  to  the  scattering  from  the  imperfections  at  both  the  air 
and  the  substrate  interfaces,  which  is  driven  by  the  relatively  high 
neutron  contrast  at  these  boundaries.  To  evaluate  this  component, 
films  of  pure  d-PS  were  spin  cast  under  identical  conditions  and 
their  scattering  measured. 

The  pure  d-PS  film  data  were  fit  with  a  simple  Debye-Bueche 
form  factor  to  obtain  parameters  for  a  roughness  term.  The  blend 
film  data  were  then  fit  by  scaling  this  roughness  term  and  adding  a 
component  obtained  from  the  Random  Phase  Approximation  (RPA) 
model.  In  the  fitting,  two  RPA  model  parameters  were  also  varied: 
the  blend  chain  radius  of  gyration  (Rg)  and  the  Flory  interaction 
parameter.  The  combined  model  is  illustrated  as  the  solid  line  in 
Fig.  1  along  with  the  roughness  term  (long  dashed  line)  and  RPA 
term  (short  dashed  line).  The  fact  that  the  RPA  term  is  very  close 
to  the  data  for  the  bulk  blend  illustrates  that  the  film  and  bulk 
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FIGURE  1 .  Plots  of  l(Q)  as  a  function  of  0  for  dPS/hPS  blends  of  Mn  = 
270k.  Filled  symbols,  bulk  sample;  open  symbols,  data  for  a  film  of  thickness 
D  =  18  nm.  The  fit  to  these  data  (solid  line)  was  obtained  by  utilizing  both 
a  roughness  term  (long  dashed  line)  as  well  as  the  standard  RPA  form  (short 
dashed  line). 
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samples  have  nearly  the  same  Rc.  Using  this  combined  model,  Rg 
was  determined  for  films  ranging  in  thickness  over  two  decades 
(0.5  Ro  <  D  <  50  RJ.  The  molecular  size  was  found  to  be  indepen- 
dent of  film  thickness  (Fig.  2).  Since  the  scattering  vector  is  primar- 
ily in  the  surface  plane,  this  conclusion  is  consistent  with  theoretical 
assumptions  and  suggests  that  chain  conformation  in  the  direction 
parallel  to  the  surfaces  are  unaffected  by  confinement. 
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FIGURE  2.  Plots  of  ratios  of  R  and  1(0),  derived  from  the  RPA  component 
(short  dashed  line  in  Fig.  1)  of  the  fits  to  the  blend  thin  film  data,  to  their 
corresponding  bulk  values.  The  ratios  are  plotted  versus  the  ratio  of  the  film 
thickness  to  the  bulk  Rg.  Data  are  displayed  for  blends  with  Mn  =  270k  (filled 
circles)  and  Mn  =  650k  (open  circles). 


Our  results  clearly  show  that,  in  the  thinnest  films,  the 
volume  pervaded  by  a  coil  is  decreased  as  compared  to  the  bulk. 
This  is  because  the  Rg  in  the  direction  parallel  to  the  surfaces  is 
unaffected,  while  the  corresponding  quantity  in  the  third  direction  is 
strongly  reduced.  In  conjunction  with  other  studies,  which  indicate 
a  thickness-independent  density  in  ultrathin  polymer  films,  these 
conclusions  indicate  decreased  intermolecular  entanglement  in  thin 
polymer  films.  Since  entanglement  density  directly  affects  the 
dynamic  properties  of  polymeric  systems,  we  contend  that  unusual 
thin  film  properties,  such  as  the  anomalous  thickness  dependence  of 
diffusion  coefficients  and  glass  transition  temperatures,  are  caused 
by  this  reduced  entanglement  density  near  a  surface. 

With  thin  films  as  a  model  system,  and  continuing  increas- 
es in  cold  neutron  flux,  SANS  at  the  NCNR  is  now  an  appropriate 
tool  to  study  a  host  of  problems  involving  interfacial  structure, 
finite  size  phase  behavior,  and  nano-patterning  in  systems  as  far 
ranging  as  engineering  thermoplastics  to  biological  systems.  These, 
and  related  problems,  are  the  focus  of  investigation  in  our  research 
groups. 
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PROBING  THE  LOCATION  OF  THE 
TERMINAL  GROUPS  OF  A  DENDRIMER 


Dendrimers  represent  a  new  class  of  macromolecules  developed 
in  recent  years.  Typically,  a  dendrimer  structure  has  a  tri-  or 
tetrafunctional  core,  which  is  surrounded  by  several  "generations" 
of  stepwise  added  trifunctional  monomers,  leaving  the  last  genera- 
tion with  a  large  number  of  terminal  units  as  shown  in  Fig.  1 . 
The  molecular  weight  doubles  with  each  generation,  leading  to 
high  molecular  weights,  and  causing  the  dendrimer  to  become  very 
compact  and  crowded. 

Many  of  the  potential  technological  applications  of  den- 
drimers depend  on  their  segment  density  distribution.  Previous  scat- 
tering studies  have  shown  that  dendrimers  have  uniform  interiors 
and  are  quite  spherelike  in  their  shape  [1].  The  location  of  the 
terminal  groups  is  also  of  importance,  since  they  are  usually  differ- 
ent chemically  from  the  rest  of  the  dendrimer.  This  invites  a  number 
of  applications  such  as  the  support  of  catalysts  or  drugs  or  their 
use  as  hyperfunctional  crosslink  sites.  The  accessibility  of  these 
terminal  groups  depends  on  their  location  compared  to  the  other 
dendrimer  units. 


FIGURE  1 .  Dendrimer  structure  with  labeled  terminal  units. 


The  location  of  the  terminal  units  can  be  determined  by  label- 
ing the  last  generation  of  the  dendrimer  with  deuterium  and  using 
contrast  matching  techniques.  Figure  1  shows  the  labeled  groups  in 
red  and  the  rest  of  the  dendrimer  in  blue.  By  choosing  the  proper 
mix  of  h-  and  d-  solvents,  the  interior  of  the  dendrimer  will  be 
matched,  making  only  the  labeled  end  groups  visible  in  small  angle 
neutron  scattering  (SANS). 

A  sixth  generation  polyamidoamine  (PAMAM)  dendrimer 
was  reacted  with  acrylonitrile  (vinyl-d3)  to  give  the  deuterium 
labeling  for  the  SANS.  Ethylene  diamine  was  reacted  with  the 
dendrimer  to  give  a  labeled  seventh  generation  dendrimer.  A  similar 
reaction  was  used  to  make  a  seventh  and  eighth  generation  den- 
drimer without  labeling.  Solutions  of  unlabeled  eighth  generation 
dendrimer  were  made  in  mixtures  of  CH?OH  and  CD^OH  for  deter- 
mination of  the  match  point.  Three  samples  were  analyzed,  an 
unlabeled  dendrimer  in  CD,OH  (high  contrast),  an  unlabeled  den- 
drimer in  the  match  mixture  (dendrimer  matched),  and  the  labeled 
dendrimer  in  the  match  mixture  (interior  matched). 
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FIGURE  2.  SANS  from  G8  dendrimer  in  CD,0H/CH30H  mixtures. 
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SANS  was  performed  at  the  30  m  facilities  at  NIST.  The 
spectrometers  were  operated  at  a  wavelength  of  X  -  6  A,  and  a 
wavelength  spread  of  AA./A.  =  0.15. 

Figure  2  is  a  plot  of  SANS  of  dendrimer  solutions  with  differ- 
ent CD,OH  contents.  The  intensity  is  the  strongest  in  pure  CD,OH, 
weakens  as  CH,OH  is  added  and  increases  again  when  pure  CH,OH 
is  used.  The  coherent  scattering  intensity  varies  as  I  «  (bs  -  bD)2 
where  bD  is  the  contrast  of  the  dendrimer  and  bs  is  the  average 
contrast  of  the  solvent  mixture. 

Figure  3  is  a  plot  of  the  square  root  of  the  scattered  intensity 
versus  solvent  composition  with  the  values  to  the  right  made  nega- 
tive so  that  a  straight  line  can  be  put  through  all  of  the  data  points. 
The  zero  intersection  is  at  a  mass  fraction  of  60.5  CH,OH  which 
was  the  composition  used  in  the  matching  experiments. 

Figure  4  is  a  plot  of  the  SANS  of  the  three  G7  samples. 
The  circles  give  the  scattering  from  the  high  contrast  sample,  show- 
ing strong  scattering  typical  of  large  spherical  dendrimers.  The 
diamonds  show  the  SANS  of  the  same  dendrimer,  but  under  match 
conditions.  This  sample  has  no  measurable  coherent  scattering  sig- 


nal, demonstrating  that  the  match  conditions  have  been  achieved. 
The  labeled  dendrimer  SANS  is  given  by  the  squares.  The  scattering 
is  weak  because  only  the  labeled  terminal  groups  scatter. 

A  Guinier  analysis  of  the  scattering  of  the  high  contrast 
sample  gives  the  radius  of  gyration  (Rg)  of  the  whole  dendrimer,  and 
the  labeled  -  contrast  match  sample  gives  the  Ro  of  only  the  terminal 
groups.  The  Roof  the  whole  dendrimer  is  (34.2  ±  0.2)  A,  while  the 
Rg  of  the  terminal  groups  is  (39.3  ±  1.0)  A. 

The  terminal  groups  of  a  seventh  generation  PAMAM  den- 
drimer are  15  %  larger  than  the  average  of  all  of  the  units. 
Therefore,  the  terminal  units  of  a  dendrimer  are  concentrated  in  the 
outer  shell  of  a  dendrimer. 
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FIGURE  3.  Location  of  the  match  point. 


FIGURE  4.  SANS  of  labeled  and  unlabeled  G7  dendrimer  in  match  and 
high  contrast  solvents. 
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SELF-ASSEMBLIES  OF 
MEMBRANE  ACTIVE  PEPTIDES 


I  n  recent  years  we  have  made  progress  by  neutron  diffraction 
■  in  a  major  structural  challenge,  namely  in  detecting  and  analyz- 
ing the  structures  of  peptide  assemblies  in  fluid  membranes  [1,2]. 
These  experiments  were  performed  with  membranes  in  the  form 
of  oriented  multilayers.  Originally  the  samples  were  investigated  in 
full  hydration  so  that  the  physical  properties  of  the  lipid  bilayers 
were  close  to  those  at  physiological  conditions.  However,  it  was 
soon  realized  that  new  phenomena  involving  peptide-lipid  inter- 
actions occur  when  the  sample  hydration  is  varied.  In  general, 
in  full  hydration,  the  peptide  organization  in  each  membrane 
is  uncorrected  to  the  neighboring  membranes.  As  the  hydration 
level  decreases,  the  peptides  become  correlated  between  bilayers, 
even  though  the  membranes  are  still  fluid.  In  many  cases,  further 
dehydration  strengthens  the  correlation  such  that  the  peptide  orga- 
nization in  the  multilayers  crystallizes  [2].  The  crystallization 
provides  the  possibility  for  high-resolution  diffractional  studies. 
Investigations  along  this  line  might  also  lead  to  useful  information 
for  crystallization  of  membrane  proteins. 

Antimicrobial  peptides  are  inducible  innate  host  defense  mol- 
ecules found  in  all  multicellular  organisms,  including  humans  and 
plants.  These  peptides  have  a  folded  size  comparable  to  the  mem- 
brane thickness.  All  evidence  indicates  that  antimicrobial  peptides 
act  by  permeabilizing  the  cell  membranes  of  microorganisms.  But 
the  molecular  mechanisms  of  their  actions  are  still  not  clear.  We 
have  found  that  all  peptides,  when  they  are  bound  to  lipid  bilayers, 
exhibit  two  distinct  oriented  circular  dichroism  spectra,  one  at  low 
peptide-to-lipid  ratios  (P/L)  and  another  at  high  P/L.  This  indicates 
that  each  peptide  has  two  different  physical  states  of  binding  to  a 
membrane. 

The  transition  from  the  low  to  the  high  P/L  spectrum  occurs 
over  a  narrow  range  of  P/L  as  if  there  is  a  threshold  concentration, 
called  P/L*.  At  concentrations  below  P/L*,  the  peptides  are  embed- 
ded in  the  headgroup  region,  as  suggested  by  the  peptide  orientation 
and  the  membrane  thinning  effect.  At  concentrations  above  P/L*, 
neutron  in-plane  scattering  showed  that  the  peptides  form  pores 
in  the  membranes,  while  no  pores  were  detected  below  P/L*.  The 
detection  was  achieved  by  exploiting  the  sensitivity  of  neutrons 
to  D,0,  which  had  replaced  the  water  in  the  membrane  pores. 
As  an  example,  Fig.  1  shows  neutron  in-plane  scattering  taken  at 
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FIGURE  1.  In-piane  SANS  from  alamethicin  in  DLPC  lipid  bilayers  at  a  high 
peptide-to-lipid  ratio  where  the  peptides  self-assemble  to  form  channels 
through  the  lipid  bilayer  as  depicted  in  Fig.  2.  In  the  upper  panel,  the  channels 
were  filled  with  either  D20  (+),  giving  strong  scattering  contrast,  or  H20  (o).  In 
the  lower  panel  the  lipid  was  deuterated,  providing  stronger  scattering  contrast 
with  H20  in  the  channels  (+)  compared  with  D20  (o).  These  data  were  taken 
atANL. 


Argonne  National  Laboratory  (ANL)  of  alamethicin  in  protonated 
lipid  bilayers  (Fig.  la)  and  in  deuterated  lipid  bilayers  (Fig.  lb), 
with  D,0  or  H,0  filling  the  pore  channels.  The  peak  in  these  data 
is  due  to  the  fairly  regular  pore  spacing  and  is  most  pronounced 
when  there  is  a  strong  contrast,  either  between  D20  and  protonated 
lipid,  or  between  H,0  and  deuterated  lipid.  All  four  sets  of  data  are 
consistent  with  the  model  shown  in  Fig.  2  once  the  differences  in 
contrast  are  taken  into  account  (solid  curves).  Thus  we  concluded 
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FIGURE  2.  Model  for  the  channels  formed  by  octamers  of  alamethicin  in  DLPC 
bi layers  that  is  consistent  with  the  in-plane  SANS  data  shown  in  Fig.  1. 
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FIGURE  3.  Examples  of  SANS  patterns  from  channel-forming  peptides  in  lipid 
bilayers  at  various  temperatures  and  stages  of  dehydration.  The  upper  left 
panel  is  for  magainin  pores  in  fully  hydrated  fluid  bilayers.  As  this  and  other 
peptide/lipid  systems  studied  are  dehydrated,  the  pores  in  adjacent  bilayers 
become  correlated  and  eventually  crystallize. 


that  alamethicin  in  DLPC  bilayers  forms  octameric  pores  in  the 
barrel-stave  fashion. 

Interestingly,  the  barrel-stave  model  is  not  the  only  possible 
pore  formation.  We  have  detected  another  type,  called  toroidal 
pores,  in  which  the  lipid  monolayer  bends  continuously  from  one 
leaflet  to  another  like  the  inside  of  a  torus  [1].  However,  while 
the  evidence  for  the  pores  is  clear  by  the  detection  of  the  water 
(D,0)  channels  through  the  lipid  bilayers,  the  evidence  for  the  pore 
structures  is  indirect.  Thus  the  discovery  of  the  crystalline  phases 
is  an  important  new  development  for  the  field  of  antimicrobial 
peptides. 

We  developed  a  method  of  off-plane  scattering  [2]  to  record 
the  diffraction  pattern  on  a  SANS  instrument  that  includes  both  the 
in-plane  and  out-of-plane  momentum  components.  Figure  3  exhibits 
some  typical  diffraction  patterns  as  recorded  on  the  NG-3  30  m 
SANS  instrument's  detector  by  this  method.  The  top  left  panel 
shows  the  diffraction  pattern  of  magainin  pores  in  fully  hydrated 
fluid  bilayers.  When  the  sample  was  slightly  dehydrated,  the  pattern 
changed  to  the  top  middle  panel.  Our  analysis  [2]  showed  that  the 
positions  of  the  magainin  pores  in  each  bilayer  become  correlated 
with  the  pores  in  adjacent  bilayers,  even  though  the  bilayers  are  still 
in  the  fluid  phase.  The  cause  of  this  correlation  was  hypothesized  to 
be  due  to  the  hydration  force.  Upon  further  dehydration  or  cooling, 
the  pore  arrangement  crystallized  into  a  lattice  (the  left  panel  of 
the  middle  row)  having  ABCABC  stacking  of  hexagonally  ordered 
planes. 

The  other  crstalline  patterns  shown  in  Fig.  3  were  from 
magainin  and  other  peptides  in  various  lipids  at  different  hydration 
levels  and  temperatures.  These  crystalline  phases  offer  a  new  way 
of  studying  peptide-lipid  interactions  that  will  help  us  to  understand 
the  molecular  mechanisms  of  membrane  active  peptides  in  cell 
biology. 
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SHEAR  ORIENTATION  OF 

VISCOELASTIC  POLYMER-CLAY  SOLUTIONS 


Shear-induced  structural  changes  in  complex  fluids  of  anisotro- 
pic species  are  very  general  phenomena,  occurring  in  polymer 
solutions,  liquid  crystalline  materials  and  block  copolymer  melts. 
The  purpose  of  our  work  is  to  investigate  the  influence  of  shear 
on  the  structure  of  a  highly  viscoelastic,  aqueous  clay-polymer  solu- 
tion. Many  structural  models  have  been  proposed  for  such  solutions 
[1-3], but  little  is  definitively  known  about  mesoscopic  properties  or 
shear  behavior.  This  information  is  important  in  the  production  of 
nanocomposite  materials  [4]. 

In  our  work,  we  use  small-angle  neutron  scattering  (SANS)  to 
study  a  solution  of  the  synthetic  hectorite  type  clay,  Laponite  LRD 
(Laporte  Industries  Ltd.),  and  poly(efhylene-oxide)  (PEO)  (Mw=  106 
g/mol).  The  results  reported  here  are  for  a  highly  viscoelastic  solu- 
tion containing  a  mass  fraction  of  3  %  LRD  and  2  %  PEO  at  room 
temperature.  The  clay  particles  produce  transparent  dispersions  of 
disk  shaped  particles  ca.  300  A  in  diameter  and  ca.  10  A  thick  [5,6]. 
The  pH  and  ionic  strength  of  the  solutions  were  controlled  by  the 
addition  of  NaOH  and  NaCl,  respectively. 

Figure  1  shows  the  shear  rate  dependence  of  the  birefringence 
of  the  clay-polymer  solution.  A  distinct  minimum  in  the  birefrin- 
gence is  observed  at  a  critical  shear  rate  of  approximately  40  s1. 
The  source  of  the  shear  dependence  of  the  birefringence  is  due  to 
the  alignment  of  the  clay  particles  and  the  PEO.  Previous  work 
demonstrated  that  the  sign  of  the  birefringence  of  the  clay  particles 
oriented  along  a  flow  field  is  negative,  therefore  at  low  shear  rates, 
the  orientation  of  the  clay  dominates  the  birefringence.  Above  the 
critical  shear  rate,  the  birefringence  due  to  the  orientation  of  the 
polymer  chains  dominates. 

A  double  logarithmic  plot  of  viscosity,  v\,  versus  shear  rate 
shows  that  the  solution  is  shear-thinning  over  the  entire  range 
according  to  a  power  law  with  exponent  m  =  -0.65.  No  signature  of 
the  critical  shear  rate  is  observed  in  the  viscosity  behavior. 

The  SANS  shear  cell  utilized  has  been  described  previ- 
ously [2].  It  consists  of  a  cylinder  that  rotates  within  an  outer 
cylinder  with  the  sample  in  the  gap  between  them.  The  instrument 
was  configured  in  both  "radial"  (incident  beam  parallel  to  the 
shear  gradient  along  the  cylinder  diameter)  and  "tangential"  (inciT 
dent  beam  passing  between  the  cylinders,  parallel  to  the  flow  direc- 
tion) geometries.  Using  9  A  wavelength  neutrons  gives  a  Q  range 
between  0.0027  A"1  and  0.0199  A1.  The  primary  contrast  in  the 
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FIGURE  1 .  Optical  birefringence  as  a  function  of  shear  rate.  The  arrow 
indicates  the  shear  rate  where  the  minimum  in  the  birefringence  occurs. 

SANS  experiment  used  to  detect  the  orientation  of  the  clay  platelets 
and  polymer  chains  under  shear  is  between  D20  and  the  other 
solution  components. 

The  results  obtained  from  the  polymer-clay  solutions  in  the 
"radial"  and  "tangential"  beam  configurations  are  summarized  in 
Fig.  2.  At  low  shear  rates,  a  diffuse  isotropic  ring  of  SANS  intensity 
is  observed  (Fig.  2a).  The  diffuse  ring  corresponds  to  an  average 
spacing  between  platelets  of  800  A  to  1100  A.  With  increasing  shear 
rate,  the  ring  becomes  more  diffuse  (Fig.  2b)  and  an  anisotropic 
streak  develops  parallel  to  the  vorticity  axis  of  the  flow  field  (the 
cylinder  axis).  If  we  neglect  the  main  reflected  beam  which  appears 
as  a  background  streak  in  the  gradient  direction  for  tangential  beam 
measurements  (Fig.  2d),  the  anisotropic  streak  becomes  the  domi- 
nant feature  in  both  scattering  geometries  with  increasing  shear  rate 
After  cessation  of  shear,  the  streaks  relaxed  to  an  isotropic  state  in  | 
less  than  2  min. 

To  account  for  the  SANS  and  birefringence  results,  our  cur 
rent  understanding  is  that  the  polymer  chains  are  in  a  dynamic 
adsorption/desorption  equilibrium  with  the  clay  particles  to  form 
a  network.  The  peak  position  in  the  quiescent  scattering  pattern 
in  Fig.  2a  is  an  indication  of  the  mesh  size  of  this  network 
(<=»  1000  A).  A  2  %  solution  of  only  PEO,  at  the  same  pH,  polymer 
and  salt  concentration  showed  no  anisotropic  SANS  scattering  at 
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shear  rates  up  to  100  s1.  Similarly  a  3  %  aqueous  clay  solution 
shows  no  evidence  of  an  anisotropic  SANS  pattern.  Therefore,  we 
can  conclude  that  the  anisotropic  SANS  pattern  observed  in  the 
clay-polymer  solutions  is  due  to  this  coupling  between  clay  platelets 
and  polymer,  allowing  a  higher  orientation  than  either  single  com- 
ponent in  solution  can  produce.  From  the  birefringence  data,  the 
clay  particles  orient  at  low  shear  rates,  while  strong  orientation  of 
the  PEO  does  not  occur  until  the  critical  shear  rate  is  exceeded. 
Since  the  clay  platelets  and  the  PEO  chains  are  of  comparable  size 
(both  about  300  A),  the  lack  of  internal  flexibility  of  the  rigid  clay 
particles  makes  them  much  easier  to  align  than  the  flexible  polymer 
chains. 

According  to  SANS  patterns  from  both  beam  configurations 
(Fig.  2)  the  shear  flow  results  in  an  alignment  of  clay  platelets 
orienting  with  their  surface  normals  in  the  vorticity  direction.  One 
would  expect  the  surface  normals  to  orient  along  the  gradient  direc- 


a)  radial  beam  b)  radial  beam 
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FIGURE  2.  SANS  patterns  of  the  clay-polymer  solutions  as  a  function  of  shear 
rate  in  the  radial  geometry  (a-c)  and  tangential  geometry  (d)  at  shear  rates  of 
a)  0.5  s"1 ,  b)  20  s 1 ,  c)  90  s 1 ,  d)  90  s 1 . 


tion  of  the  flow  field,  however,  the  type  of  orientation  observed 
in  these  clay-polymer  solutions  is  also  observed  in  some  liquid 
crystalline  lamellar  phases,  block  copolymer  solutions,  and  melts. 
The  critical  shear  rate  is  the  shear  rate  at  which  the  rate  of  chain 
desorption  is  slower  than  the  terminal  relaxation  time  of  the  chain, 
hence  chain  extension  is  observed  in  the  birefringence. 

On  cessation  of  shear,  the  stress  on  the  network  decays 
almost  immediately,  and  the  recovery  of  the  isotropic  structure  is 
controlled  by  the  relaxation  of  the  stretched  chains.  As  the  chains 
retract,  the  coupling  of  the  chains  to  the  clay  allows  the  platelets 
to  randomize  in  orientation  in  the  local  viscous  environment.  The 
recovery  from  anisotropy  is  much  faster  than  expected  from  simple 
Brownian  motion  of  only  the  clay  particles  in  a  medium  of  the 
same  viscosity  as  the  clay-polymer  solution  exhibited  macroscopi- 
cally,  and  is  indicative  of  the  dynamic  coupling  of  the  polymer 
chains  to  the  clay.  Future  work  will  compare  the  relative  rates  of 
the  relaxation  in  the  PEO  and  clay  with  the  cooperative  adsorption/ 
desorption  kinetics  which  occur  during  deformation. 
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UNIQUE  DETERMINATION  OF 
BIOMIMETIC  MEMBRANE  PROFILES  BY 
NEUTRON  REFLECTIVITY 


ew  biomimetic  membrane  materials,  of  fundamental  impor- 
tance in  understanding  such  key  biological  processes  as 


molecular  recognition,  conformational  changes,  and  molecular  self- 
assembly,  can  be  characterized  using  neutron  reflectometry.  In  par- 
ticular, scattering  length  density  (SLD)  depth  profiles  along  the 
normal  to  the  surface  of  a  model  biological  bilayer,  which  mimics 
the  structure  and  function  of  a  genuine  cell  membrane,  can  be 
deduced  from  specular  neutron  reflectivity  data  collected  as  a  func- 
tion of  wavevector  transfer  Q.  Specifically,  this  depth  profile  can  be 
obtained  by  numerically  fitting  a  computed  to  a  measured  reflectiv- 
ity. The  profile  generating  the  best  fitting  reflectivity  curve  can 
then  be  compared  to  cross-sectional  slices  of  the  film's  chemical 
composition  predicted,  for  example,  by  molecular  dynamics  simula- 
tions [1].  However,  the  uniqueness  of  a  profile  obtained  by  conven- 
tional analysis  of  the  film's  reflectivity  alone  cannot  be  established 
definitively  without  additional  information.  In  practice,  significantly 
different  SLD  profiles  have  been  shown  to  yield  calculated  reflectiv- 
ity curves  with  essentially  equivalent  goodness-of-fit  to  measured 
data  [2],  as  illustrated  in  Fig.  1. 


The  existence  of  multiple  solutions,  only  one  of  which  can  be 
physical,  is  especially  problematic  in  cases  where  a  key  additional 
piece  of  structural  or  compositional  information  is  lacking  as  can 
happen  in  the  investigation  of  these  biological  membrane  systems. 

Why  this  inherent  uncertainty?  The  neutron  specular  reflection 
amplitude  for  a  model  SLD  can  be  computed  exactly  from  first 
principles;  the  square  of  its  modulus  gives  the  measurable  reflectiv- 
ity. It  is  firmly  established,  however,  that  the  complex  amplitude 
is  necessary  and  sufficient  for  a  unique  solution  of  the  inverse 
problem,  that  of  recovering  the  SLD  from  reflection  measurements. 
Unambiguous  inversion  requires  both  the  magnitude  and  phase  of 
reflection.  Once  these  are  known,  practical  methods  [3]  exist  for 
extracting  the  desired  SLD. 

In  fact,  considerable  efforts  were  made  about  a  quarter  century 
ago  to  solve  the  analogous  "phase  problem"  in  X-ray  crystallography 
using  known  constraints  on  the  scattering  electron  density  [4]  and  by 
the  technique  of  isomorphic  substitution  [5].  Variations  of  the  latter 
approach  have  been  applied  to  reflectivity,  using  a  known  reference 
layer  in  a  composite  film  in  place  of  atomic  substitutions.  These 


2- 


0.1        .  0.2 
Q(A"1) 


FIGURE  1 .  Family  of  scattering  length  density  profiles  obtained  by  model- 
independent  fitting  of  the  reflectivity  data  in  the  inset.  The  profile  represented 
by  the  blue  dashed  line  is  unphysical  for  this  Ti/TiO  film  system  yet  generates 
a  reflectivity  curve  that  fits  the  data  with  essentially  equivalent  goodness-of-fit 
(all  the  reflectivity  curves  corresponding  to  the  SLD's  shown  are  plotted  in  the 
inset  but  are  practically  indistinguishable  from  one  another). 


FIGURE  2.  Reflectivity  curves  for  the  thin  film  system  depicted  schematically  in 
the  inset,  one  for  a  Si  fronting  (red  triangles),  the  other  for  Al203  (black  circles). 
The  curve  in  the  lower  part  of  the  figure  (blue  squares)  is  the  real  part  of  the 
complex  reflection  amplitude  for  the  films  obtained  from  the  reflectivity  curves 
by  the  method  described  in  the  text. 
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solution  methods,  however,  were  tied  to  the  Born  approximation, 
which  generally  is  valid  in  crystal  structure  determination  but  which 
fails  catastrophic  ally  at  low  Q  (low  glancing  angles)  in  reflection 
from  slab-shaped  samples  such  as  thin  films.  Exact  inversion 
requires  accurate  knowledge  of  the  reflection  amplitude  over  the 
entire  Q-range,  especially  at  low  Q. 

In  this  decade  the  reflection  phase  problem  has  been  exactly 
solved  using  a  protocol  of  three  reflectivity  measurements  on  com- 
posite films  consisting  of  the  film  of  interest  in  intimate  contact  with 
each  of  three  known  reference  layers  [6,  7].  Subsequently,  variations 
using  only  two  measurements  have  been  shown  to  partially  solve 
the  phase  problem,  an  additional  procedure  being  required  to  choose 
between  two  solution  branches,  only  one  of  which  is  physical  [8, 
9].  In  the  past  year  [10],  an  exact  solution  has  been  found  for  a 
two  measurement  strategy  in  which  the  film  surround,  either  the 
fronting  (incident)  or  backing  (transmitting)  medium,  is  varied.  This 
new  approach  is  simpler  to  apply  than  reference  layer  methods 
and  is  adaptable  to  many  experiments.  Surround  variation  neutron 


FIGURE  3.  SLD  profile  (red  line)  resulting  from  a  direct  inversion  of  the  Re  r  of 
Fig.  2  compared  with  that  predicted  by  a  molecular  dynamics  simulation  (white 
line)  as  discussed  in  the  text.  The  headgroup  for  the  Self-Assembled-Monolayer 
(SAM)  at  the  Au  surface  in  the  actual  experiment  was  ethylene  oxide  and  was 
not  included  in  the  simulation  but,  rather,  modelled  separately  as  part  of  the 
Au.  Also,  the  Cr-Au  layer  used  in  the  model  happened  to  be  20  A  thicker  than 
that  actually  measured  in  the  experiment. 


reflectometry  has  been  successfully  applied  to  the  challenging  type 
of  biological  membrane  depth  profiling  described  earlier. 

In  Fig.  2  are  plotted  a  pair  of  neutron  reflectivity  curves 
measured  for  the  layered  film  structure  schematically  depicted  in 
the  upper  right  inset,  one  with  Si  and  the  other  with  A1203  as  the 
fronting  medium.  The  lower  part  of  Fig.  2  shows  the  real  part  of 
the  complex  reflection  amplitude  for  the  multilayer  as  extracted  from 
the  reflectivity  data,  according  to  the  method  described  above,  and 
which  was  subsequently  used  to  perform  the  inversion  to  obtain 
the  SLD  shown  in  Fig.  3.  For  comparison,  the  SLD  predicted  by 
a  molecular  dynamics  simulation  is  also  shown  in  Fig.  3,  in  a 
slightly  distorted  version,  corresponding  to  a  truncated  reflectivity 
data  set,  which  indicates  the  spatial  resolution  of  an  SLD  obtainable 
in  practice.  This  latter  SLD  was  obtained  by  inversion  of  the  reflec- 
tion amplitude  computed  for  the  exact  model  SLD,  but  using  values 
only  up  to  the  same  maximum  Q  value  (0.3  A1)  over  which 
the  actual  reflectivity  data  sets  were  collected.  Overall,  agreement 
between  the  experimentally  determined  profile  and  the  theoretical 
prediction  is  remarkable,  essentially  limited  only  by  the  Q-range  of 
the  measurement.  Surround  variation  neutron  reflectivity  thus  makes 
it  possible  to  measure  complicated  thin  film  structures  without  the 
ambiguity  associated  with  curve  fitting.  The  veridical  SLD  profile  is 
obtained  directly  by  a  first  principles  inversion. 
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MULTI-TECHNIQUE  STUDIES  OF 
ULTRATHIN  SiQ2  FILMS 


f  I 


Cun"ent  gate  dielectrics  in  silicon  based  devices  are  only  a 
few  nm  thick.  Optical  techniques  such  as  ellipsometry  are 
used  to  monitor  film  thicknesses  and  optical  properties  in  produc- 
tion. However,  for  the  current  integrated  circuit  (IC)  generation 
the  accuracy  of  ellipsometry  degrades  because  parameters  such 
as  thickness  and  index  of  refraction  (which  reflects  the  composi- 
tion) become  strongly  correlated.  Thus,  it  is  difficult  to  unambigu- 
ously determine  these  parameters  simultaneously,  and  the  accuracy 
of  ellipsometry  would  benefit  from  an  independent  calibration.  In 
reflectometry  techniques,  on  the  other  hand,  these  parameters  are 
nearly  decoupled.  The  thickness  of  a  layer  is  approximately  inverse- 
ly proportional  to  the  oscillation  period  of  the  reflected  intensity, 
whereas  the  differences  in  scattering  length  density  SLD  (also  an 
indicator  of  composition)  between  the  layers  is  related  to  the  ampli- 
tude of  the  oscillations.  Neutron  reflectometry  (NR)  is  better  suited 
than  X-ray  reflectometry  (XR)  for  the  study  of  the  Si02/Si  system 
because  there  is  a  relatively  large  contrast  (or  difference  in  SLD) 
between  the  scattering  length  densities  of  the  two  materials:  65  %, 
vs.  7.6  %  for  X-rays. 

Consider  as  an  example  a  sample  with  a  nominally  10  nm 
thick  thermal  oxide  film  on  silicon.  This  moderate  thickness  was 
chosen  to  increase  our  confidence  in  the  results  of  the  various 


characterization  methods,  while  remaining  thin  enough  that  the 
results  are  relevant  to  film  of  technological  interest.  Figure  1  shows 
spectroscopic  ellipsometry  (SE)  data  and  corresponding  best  fits  for 
the  sample  with  surface  contamination  and  after  an  organic  cleans- 
ing. Nominally,  the  only  change  is  a  decrease  in  the  thickness  of  the 
contamination  layer  [1]. 

In  XR  data  (Fig.  2)  two  oscillation  periods  are  observed 
for  the  contaminated  sample.  The  high  frequency  oscillation  corre- 
sponds to  the  SiO,  film,  whereas  the  low  frequency  modulation  is 
due  to  the  thinner  contamination  layer  (which  is  not  present  after 
cleaning,  indicating  removal  of  the  contamination.) 

The  NR  measurements  (Fig.  3)  were  done  in  a  vacuum  to 
reduce  the  air  scattering  background.  This  allowed  us  to  achieve 
a  very  large  range  in  reflectivity,  over  108,  which  is  among  the 
best  examples  in  NR  measurements  to  date.  A  slightly  thinner 
contamination  layer  in  NR  is  consistent  with  the  fact  that  the  XR 
was  done  in  air,  during  which  the  contamination  was  growing.  This 
was  confirmed  by  changes  in  XR  scans  immediately  following  those 
in  Fig.  2. 

The  average  of  the  5  measurements  of  the  SiO,  film  thickness 
was  10.27  ±0.13  nm.  The  excellent  agreement  among  the  results 
for  the  three  different  techniques  increases  our  confidence  in  the 
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FIGURE  1.  Comparison  of  spectroscopic  ellipsometry  experimental  data,  x\' 
and  A,  to  the  fits  (solid  line)  for  the  clean  and  surface  contaminated  sample. 
The  inset  shows  <  e,>  and  <s2>  as  a  function  of  depth  determined  by  the  fit. 


FIGURE  2.  X-ray  reflectivity  and  best  fits  for  the  clean  and  surface 
contaminated  sample.  The  inset  shows  the  scattering  length  density  profile 
determined  by  the  fits. 
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FIGURE  3.  Neutron  Reflectivity  data  and  best  fit  for  a  sample  with  surface 
contamination.  The  inset  is  the  scattering  length  density  profile  determined 
by  the  fit. 

parameters  extracted  via  these  models.  Thus  the  XR  and  particularly 
the  NR  corroborate  the  correct  analysis  required  in  SE  (which  is  the 
technique  most  practical  in  monitoring  production). 

To  further  investigate  the  applicability  of  these  techniques  to 
thinner  SiO,  layers,  we  simulated  the  SE  and  the  XR  and  NR  curves 
for  6  nm,  4  nm,  and  2  nm  thick  layers  of  SiO,.  For  a  realistic  and 
consistent  set  of  roughness  parameters  in  the  models,  we  used  the 
average  values  obtained  from  the  actual  measurements  previously 
discussed.  The  SE  simulation,  Fig.  4a,  shows  distinct  differences 
in  both  the  magnitude  and  shape  of  A  among  the  three  thicknesses 
shown.  In  models  of  XR,  shown  in  Fig.  4b,  only  very  weak  oscilla- 
tions are  seen  for  even  the  thickest  of  the  SiO,  layers  because  of 
the  low  contrast  between  SiO,  and  Si.  However,  in  NR,  for  SiO, 
layers  as  thin  as  2  nm  strong  oscillations  are  clearly  seen  above 
the  10"8  lower  limit,  demonstrated  in  Fig.  3.  Therefore  both  NR  and 
SE  are  well  suited  for  the  study  of  SiO,  films  as  thin  as  2  nm. 
Encouraged  by  these  models,  we  obtained  NR  data  for  a  thinner, 
2.4  nm,  sample.  While  these  data  are  not  yet  fit  to  a  model  curve, 
we  note  that  both  the  reflected  intensities  and  oscillation  amplitude 
are  similar  to  those  of  the  2.0  nm  model,  indicating  similar  interface 
widths. 
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FIGURE  4.  Model  calculations  for  thin  S1O2  Alms  of  thickness  indicated  in  the 
figure,  on  Si.  a)  Spectroscopic  ellipsometry,  and  b)  unless  otherwise  noted 
solid  lines  are  neutron  reflectometry  models.  For  clarity  the  6  nm  (4  nm) 
films  are  shifted  up  by  2(1)  orders  of  magnitude.  The  data  points  are  for  a 
==  2.4  nm  film. 

We  have  shown  that  three  different  techniques  can  offer  com- 
plimentary information  on  the  structure  of  thin  SiO,  films  on  Si. 
All  offer  a  significant  degree  of  sub-monolayer  thickness  sensitivity, 
although  in  NR  there  is  a  much  higher  contrast  between  SiO,  and 
Si  than  in  XR. 
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CERTIFICATION  OF  AN  ION-IMPLANTED 

ARSENIC  IN  SILICON 

STANDARD  REFERENCE  MATERIAL 


The  secondary-ion  mass  spectrometry  (SIMS)  community  in 
the  United  States  recently  undertook  a  round-robin  study  to 
calibrate  the  implanted  dose  of  arsenic  in  silicon  by  consensus. 
"Variations  in  dose  determination  among  laboratories  were  as  high  as 
30  %,  reflecting  primarily  the  errors  of  the  respective  in-house  stan- 
dards. By  contrast,  in  an  international  round-robin  exercise  spon- 
sored by  the  International  Standards  Organization  (ISO)  the  U.S. 
laboratories  achieved  a  relative  standard  deviation  of  4  %  for  nine 
independent  determinations  of  the  boron  content  in  an  unknown 
boron-doped  silicon  sample.  This  level  of  agreement  was  only  pos- 
sible because  all  the  laboratories  used  SRM  2137  Boron  Implant  in 
Silicon  as  a  common  reference  material.  These  results  demonstrate 
the  need  for  a  common  arsenic  reference  material  to  improve  inter- 
laboratory  reproducibility.  Furthermore,  SEMATECH  (a  consortium 
of  semiconductor  manufacturers)  recently  listed  SRM  implants  of 
phosphorus  and  arsenic  in  silicon  as  high  priority  industrial  needs. 
Consequently,  a  Standard  Reference  Material  (SRM  2134  Ion- 
Implanted  Arsenic  in  Silicon)  was  produced  at  NIST,  using  a  wafer 
from  the  SIMS  intercomparison.  The  material  for  this  SRM  was 
provided  by  a  major  ion  implanter  manufacturer  who  supplied  three 
200  mm  diameter  wafers  that  had  been  implanted  with  arsenic  at 
an  energy  of  100  keV  (in  the  same  batch).  One  of  the  wafers  was 
diced  and  distributed  to  12  participating  laboratories  for  the  arsenic 
round-robin  study  described  above.  Each  of  the  remaining  wafers 
could  provide  221  SRM  units,  enough  for  an  estimated  10  year  sup- 
ply. One  of  the  remaining  wafers  was  therefore  diced  into  1  cm  x  1 
cm  pieces  with  a  wafer  saw  for  use  as  SRM  2134.  Because  of  the 
specificity  and  matrix  independence  of  instrumental  neutron  activa- 
tion analysis  (INAA),  this  technique  was  chosen  as  the  primary 
method  for  certification  of  the  arsenic  implanted  dose. 

Ten  analytical  specimens  of  silicon  were  selected  from 
among  the  22 1  pieces  cut  from  the  wafer  to  determine  arsenic 
content  and  test  for  homogeneity.  Two  blank  chips  of  silicon  were 
analyzed  in  parallel.  Pre-irradiation  handling  was  done  on  a  class 
100  clean  bench.  In  a  Teflon  TFE  beaker,  the  specimens  and  blanks 
were  rinsed  together  in  ultrapure  (UP)  water,  rinsed  with  95  % 
ethanol,  cleaned  ultrasonically  for  one  minute  in  ethanol,  and  rinsed 
twice  with  UP  water.  The  samples  and  blanks  were  then  agitated  for 
one  minute  in  UP  water,  blotted  on  Whatman  41  filter  paper,  and 
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FIGURE  1 .  Region  of  the  gamma-ray  spectrum  for  analysis  of  As  in  SRM  2134. 
Limits  of  peak  and  baseline  areas  are  indicated  by  vertical  lines. 


air-dried.  Each  was  weighed  and  heat-sealed  in  a  polyethylene  bag 
that  had  been  cleaned  in  high  purity  nitric  acid. 

A  dilute  solution  standard  was  prepared  by  gravimetrically 
diluting  a  10  mg/g  solution  of  As^O.,  in  dilute  ammonia  to  about  2 
mg/kg.  Another  solution  was  prepared  similarly  from  SRM  3103a 
Arsenic  Spectrometric  Solution.  Aliquants  of  both  solutions  were 
weighed  from  a  polyethylene  pipette  onto  1-cm  disks  of  acid- 
leached  Whatman  41  filter  paper  and  dried  to  produce  standards 
with  about  200  ng  arsenic  each.  Two  standards  from  each  solution 
were  irradiated  with  the  samples.  A  blank  filter  paper  was  shown  in 
a  separate  experiment  to  contain  no  significant  arsenic  («  0.01  % 
of  the  amount  in  the  standards). 

Samples,  standards,  and  blanks  were  stacked  in  the  irra- 
diation container  in  the  following  sequence:  standard,  blank,  five 
samples,  two  standards,  five  samples,  blank,  and  standard.  The 
stack  was  irradiated  for  four  hours  in  the  RT- 1  pneumatic  tube  of 
the  NIST  research  reactor,  at  a  thermal  fluence  rate  of  7.7  x  10'3 
cnr2s'.  The  container  was  inverted  halfway  through  the  irradiation 
to  equalize  to  the  first  order  for  the  axial  fluence  rate  gradient  along 
the  container. 

After  irradiation  the  silicon  samples  were  cleaned  ultra- 
sonically for  one  minute  in  de-ionized  water,  blotted  dry,  and 
heat-sealed  in  clean  conventional  polyethylene  bags  for  counting. 
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Gamma-ray  activity  was  assayed  at  10  cm  from  a  germanium  detec- 
tor with  a  1.71  keV  resolution  and  40  %  relative  efficiency  at  1333 
keV.  Each  sample  and  standard  was  counted  at  least  twice,  after 
the  2.62  h  31Si  matrix  activity  had  substantially  decayed.  The  559 
+  563  keV  doublet  of  76As  in  each  spectrum  was  integrated  with 
the  SUM  computer  program  [1].  See  Fig.  1.  No  evidence  indicating 
significant  heterogeneity  among  samples  could  be  seen  when  the 
observed  sample-to-sample  precision  was  compared  to  what  was 
expected  from  counting  statistics  combined  with  other  sources  of 
analytical  variability. 

NIST  currently  certifies  elemental  concentrations  in  SRMs 
using  one  of  three  modes:  (1)  a  primary  method  at  NIST  with 
confirmation  by  other  method(s);  (2)  two  independent  critically 
evaluated  methods  at  NIST,  and  (3)  one  method  at  NIST  and  differ- 
ent methods  by  outside  collaborating  laboratories  [2]. 

Table  1.  Individual  Uncertainty  Components  for  Determination  of 
Arsenic  in  SRM  2134. 


Source  of  Uncertainty  Individual  Uncert.  (1s)  % 


Basis  Measure  -  Sample  Area  /  Micrometer  Accuracy 

0.012 

Basis  Measure  -  Sample  Area  /  Precision  (n  =  26)1 

0.0089 

Concentrations  of  Comparators  (Standards) 

0.106 

Mass  Determination  -  Comparators 

0.041 

Isotopic  Variability 

0 

Blank  and  Blank  Correction  (n  =  102)1 

0.013 

Irradiation  Geometry  Differences 

0.039 

Neutron  Self  -  Shielding/Scattering  Differences 

0.024 

Timing 

0 

Irradiation  Interferences 

Negligible 

Measurement  Replication  (n  =  9)1 

0.081 

Counting  Statistics  -  Standards  (n  =  105)1 

0.073 

Counting  Geometry  Differences 

0.009 

Pulse-Pileup 

0.066 

Dead-Time  Effects  (inadequacy  of  LT  extension) 

Negligible 

Decay  Timing  Effects  (Half  life  related) 

0.053 

Gamma-Ray  Self  Shielding 

0.004 

Gamma-Ray  Interferences 

0.00004 

Peak  Integration  Method 

0.033 

Overall 

0.189 

'  '^"indicates  degrees  of  freedom  which  are  listed  for  Type  A  sources  of  uncertainty. 


Certification  using  a  primary  method  is  only  possible  when 
all  potentially  significant  sources  of  uncertainty  have  been  evaluated 
explicitly  for  the  application  of  the  method  and  the  matrix  under 
investigation.  In  addition,  confirmation  of  measurements  by  a  pri- 
mary NIST  method  may  be  accomplished  by  one  or  more  of  the 
following:  determination  of  certified  constituents  in  other  SRM(s) 
or  CRM(s)  of  similar  matrix  and  constituent  concentration  range; 
a  second  NIST  technique  with  appropriate  controls;  or  results  of 
measurements  from  selected  outside  collaborating  laboratories  with 
appropriate  experience. 

In  order  to  certify  the  arsenic  concentration  in  SRM  2134 
using  INAA  as  a  primary  method  it  is  therefore  necessary  to  evalu- 
ate all  significant  sources  of  uncertainty  explicitly.  For  this  set  of 
measurements,  we  considered  sources  of  uncertainty  greater  than 
0.01  %  relative  to  be  significant.  The  results  of  a  complete  evalua- 
tion of  all  sources  of  uncertainty  are  listed  in  Table  1.  This  evalua- 
tion yielded  an  expanded  relative  uncertainty  of  0.38  %  (as  defined 
by  ISO  and  NIST)  and  gives  an  approximate  level  of  confidence 
of  95  %.  The  arsenic  concentration  observed  was  91.20  ng/cnr  ± 
0.35  ng/cm2. 

In  conclusion,  we  have  successfully  applied  INAA  as  a 
primary  method  for  the  certification  of  this  new  SRM.  The  observed 
relative  expanded  uncertainty  of  0.38  %  is  considerably  smaller 
than  the  1  %  value  desired  by  the  semiconductor  industry.  This 
new  SRM  should  greatly  enhance  the  U.S.  semiconductor  industry's 
ability  to  achieve  accurate  and  reproducible  analytical  results  for 
this  key  dopant  in  silicon. 
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THE  ROTATIONAL  DYNAMICS  OF  H2  IN 
POROUS  VYCOR  GLASS 


Porous  materials  play  an  increasingly  important  role  in  the 
exploration  of  fundamental  scientific  issues  related  to  diverse 
technological  applications  including  adhesion,  lubrication,  tribology 
and  the  engineering  of  materials.  Surface  interactions  and  finite  size 
effects  both  play  a  key  role  in  the  qualitative  modification  of  the 
properties  of  the  materials  contained  within  the  porous  host.  Due  to 
its  light  mass  and  weak  electronic  interactions,  molecular  hydrogen 
is  an  ideal  system  to  probe  the  effects  of  surface  interactions.  H2  is 
well  described  as  a  quantum  rigid  rotor  with  discrete  energy  levels 
labeled  by  the  rotational  quantum  number,  J,  and  energies  given  by 
Ej  =  BJ(J  +  1 )  where  B  is  the  rotational  constant  which  is  equal 
to  7.35  meV  for  H,.  These  energy  levels  can  be  studied  by  using 
neutrons  to  stimulate  the  J  =  1  to  J  =  0  transition  of  the  H,  molecule, 
a  process  that  is  normally  doubly  forbidden  due  to  the  quantum  sta- 
tistics obeyed  by  H,.  When  neutron  stimulated  conversion  occurs, 
the  neutron  gains  an  amount  of  energy  equal  to  the  rotational  energy 
of  the  molecule.  This  results  in  an  extremely  clean  signal  since 
all  other  processes  are  frozen  out  at  the  low  temperatures  of  these 
measurements. 

Studies  of  the  rotational  levels  of  H,  adsorbed  in  \y  cor, 
a  commercially  available  porous  glass,  were  carried  out  for  pore 
fillings  ranging  from  0.10  (corresponding  to  less  than  a  single 
monolayer  on  the  pore  surface)  to  0.92  (corresponding  to  nearly  full 
pores).  As  can  be  seen  in  Fig.  1,  at  low  filling  fractions  only  a  single 
broad  peak  at  10  meV  is  present  which  can  be  attributed  to  the 
scattering  from  an  adsorbed  layer  of  H,  strongly  bound  to  the  pore 
surface.  The  appearance  of  scattering  at  14  meV,  corresponding  to 
the  free  molecule  transition,  for  fillings  above  0.45  is  associated 
with  the  appearance  of  "bulk"-like  material  in  the  pore  center. 
The  clear  separation  of  the  scattering  from  molecules  adsorbed 
on  the  walls  and  in  the  pore  center  suggests  that  these  are  quite 
distinct  states  which  we  refer  to  as  the  bound  and  bulk-like  states, 
respectively. 

The  shift  in  the  energy  of  the  rotational  transition  at  low 
fillings  can  be  directly  related  to  the  interaction  of  the  molecules 
with  the  surface  [2].  Figure  2  shows  the  shift  of  the  J  =  1  to 
J  =  0  transition  as  a  function  of  the  rotational  energy  barrier  for 
the  molecule.  From  Fig.  2  one  finds  that  the  observed  rotational 
energy  of  10  meV  corresponds  to  an  orientational  potential  with 
a  barrier  height  of  VB/BH  =  2.7  (VB  =  19.8  meV).  The  magnitude 


FIGURE  1 .  The  top  panel  shows  a  contour  plot  of  the  inelastic  scattering  from 
H2  in  Vycor  at  various  pore  fillings  at  6  K.  The  bottom  panel  is  a  cut  through  the 
contour  plot  at  f  =  0.65  (indicated  by  the  line). 

of  this  potential  is  quite  large,  comparable  to  that  seen  in  alpha 
alumina,  a  catalyst  used  for  ortho-H,  enrichment  [3],  and  on  leached 
glass  [4]. 

Careful  inspection  of  the  data  reveals  that  the  peak  due  to 
the  bound  state  actually  shifts  from  9.1  meV  to  10.1  meV  between 
filling  fractions,  f,  of  0.20  and  0.27.  Thus  for  f  between  0.27  and 
0.45,  the  bound  state  may  be  separated  into  a  layer  in  direct  contact 
with  the  pore  wall  and  a  layer  that  feels  a  weaker  orientational 
potential.  The  scattering  from  these  two  layers  can  be  separately 
determined  assuming  the  intensity  from  the  first  layer  varies  linearly 
with  filling. 

The  scattering  from  the  second  surface  layer  can  be  obtained 
from  the  difference  of  the  monolayer  and  bilayer  spectra.  The  scat- 
tering from  the  second  layer  is  centered  at  E  =  11.3  meV,  compared 
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FIGURE  2.  Shift  of  the  rotational  energy  levels  of  H2  as  a  function  of 
barrier  height. 

to  the  first  layer  which  is  centered  at  9.1  meV,  indicating  that  the 
second  layer  interacts  less  strongly  with  the  wall. 

The  width  of  the  peak  we  associate  with  the  bulk-like  mol- 
ecules is  due  entirely  to  instrumental  resolution.  This  is  consistent 
with  the  view  of  these  molecules  as  being  free  rotors  with  a 
single  well  defined  value  for  B.  On  the  other  hand,  the  width  of 
the  peak  which  we  attribute  to  the  bound  layer  is  much  broader 
than  the  instrumental  resolution.  A  single  well  defined  orientational 
potential  would  yield  a  shift  in  the  peak  location,  as  observed,  but 
no  additional  broadening,  since  the  energy  levels  would  still  be  well 
defined.  Thus,  we  attribute  this  additional  width  to  a  distribution  of 
orientational  potentials  which  can  also  be  directly  extracted  using 
the  model  of  White  and  Lassettre  [2]  (Fig.  2).  After  correcting 
for  instrumental  effects,  the  strength  of  the  scattering  at  a  given 
energy  is  directly  proportional  to  the  number  of  molecules  which 
feel  the  corresponding  orientational  potential.  Thus,  plotting  the 
scattering  as  a  function  of  V/BH,  directly  yields  the  distribution 
of  orientational  potentials  shown  in  Fig.  3.  As  can  be  seen,  the 


— I  

First  Layer 
Second  Layer - 


-4  -2 


4 

V/B 


10  12 


FIGURE  3.  Distribution  of  potential  barriers  felt  by  the  molecules  in  the  first 
(red)  and  second  (blue)  layers. 

orientational  potential  for  the  first  layer  is  broad,  asymmetric,  and 
centered  at  VB=  3.5BH,  whereas  the  potential  of  the  second  layer  is 
relatively  narrow  and  centered  at  VB  =  2.0BH,  indicating  a  weaker 
interaction  with  the  surface. 

The  data  allow  us  to  draw  the  following  picture  of  the  adsorp- 
tion of  molecular  hydrogen  on  a  porous  glass  surface.  The  first 
monolayer  is  bound  tightly  to  the  rough  pore  surface.  The  rough 
surface  prevents  free  rotation  of  the  molecule,  altering  the  rotational 
energy  states.  The  second  layer  perceives  the  surface  roughness 
smoothed  by  the  presence  of  the  first  layer,  and  accordingly,  the 
rotational  transitions  are  affected  to  a  lesser  degree.  Subsequent  to 
the  completion  of  the  second  monolayer,  the  H,  molecules  sense 
no  significant  orientational  interaction  from  the  glass  surface.  This 
is  in  agreement  with  the  picture  presented  by  Katsaros  et  al.  [5], 
in  which  surface  roughness  is  created  by  dangling  bonds  on  the 
pore  surface. 
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FIRST-PRINCIPLES  COMPUTATIONAL  AND 
NEUTRON  SCATTERING  STUDY  OF 
PROTONIC  CONDUCTORS 


Fuel  cells,  which  produce  electricity  via  hydrogen  oxidation  to 
water,  have  emerged  in  the  last  decade  as  one  of  the  key 
technologies  for  meeting  the  world's  energy  needs  well  into  the 
21st  century.  Solid-oxide  fuel  cells  (SOFCs)  are  the  most  promising 
among  the  many  different  types  of  fuel  cells  being  developed. 
However,  the  required  high  operating  temperature  of  an  SOFC 
places  stringent  requirements  on  its  component  materials.  Potential 
candidates  to  replace  the  oxygen-conducting  electrolyte  used  in  cur- 
rent SOFCs  include  the  perovskite-based  high-temperature  protonic 
conductors  (HTPCs),  since  switching  from  anionic  to  protonic  con- 
duction can  lower  the  practical  operating  temperature  from  1300  K 
to  around  1000  K. 

We  undertook  neutron  vibrational  spectroscopy  (NVS)  mea- 
surements of  various  HTPCs  (namely,  SrCeO,  and  SrZr03  aliova- 
lently  doped  with  different  rare  earth  cations)  in  order  to  probe  the 
bonding  potentials  experienced  by  the  residual  protons  incorporated 
as  OH"  species  via  steaming.  At  high  temperatures,  these  protons 
migrate  from  site  to  site,  giving  rise  to  protonic  conduction.  The 
NVS  measurements  revealed  the  existence  of  dopant-related  pertur- 
bations to  the  bending-mode  energies  of  OH  in  SrCe„„sM„„ M  O,  , 

°  °  0.95     0.05    x  3-6 

(M  =  Sc,  Ho,  and  Nd),  confirming  the  trapping  effects  of  the 


dopants.  The  bending-mode  energy  could  be  correlated  with  the 
size  of  the  dopant  cation,  generally  increasing  for  smaller  cations. 
NVS  measurements  of  SrZr095M005H  0,_5  (M  =  Sc,  Y,  and  Nd) 
indicated  differences  in  the  OH  bending-mode  energies  between 
the  cerates  and  zirconates.  These  differences  reflect  changes  in  the 
lattice  potential  experienced  by  the  protons,  which  ultimately  effects 
the  proton  jump  rates  and  therefore  the  performance  of  these  materi- 
als for  use  in  fuel  cells. 

An  understanding  of  the  dynamics  of  the  undoped 
perovskites  is  the  first  step  towards  a  comprehensive  picture  of 
the  protonated  doped  materials.  To  explore  this  in  more  detail,  we 
performed  ab  initio  total-energy  calculations  using  the  CAmbridge 
Serial  Total  Energy  Package  (CASTEP).  We  calculated  the 
Q  =  0  phonons  in  the  primitive  unit  cell  (i.e.,  Zr4Sr4012)  of  undoped 
strontium  zirconate  and  compared  with  the  experimental  density  of 
states  (DOS)  at  room  temperature,  as  shown  in  Fig.  1.  While  at  low 
energies,  the  calculated  modes  are  almost  at  the  same  energies  as 
observed  modes,  at  high  energies,  we  overestimate  the  energies  of 
the  modes  by  about  10  %.  However  there  is  still  a  one-to-one  cor- 
respondence between  the  main  experimental  and  calculated  spectral 
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FIGURE  1.  LEFT:  A  view  of  the  structure  of  undoped  SrZr03  in  the  [001]  direc- 
tion. RIGHT:  Vibrational  spectra  for  SrZr03  measured  by  neutron  time-of-flight 
spectroscopy.  The  space  group  and  the  lattice  parameters  (in  A)  are  also 
shown.  Numbers  in  parentheses  are  the  theoretical  values.  The  lower  solid 
line  is  the  DOS  of  the  Q  =  0  phonons  calculated  from  first  principles. 
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features,  allowing  us  to  identify  the  modes  observed  by  our  NVS 
measurements. 

In  order  to  model  the  dynamics  of  protons  trapped  in 
Sc-doped  SrZrO,,  we  replaced  one  of  the  Zr  atoms  in  the 
\flX\flX  1  supercell  by  (Sc  +  H),  which  yields  a  cell  formula 
SrsZr.ScHO,4  (See  Fig.  2).  We  performed  calculations  for  protons 
at  either  the  "undoped"  (U)  or  "doped"  (D)  sites.  Even  though  the 
MOfe  octahedra  are  quite  rigid,  the  distortions  due  to  the  presence  of 
the  proton  at  these  sites  are  quite  large. 

These  distortions  are  also  reflected  by  the  vibrational  spec- 
trum of  the  proton.  The  two  tangential  OH-  bending  modes  depend 
strongly  on  the  proton  siting.  At  the  Sc  site,  the  lowest  tangential 
mode  is  found  at  122.9  meV  However  at  the  Zr  site,  the  mode  is 
much  softer  at  88.5  meV  Interestingly,  we  observe  new  features 
in  our  NVS  spectra  at  these  energies  upon  proton  addition.  To 
further  investigate  this  effect,  we  performed  "embedded  cluster" 
calculations  in  which  the  vibrational  spectrum  of  the  H-M06  cluster 
is  calculated  while  all  other  atoms  are  kept  at  their  equilibrium 
positions.  The  similarity  between  experiment  and  the  calculated 
spectrum  from  H-Zr06  +  H-Sc06  clusters  suggests  that  the  mode 
observed  near  120  meV  is  due  to  protons  trapped  at  the  Sc  sites 


:  H  at  the  undoped  site 


while  some  portion  of  the  peak  near  80  meV  is  due  to  protons  at 
the  undoped  sites.  Thus  NVS  can  be  used  to  determine  the  hydrogen 
occupancy  of  the  various  sites. 

For  a  better  understanding  of  the  protonic  conduction  in 
these  materials,  we  calculated  the  total  energy  of  the  system  as  the 
proton  migrates  from  Zr  to  Sc  sites.  We  found  that  the  doped  site 
has  a  much  lower  energy  (-1.13  eV)  than  the  undoped  site.  The 
two  sites  are  separated  by  an  energy  barrier  of  1.5  eV.  Hence,  we 
expect  that  most  of  the  protons  are  trapped  at  the  dopant  site.  This 
is  primarily  because  the  proton  prefers  to  be  closer  to  the  Sc  cation, 
which  has  a  charge  of  +3  compared  to  +4  for  Zr.  Yet,  one  also 
has  to  consider  steric  effects;  i.e.,  if  the  dopant  cation  has  a  larger 
radius,  then  the  proton  may  prefer  to  occupy  the  Zr  site,  despite 
the  larger  Coulomb  interaction.  Thus  these  calculations  indicate 
that  the  protonic  transport  is  sensitive  to  the  competition  between 
short-range  repulsive  and  Coulomb  interactions,  and  suggest  that 
the  use  of  a  large  dopant  cation  is  one  important  step  in  developing 
a  HTPC  with  increased  protonic  conductivity.  We  are  currently 
testing  the  usefulness  of  the  calculated  proton  potentials  for  predict- 
ing the  protonic  diffusional  motions  observed  experimentally  via 
quasielastic  scattering  measurements. 
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FIGURE  2.  LEFT:  Optimized  structures  when  the  proton  is  trapped  at  the  undoped 
site  (top)  and  at  the  doped  site  (bottom).  MIDDLE:  Comparison  of  the  neutron 
vibrational  spectrum  of  SrCe^Sc.HOj  (top)  and  SrZr,  Sc^Oj  (calculated). 
Dashed  and  dotted  lines  shown  at  the  bottom  are  the  contributions  from  the 
H-M0S  clusters,  where  M  =  Zr  (undoped  site)  and  Sc  (doped  site),  respectively. 
RIGHT:  Potential  energy  of  the  crystal  as  the  proton  migrates  from  the  undoped 
site  (U)  to  the  doped  site  (D). 
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MAGNETIC  TRAPPING  OF 
ULTRACOLD  NEUTRONS 


The  demonstration  of  three-dimensional  magnetic  trapping  of 
neutrons  was  performed  this  past  year  at  NIST  [1].  The  tech- 
niques developed  should  lead  to  improved  precision  in  the  measure- 
ment of  the  neutron  beta-decay  lifetime,  thereby  expanding  our 
knowledge  of  the  weak  nuclear  force  and  our  understanding  of  the 
creation  of  matter  during  the  Big  Bang. 

Magnetic  traps  are  formed  by  creating  a  magnetic  field  mini- 
mum in  free  space.  The  confining  potential  depth  is  determined  by 
the  magnetic  moment  of  the  neutron  and  the  difference  between  the 
magnitude  of  the  field  at  the  edge  of  the  trap  and  at  the  minimum.  A 
neutron  in  a  low-field-seeking  state  (one  with  its  magnetic  moment 
anti-aligned  with  the  local  magnetic  field  vector)  feels  a  force  push- 
ing it  towards  the  trap  minimum  and  will  remain  confined  within 
the  trapping  region. 

In  order  to  load  a  neutron  into  a  static  conservative  trap,  its 
energy  must  be  lowered  while  it  is  in  the  potential  well.  We  rely 
on  a  loading  technique  that  employs  the  "superthermal  process"  [2]. 
Superfluid  helium  fills  the  trapping  region  and  serves  as  a  neutron 


scattering  medium.  A  neutron  with  kinetic  energy  near  1 1  K  (where 
the  free  neutron  and  Landau-Fey nman  dispersion  curves  cross)  that 
passes  through  the  helium  can  lose  nearly  all  of  its  energy  in  a 
single  scattering  event.  Neutrons  that  scatter  to  energies  less  than 
the  trap  depth  ( 1  mK)  and  in  the  appropriate  spin  state  are  trapped. 
Neutrons  in  this  energy  range  are  called  ultracold  neutrons  (UCN). 

Isotopically  pure  superfluid  helium  is  contained  in  a  tube 
located  inside  the  superconducting  magnet  and  centered  axially 
within  its  trapping  field  (see  Fig.  1).  The  superconducting  magnet, 
trapping  region,  and  other  key  parts  of  the  apparatus  reside  within 
a  cryogenic  dewar.  The  incident  neutron  beam  is  collimated,  passes 
through  the  trapping  region,  and  is  absorbed  by  the  beam  stop.  As 
the  beam  traverses  the  trapping  region,  about  1  %  of  the  neutrons 
scatter  in  the  helium.  Some  of  these  neutrons  are  trapped  and  the 
remainder  are  absorbed  by  shielding  materials  that  surround  the 
helium.  Low-field-seeking  UCN  are  trapped  and  remain  in  the  trap- 
ping region  until  they  decay. 
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FIGURE  1 .  Half-section  view  of  the  neutron  trapping  apparatus.  A  beam 
of  cold  neutrons  enters  from  the  left,  is  collimated,  passes  through  the 
trapping  region  and  is  absorbed  at  the  rear.  Scattered  neutrons  (yellow)  in 
the  low-field-seeking  spin  state  and  with  energy  below  the  trap  depth  are 
magnetically  confined.  Electrons  from  neutron  beta-decay  create  EUV 
scintillations  in  the  superfluid  helium  which  are  wavelength  shifted  to  the 
visible  and  transported  to  the  photomultiplier  (to  the  right). 
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The  trapped  neutrons  are  detected  when  they  decay.  When  a 
trapped  neutron  decays  (into  an  electron,  proton  and  anti-neutrino), 
the  resulting  high-energy  electron  travels  through  the  helium  leav- 
ing a  trail  of  ionized  helium  atoms.  These  ions  quickly  combine 
with  neutral  helium  atoms  to  form  metastable  diatomic  molecules. 
Most  of  these  molecules  decay  within  10  ns,  emitting  a  pulse  of 
light  in  the  extreme  ultraviolet  (EUV),  X  «=  70  nm  to  90  nm.  This 
pulse  of  scintillation  light  is  the  signal  of  a  neutron  decaying  in 
the  trap. 

In  an  experimental  run,  the  cold  neutron  beam  is  allowed  to 
pass  through  the  trapping  region,  after  which  the  beam  is  blocked 
and  pulses  of  light  are  counted.  A  background  signal,  with  both 
constant  and  time-varying  components,  obscures  the  trapped  neu- 
tron signal.  These  backgrounds  are  subtracted  by  collecting  data 
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FIGURE  2.  Counting  rate  as  a  function  of  time  after  the  neutron  beam  is 
turned  off  (pooled  background  subtracted  data),  (a)  Trapping  data  set  I, 
N.  =  560  ±  160,  (b)  Trapping  data  set  II,  NM  =  240  ±  65,  (c)  Combined  trapping 
signal,  t  =  750:^  s,  (d)  Combined  3He  data,  N3He=  53  ±  63. 
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where  the  magnetic  field  is  on  during  the  initial  loading  phase,  so 
that  UCN  are  confined  by  the  trap  and  subtracting  data  where  the 
magnetic  field  is  off  initially,  so  that  no  neutrons  are  trapped.  Equal 
numbers  of  each  data  set  were  taken,  pooled  and  subtracted  to  give 
the  background  subtracted  data. 

Two  sets  of  background  subtracted  data  were  collected:  set 
I  with  a  trap  depth  of  0.76  mK  (Fig.  2a)  and  set  II  with  a  lower 
trap  depth  of  0.50  mK  (Fig.  2b)  due  to  problems  with  the  magnet. 
For  each  set,  the  pooled  background  subtracted  data  are  modeled  to 
extract  the  amplitude  and  lifetime  of  the  decaying  neutron  signal. 
The  best  fit  values  for  the  initial  counting  rates  combined  with  the 
measured  detection  efficiency  gives  N  =  560  ±  160  and  Nn  =  240 
±  65.  Calculations  using  the  known  beam  flux,  trap  geometry  and 
the  theory  of  the  superthermal  process  predict  N,  =  480  ±  100  and 
N  =  255  ±  50,  in  good  agreement  with  the  measured  values.  The 
best  fit  value  for  the  trap  lifetime,  t  =  750  loo  s,  is  consistent  with 
the  presently  accepted  value  of  the  neutron  beta-decay  lifetime  of 
886.7  ±  1.9  s  [3]. 

To  verify  that  our  signal  is  in  fact  due  to  trapped  neutrons,  we 
doped  the  isotopically  pure  4He  with  'He  at  a  concentration  of  2  x 
10"7  3He/4He.  This  amount  of  'He  absorbs  the  trapped  neutrons  in 
less  than  1  s  without  affecting  anything  else  in  the  experiment  (less 
than  1  %  of  the  cold  neutron  beam  is  absorbed  by  3He).  The  data 
(Fig.  2d)  is  modeled  with  the  lifetime  fixed  at  =  750  s,  yielding 
N,„  =  53  ±  63,  consistent  with  zero.  This  confirms  that  our  signal 

3He 

is  due  to  trapped  neutrons. 

Magnetic  trapping  of  neutrons  is  a  new  technique  that  should 
allow  a  higher  precision  measurement  of  the  neutron  lifetime,  and 
offers  the  prospect  of  precision  much  greater  than  the  current  limit 
of  one  part  in  103.  In  order  to  realize  the  potential  of  this  technique, 
we  are  in  the  process  of  improving  our  apparatus  in  many  ways, 
including  increasing  the  size  and  depth  of  the  trap.  We  expect  to 
make  a  competitive  neutron  lifetime  measurement  in  the  near  future. 
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CHARGE  DISPROPORTIONATION  AND 
MAGNETIC  ORDERING  IN  CaFeO. 


Transition  metal  oxides  adopting  the  perovskite  crystal  structure 
(or  one  of  its  relatives)  have  occupied  a  place  in  the  scientific 
limelight  for  nearly  fifty  years  now.  Beginning  with  the  discovery 
of  ferroelectricity  in  BaTiOv  following  World  War  H,  the  dielectric 
properties  of  perovskites  have  been  extensively  studied  over  this 
entire  period.  Today  perovskite  dielectrics,  due  to  their  widespread 
usage  in  telecommunication,  are  still  extensively  studied.  The  dis- 
covery of  superconductivity  first  in  doped  BaBiO  ,  and  later  in 
the  cuprates,  such  as  YBa2Cu30?  ,  triggered  an  unprecedented 
avalanche  of  scientific  activity  in  the  late  1980's  and  early 
1990's.  More  recently,  phenomena  such  as  collosal  magnetoresi- 
tance  (CMR),  charge,  orbital  and  spin  ordering,  and  phase  separa- 
tion in  the  manganate  perovskites,  (Ln(  xAx)Mn03  (Ln  =  lanthanide 
ion,  A  =  alkaline  earth  ion),  have  captured  the  imagination  of  the 
condensed  matter  scientific  community. 

Ultimately  the  electronic  and  magnetic  properties  of  a  mate- 
rial depend  upon  the  behavior  of  the  outermost  or  valence  electrons. 
In  a  first  row  transition  metal  oxide  the  valence  electrons  are  shared, 
though  not  equally,  between  the  3d  orbitals  of  the  transition  metal 
ion  and  the  2p  orbitals  of  oxygen.  In  compounds  with  the  perovskite 
structure  the  oxygen  atoms  are  arranged  in  an  octahedral  geometry 
about  the  transitional  metal  ion;  each  oxygen  is  then  shared  by 
two  transition  metal  ions  to  form  a  three  dimensional  network  of 
corner  sharing  M06  octahedra.  The  octahedral  coordination  removes 
the  energetic  degeneracy  of  the  3d  orbitals,  forming  the  familiar 
triply  degenerate  t2g  and  doubly  degenerate  eg  set  of  orbitals.  The 
eo  orbitals  point  directly  at  the  oxygen  ligands  to  form  a  strongly 
antibonding  a*  bond.  In  contrast,  the  t2g  orbitals  have  a  smaller 
overlap  with  the  oxygen  2p  orbitals,  which  leads  to  a  weakly 
antibonding  it*  band.  If  this  covalent  bonding  interaction  between 
the  transition  metal  and  oxygen  is  weak,  the  valence  electrons 
are  localized  and  a  magnetic  insulator  is  typically  observed.  If 
we  increase  the  strength  of  the  interaction  sufficiently,  either  by 
increasing  the  electronegativity  of  the  transition  metal  ion  or  the 
spatial  overlap  of  the  metal  3d  and  oxygen  2p  orbitals,  partial 
derealization  of  the  valence  electrons  can  occur.  Frequently,  this 
leads  to  metallic  conductivity.  In  a  number  of  transition  metal 
oxides  the  metal-oxygen  interaction  strength  is  such  that  delocalized 


FIGURE  1 .  The  low  temperature  structure  of  CaFe03.  The  green  spheres  depict 
the  calcium  ions.  The  other  shaded  objects  represent  iron  centered  octahedra, 
where  an  oxygen  atom  can  be  found  at  each  vertex  and  an  iron  atom  at 
the  center  of  each  octahedron.  The  purple  and  blue  objects  represent  Fe3* 
centered  and  Fe5+  centered  octahedra  respectively. 


electrons  and  metallic  conductivity  are  seen  at  higher  temperatures 
but  give  way  to  a  magnetically  ordered,  insulating  ground  state 
upon  cooling. 

LaMn03,  SrFe03  and  CaFe03  are  isolectronic;  they  each  have 
an  electronic  configuration  of  t2g3eg'.  According  to  the  Jahn-Teller 
theorem,  the  presence  of  a  single  localized  electron  in  the  doubly 
degenerate  eg  set  of  orbitals  is  not  a  stable  situation.  LaMn03 
responds  to  this  instability  by  undergoing  a  cooperative  Jahn-Teller 
distortion  of  the  Mn06  octahedra,  producing  two  long  and  four  short 
Mn-0  bonds,  Thereby,  removing  the  degeneracy  of  the  eg  orbitals. 
SrFe03  takes  a  different  approach,  by  delocalizing  its  eg  electrons 
to  form  a  a*  band.  The  contrasting  behavior  of  these  two  materials 
can  be  understood  in  terms  of  the  strength  of  the  metal-oxygen 
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interaction.  Contrasting  LaMnO,  to  SrFeO,,  the  oxidation  state  of 
the  transition  metal  increases  from  +3  to  +4,  which  leads  to  an 
increase  in  the  metal-oxygen  interaction  strength,  increasing  the 
width  of  the  a*  band  and  stabilizing  a  metallic  ground  state.  Since 
they  both  contain  Fe4*,  one  might  expect  that  CaFeO,  and  SrFeO, 
would  behave  very  much  alike.  However,  Ca2+  is  smaller  than  Sr2"1", 
which  causes  a  tilting  of  the  Fe06  octahedra  to  satisfy  the  valence 
requirements  of  calcium.  This  seemingly  subtle  distortion  has  sev- 
eral important  consequences.  The  Fe-O-Fe  bond  is  distorted  away 
from  the  linear  geometry  observed  in  SrFeO,  (the  Fe-O-Fe  angle 
is  158°  in  CaFe03).  This  reduces  the  spatial  overlap  of  the  Fe  eg 
and  0  2p  orbitals,  and  the  width  of  the  a*  band  decreases.  The 
reduction  in  bandwidth  triggers  an  electron  localization  that  occurs 
just  below  room  temperature  (290  K).  Once  electron  localization 
occurs,  a  cooperative  Jahn-Teller  distortion  to  the  LaMn03  crystal 
structure  is  expected. 

However,  CaFeO,  refuses  to  conform  with  expectations. 
Instead  CaFeO,  undergoes  a  charge  disproportionation  (CD),  IFe4* 
(t2  V)  ->  Fe3+  (t  Jee2)  +  Fe5+  (t2g3e  °).  Evidence  for  a  CD  in  CaFe03 
was  first  proposed  over  20  years  ago,  based  on  Mossbauer  studies 
[1].  However,  crystallographic  confirmation  of  this  rare  phenom- 
enon has  proven  elusive.  Through  the  combined  use  of  synchrotron 
x-ray  (X7a-NSLS)  and  neutron  powder  diffraction  (BT-1  at  NCNR), 
we  have  elucidated  the  crystal  structure  of  CaFe03  in  its  CD  state 
for  the  first  time.  The  resulting  structure  (Fig.  1)  clearly  shows 
the  presence  of  two  chemically  and  crystallographically  distinct  Fe 
sites.  The  average  Fe-0  bond  length  about  the  "Fe5+"  site  is  1.872(6) 
A,  while  the  same  distance  about  the  "Fe3+"  site  is  1.974(6)  A. 
The  ordered  arrangement  of  Fe3+/Fe5+  is  such  that  each  Fe3+  is  sur- 
rounded by  six  Fe5+  ions,  and  vice  versa  (NaCl  or  G-type  ordering), 
optimizing  the  Coulomb  stabilization  of  the  CD  state.  Alternately, 
the  CD  process  can  be  viewed  as  the  condensation  of  a  breathing 
phonon  mode. 

In  addition  to  electronic  properties,  there  has  long  been  both 
technological  and  fundamental  scientific  interest  in  the  magnetic 
properties  and  interactions  in  perovskites  containing  transition  metal 
ions.  Goodenough's  study  of  the  various  antiferromagnetic  struc- 


tures observed  in  the  (La1  xCax)Mn03  system  was  vital  in  the  devel- 
opment of  the  Goodenough-Kanamori  rules  of  superexchange  [2]. 
These  rules  correctly  predict  an  A-type  antiferromagnetic  structure 
for  LaMnO,.  Using  the  same  rules,  we  would  expect  CaFe03  to 
have  a  simple  ferromagnetic  ground  state.  Perhaps  not  surprisingly, 
our  investigation  showed  that  again  CaFeO,  defies  expectations. 
Low  temperature  neutron  powder  diffraction  data  reveals  instead 
an  incommensurate  antiferromagnetic  ground  state  (TN     120  K). 
Analysis  of  the  data  shows  the  magnetic  structure  to  be  either 
a  screw  spiral  or  a  sinusoidal  amplitude-modulated  structure.  In 
either  case,  it  would  appear  that  a  long  range  AFM  interaction 
(probably  between  next-neighbors)  is  present  in  addition  to  the 
nearest  neighbor  FM  superexchange  interaction. 
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BOND  VALENCE  ANALYSIS 
OF  RUTHENATES 


B 


I  ond  valence  analysis  allows  us  to  evaluate,  for  any  compound 
'whose  bonding  scheme  is  known  or  assumed,  the  bond  dis- 
tances that  the  atoms  would  form  in  an  ideal  structure  in  which 
all  atomic  valences  are  exactly  balanced.  In  a  significant  number 
of  cases  this  information  is  sufficient  to  make  accurate  predictions 
about  the  real  crystal  structure  of  the  compound. 

The  bond  valence  method  is  based  on  two  concepts  that  can 
be  stated  in  the  following  way. 

(i)  A  bond  valence  v  =  v  is  assigned  to  a  bond  between 
atoms  i  and  j,  of  valences  V(i)  and  V(j),  so  that 


n(i) 


n(j) 


2;.Vv-7(0  and 


(1) 


where  n(i)  and  n(j)  are  the  number  of  atoms  in  the  coordination 
spheres  of  i  and  j,  respectively.  This  principle  of  local  valence 
balance  is  a  generalization  of  Pauling's  principle  of  local  charge 
balance  in  ionic  crystals,  and  is  known  as  the  valence  sum  rule  [1]. 

( ii)  The  sum  of  the  bond  valences  around  any  loop  in  the 
structure,  taken  with  alternating  signs,  is  equal  to  zero 


v„  =0. 

loop  'J 


(2) 


Equation  (2)  expresses  the  mathematical  conditions  that  result 
in  the  most  regular  distribution  of  the  valences  among  the  bonds  in 
a  structure  and  is  known  as  the  equal  valence  rule  [1].  The  system 
of  equations  (1)  and  (2)  allows  us  to  evaluate  the  valences  of  all  the 
individual  bonds  if  we  know  how  the  atoms  are  bonded  together  in 
a  structure  [2].  The  description  of  atomic  bonding  in  terms  of  bond 
valences  is  useful  because  the  length  d  of  a  bond  between  atoms 
i  and  j  is  a  function  only  of  the  bond  valence  v...  The  relationship 
between  these  two  quantities  is  expressed  by  the  empirical  formula 


-0.37  lnv„ 


(3) 


where  the  bond  valence  parameter  R..  depends  on  the  nature  and  the 
oxidation  states  of  atoms  /'  and  j  forming  the  bond,  and  expresses 
the  length  of  a  bond  of  unit  valence.  Values  of/?  can  be  evaluated 
from  the  bond  distances  of  known  structures  and  are  tabulated  for 
most  chemical  species  [3,4].  By  means  of  equations  (1-3)  we  may 
evaluate  the  expected  bond  lengths  for  any  known  atomic  configu- 


(a) 


e  Ru  o •  0 


FIGURE  1 .  Schematic  representation  of  (a)  the  4-layer  structure  of  TRu03 
[T=  0.875  Ba  +  0.125  Sr);  and  (b)  the  9-layer  structure  of  BaRuG,.  For  clarity 
only  the  Ru  and  0  atoms  are  shown  in  the  figure.  The  symbols  c  and  /  indicate 
the  layers  on  which  the  Ru06  octahedra  share  corners  and  faces,  respectively. 

ration.  These  distances  satisfy  exactly  the  valence  requirements  of 
the  atoms  and,  in  general,  differ  significantly  from  those  determined 
experimentally.  The  discrepancies  are  in  some  cases  due  to  the  elec- 
tronic behavior  of  particular  cations,  which  may  cause  distortions 
not  accounted  for  by  the  bond  valence  model  (for  example,  lone 
pair  distortions  around  cations  such  as  Bi3+  and  Pb2+,  or  Jahn-Teller 
distortions  around  Mn3+  and  Cu2+).  In  the  majority  of  cases,  how- 
ever, the  bond  lengths  calculated  with  equations  (1-3)  are  incom- 
mensurate under  the  constraints  imposed  by  the  crystal  geometry, 
and  have  to  be  stretched  or  compressed  in  order  to  fit  them  into  a 
particular  configuration.  Since  these  changes  introduce  strains  into 
the  structure,  the  process  of  adapting  the  theoretical  model  to  the 
requirements  of  a  space  group  symmetry  must  be  carried  out  in 
such  a  way  that  the  violations  of  the  bond  valence  sum  rule  and 
of  the  equal  valence  rule  are  kept  as  small  as  possible.  We  have 
recently  applied  the  concepts  discussed  above  to  the  determination 
of  the  crystal  structures  of  IRu03  (T  =  0.875Ba  +  0.125Sr)  [5]  and 
BaRuO,  [6],  using  initially  only  the  information  obtained  from  the 
indexing  of  the  neutron  diffraction  patterns  of  these  materials  (i.e., 
crystal  system  symmetry  and  lattice  parameters),  and  ignoring  any 
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other  structural  details  obtained  in  the  experimental  work.  Since  T2+ 
and  Ba2+  have  ionic  radii  similar  to  that  of  O2 ,  we  may  expect  that 
the  structures  of  IRuO,  and  BaRuO,  are  built  with  some  sphere 
packing  mechanism.  This  assumption  is  corroborated  by  the  fact 
that  in  both  cases  the  a-parameter  calculated  from  the  average  r 
of  the  ionic  radii  of  A  and  O  (A  =  T,  Ba)  for  the  composition 
AO,  is  in  good  agreement  with  the  experimental  values,  and  by  the 
fact  that  the  ionic  radius  of  Ru4+  is  quite  close  to  the  radius  of 
the  octahedral  void  formed  by  the  close  packing  of  oxygen  anions. 
The  periodicity  n  of  the  stacking  sequences  of  the  AO,  layers 
in  the  vertical  direction  of  the  c-axis,  evaluated  with  the  formula 
n  =  c/(2r\/2/3),  shows  that  7TtuO,  and  BaRuO,  have  4-  and  9-layer 
structures,  respectively,  in  which  the  Ru06  octahedra  are  related  to 
one  another  as  indicated  in  Fig.  la  and  lb. 

In  order  to  fit  this  configuration,  the  theoretical  bond  lengths 
calculated  with  equations  (1-3)  have  to  be  changed,  and  in  particular 
the  A-0  bonds  have  to  be  compressed,  on  the  average,  and  the 

Table  1.  Models  of  the  Structure  of  BaRu03  (R  3  m) 


12  3  4 

Lattice  parameters  (A)  and  atomic  positions 


a 

5.754 

5.754 

5.747(1) 

0.007 

c 

21.142 

21.626 

21.602(1) 

0.024 

X 

1/6 

0.1769 

0.1769(1) 

0.0000 

z, 

1/9 

0.1087 

0.1082(1) 

0.0005 

Z2 

2/9 

0.2185 

0.2175(1) 

0.0010 

Z3 

7/18 

0.3844 

0.3829(1) 

0.0015 

Bond  distances  (A) 

Ba(1)-0(1) 

2.877 

2.877 

2.8733(1) 

0.004 

-0(2) 

2.877 

2.938 

2.926(2) 

0.012 

Ba(2)-0(1) 

2.877 

2.988 

3.002(2) 

-0.014 

-0(2) 

2.877 

2.882 

2.880(2) 

0.002 

-0(2') 

2.877 

2.957 

2.945(3) 

0.012 

Ru(1)-0(1) 

2.034 

2.001 

2.005(2) 

-0.004 

Ru(2)-0(1) 

2.034 

1.995 

1.974(1) 

0.021 

-0(2) 

2.034 

1.995 

2.007(2) 

-0.012 

Ru-0  bonds  have  to  be  stretched.  This  process  of  relaxation  of 
the  initial  model  is  carried  out  by  modifying  the  structural  param- 
eters obtained  from  sphere  packing  geometry  in  such  way  that  the 
violations  of  equations  (1)  and  (2)  are  contained  within  reasonable 
bounds  [7].  The  results  of  the  bond  valence  analysis  of  BaRu03  are 
reported  in  Table  1 ,  where  they  are  compared  with  the  correspond- 
ing values  determined  experimentally.  Similar  results,  reported  in 
reference  [5],  were  obtained  for  7RuOr  The  agreement  between 
the  observed  and  calculated  structures  is  quite  good  for  both  com- 
pounds, and  the  differences  between  bond  distances  are  well  within 
0.02  A.  This  result  proves  that,  at  least  in  favorable  cases,  the 
bond  valence  method  may  yield  an  accurate  model  of  the  structure 
without  requiring  more  information  than  that  needed  to  index  a 
powder  pattern.  More  importantly,  however,  it  shows  that  the  need 
to  satisfy  the  valence  requirements  of  the  atoms  with  an  acceptably 
regular  distribution  of  the  bond  valences  is  the  driving  force  in 
determining  the  magnitude  and  the  direction  of  the  atomic  shifts 
allowed  by  the  symmetry,  and  that,  as  a  consequence,  non-bonded 
metal-metal  and  oxygen-oxygen  interactions  do  not  play  an  impor- 
tant role  in  the  way  in  which  the  structures  of  these  ruthenates  are 
built.  In  particular,  the  shifts  that  pull  together  the  oxygen  atoms 
forming  the  shared  faces  of  the  Ru06  octahedra  (thus  providing  a 
shielding  effect  to  Ru-Ru  interactions)  are  specifically  designed  to 
improve  the  local  valence  balance  of  the  Ru  and  O  atoms  involved 
in  the  Ru-0  bonds.  Application  of  the  method  to  structural  types 
other  than  perovskites  is  now  under  consideration. 
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Note.  1,  commensurate  structure  derived  from  sphere  packing  geometry; 
2,  model  obtained  with  the  relaxation  process  discussed  in  the  text;  3, 
experimental  results;  4,  difference  between  calculated  and  observed  values. 
The  theoretical  bond  distances  calculated  with  equations  (1-3)  are: 
1-0  =  2.932  A,  Ba-0  =  2.948  A  and  Ru-0  =  1.984  A. 
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RESIDUAL  STRESSES  IN  COLD-COILED 
HELICAL  AUTOMOMOTIVE  SPRINGS 


Most  front-wheel-drive  automotive  suspension  systems  use 
helical  springs.  The  process  chosen  to  produce  these  is, 
like  any  other  engineering  dilemma,  determined  by  quality,  perfor- 
mance, price,  environmental  issues,  etc.  Ford  Motor  Company  has 
developed  a  potentially  cost  saving  cold-coiling  process  in  which  . 
less  time  is  spent  treating  spring  metal  at  elevated  temperatures. 
The  pronounced  residual  stress  pattern  within  the  as-cold-coiled 
spring  is  undesirable  for  its  unpredictable  effect  on  fatigue  and 
corrosion  behavior,  and  Ford  evaluates  these  stresses  by  X-ray 
measurements  of  the  surface  stress  field  along  with  modeling  of  the 
internal  stresses.  The  success  of  these  at-home  procedures  requires 
an  independent  verification  of  the  actual  residual  stress  field  over 
the  cross-section  of  the  original  wire  stock.  The  only  available  well- 
established  method  for  this  is  neutron  diffraction,  and  that  is  where 
NCNR  expertise  comes  into  play. 

Generally,  there  are  two  ways  to  coil  a  spring:  hot  coiling 
and  cold-coiling.  Hot  coiling  implies  that  the  spring  is  wound  from 
stock  at  or  above  the  recrystallization  temperature.  The  strength 
and  fatigue  resistance  are  controlled  afterwards  by  an  appropriate 
heat  treatment.  Cold-coiling  means  that  the  helical  winding  takes 
place  at  a  low  temperature  after  the  spring  has  been  hardened  and 
tempered.  Cold-coiling  allows  the  high  temperature  heat  treatments 
to  take  place  on  the  bar  stock,  which  is  easier  to  handle  than  the 
coiled  end-product.  The  resulting  residual  stresses  can  be  essentially 
eliminated  by  a  relatively  low  temperature  tempering  treatment  fol- 
lowing the  cold  coiling. 

The  idea  is  to  measure  the  residual  stress  field  in  a  number  of 
specimens  that  represent  various  stages  of  the  production  process. 
Using  neutron  diffraction  one  can  determine  the  effect  of  the  prior 
processing  on  the  residual  stress  state  of  that  particular  stage  in 
the  process.  Of  equal  importance,  these  measurements  can  serve  to 
verify  well  established  elasto-plastic  models  that  are  being  used  to 
predict  the  formation  of  residual  stress.  Finally  one  can  look  for  a 
way  to  correlate  the  residual  stress  at  the  surface  to  the  stress  field 
as  a  whole. 

We  have  looked  at  three  cold-coiled  springs.  The  first  spring 
is  an  as-cold  coiled  spring.  The  second  one  is  cold-coiled  followed 
by  a  relatively  low  temper.  The  third  one  is  identical  to  the  second 
one,  but  in  addition  to  being  tempered  the  spring  has  been  com- 
pressed to  the  point  where  the  length  of  the  spring  is  equal  to  the 


FIGURE  1.  Experimental  setup.  The  three  arrows  indicate  three  directions  of 
measured  lattice  spacing  with  respect  to  the  specimen:  A  =  axial,  R  =  radial, 
T  =  tangential.  In  the  current  configuration,  the  spacing  of  planes  whose 
normal  is  parallel  to  T  is  being  measured.  The  sampling  volume  is  defined  by 
the  primary  aperture  P  and  secondary  aperture  S,  respectively. 


number  of  windings  times  the  wire  thickness.  After  this  the  spring 
was  allowed  to  relax.  A  small  part  of  this  torsion  strain  is  in  the 
plastic  region,  so  this  spring  is  slightly  shorter  than  all  the  others.  In 
the  automotive  industry  this  process  is  known  as  "bulldozing". 

The  measurements  were  carried  out  on  the  Double  Axis 
system  for  Residual  stress,  Texture  and  Single  crystal  analysis 
(DARTS)  at  beam  tube  8  (BT-8)  in  the  NIST  Center  for  Neutron 
Research.  This  instrument  is  specifically  designed  for  residual  stress 
measurements,  and  to  that  effect  is  equipped  with  very  accurately 
positioned  apertures  as  is  shown  in  Fig.  1 . 

For  these  experiments  the  apertures  were  chosen  to  allow 
a  sampling  volume  of  2  x  2  x  2  mm3.  The  residual  stress  in 
the  three  springs  across  the  cross-section  of  the  originally  14  mm 
thick  wire  stock  was  determined  by  detecting  the  diffracted  mono- 
chromatic neutron  intensity  with  a  position  sensitive  detector.  The 
neutron  wavelength  was  chosen  such  that  the  [211]  Fe  planes  would 
scatter  diffracted  intensity  over  approximately  90°.  The  specimen 
was  rotated  and  translated  in  this  geometry,  such  that  the  sampling 
volume  was  scanned  across  this  cross-section  allowing  the  elastic 
(residual)  strain  to  be  determined  from  the  small  shifts  in  scattering 
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FIGURE  2.  Contour  map  of  the  residual  stress  in  the  direction  of  the  length  of 
the  coiled  bar  stock,  plotted  on  the  bar  stock  cross-section.  With  reference 
to  Fig.  1 ,  this  is  the  tangential  direction.  This  map  represents  the  as-cold- 
coiled  spring.  The  left  and  right  side  of  the  map  represent  the  convex  and 
concave  sides  respectively. 

angle.  This  was  done  in  three  mutually  perpendicular  directions. 
These  strain  measurements  allowed  us  to  calculate  the  residual 
stress  in  three  perpendicular  directions  from  the  equations  published 
by  Allen  et  al.  [1]  in  each  of  the  three  specimens.  For  every  speci- 
men, the  stress  free  lattice  parameter  was  determined  on  the  basis 
that  the  net  force  on  the  cross-section  under  investigation  had  to 
be  zero. 

From  these  experiments  a  set  of  interesting  observations  can 
be  made.  First  is  the  notion  that  the  residual  stress  pattern  across  the 
wire  stock  in  the  as-coiled  spring  is  very  pronounced  and  exactly 
matches  what  one  would  expect  when  a  cylindrical  bar  is  plastically 
bent  into  a  hoop,  a  process  much  resembling  helical  coiling.  With 
reference  to  Fig.  2  we  note  essentially  uniaxial  residual  stress  in  the 
length  direction  of  the  original  wire  stock.  Through  the  diameter  of 
the  stock,  the  stress  goes  from  highly  compressive  at  the  convex 
side  to  highly  tensile  at  the  concave  side.  On  its  way  through 
the  cross-section  the  stress  level  changes  sign  three  times,  while 
the  maximum  compressive  and  tensile  stresses  are  -600  MPa  and 
+800  MPa  respectively.  For  the  tempered  and  bulldozed  specimens 
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FIGURE  3.  Same  as  Fig.  2,  these  data  representing  the  as-tempered  spring. 
The  data  for  the  as-bulldozed  spring,  though  not  given  here,  look  essentially 
the  same. 


the  pattern  is  roughly  the  same,  albeit  at  a  much  reduced  level.  The 
stress  range  being  from  -170  to  +160  MPa  as  depicted  in  Fig.  3. 

The  uncertainties  in  these  stress  levels  are  around  ±  30  MPa. 
This  means  that  the  bulldozing  process  does  not  introduce  addi- 
tional residual  stresses,  a  fact  that  can  be  well  understood  consider- 
ing that  the  plastic  deformation  under  pure  torsion  is  essentially 
uniform  and  thus  cannot  contribute  to  the  residual  stress  state. 

These  results  will  allow  Ford  to  correlate  their  model  pre- 
dictions and  X-ray  residual  stress  measurements  to  the  complete 
residual  stress  field.  This  constitutes  a  powerful  tool  in  optimizing 
the  parameters  of  the  spring  manufacturing  process. 

[1]  A.  J.  Allen,  M.  T.  Hutchings,  C.  G.  Windsor,  and  C.  Andreani:  "Neutron 
Diffraction  Methods  for  the  Study  of  Residual  Stress  Fields,"  Advances  in 
Physics,  34,  445-473  (1985). 
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THE  NCNR  NEUTRON 

SPIN  ECHO  SPECTROMETER 


The  first  experiments  have  been  performed  on  the  NCNR's 
Neutron  Spin  Echo  (NSE)  Spectrometer,  which  is  the  only 
spectrometer  of  its  kind  in  the  United  States.  This  cold  neutron 
spectrometer  allows  studies  of  dynamic  processes  in  macromolecu- 
lar  systems  that  are  relevant  to,  among  others,  polymer  [1]  and 
biomedical  [2]  sciences.  It  covers  a  time-scale  of  0.01  ns  to  100  ns 
and  a  length-scale  of  2  A  to  200  A  with  the  equivalent  of  extremely 
high  energy  resolution  and  moderate  Q-resolution.  Unlike  the  other 
high-resolution  inelastic  instruments  at  the  NCNR,  which  measure 
in  the  energy  domain,  the  NSE  spectrometer  measures,  in  the  time 
domain,  the  real  part  of  the  intermediate  scattering  function,  I(Q,t). 
This  is  done  by  using  the  neutron's  spin  precession  in  a  magnetic 
field  as  a  clock  to  determine  the  energy  transfer  in  the  scattering 
process. 

The  NSE  spectrometer,  developed  in  partnership  with  Exxon 
Research  and  Engineering  and  the  Forschungszentrum  Julich  in 
Germany  [3],  is  located  at  the  end  position  of  the  NG-5  guide. 
The  guide  is  tapered  horizontally  and  then  deviated  so  that  the 
spectrometer  is  out  of  the  direct  line  of  sight  of  the  reactor  core, 
thereby  reducing  the  background  and  the  radiation  load  on  the 
sample  region  [4].  The  taper  is  followed  by  a  Neutron  Velocity 
Selector  (NVS),  which  transmits  a  A\/\  =  10  %  FWHM  band 
of  neutrons  to  the  spectrometer.  The  last  element  of  the  guide 
is  a  transmission  polarizer  [5],  which  produces  a  polarized  beam 
of  neutrons  with  spin  anti-parallel  to  the  beam  direction  for  wave- 


lengths greater  than  5  A.  The  neutrons  precess,  from  the  first  tt/2 
flipper  (1)  to  the  sample,  through  a  phase  angle  that  is  determined 
by  the  time  that  the  neutron  spends  in  the  solenoidal  magnetic  field 
(4)  and  the  field  integral  along  the  neutron  path.  Near  the  sample,  a 
it  flipper  (2)  rotates  the  spin  by  180°  around  the  vertical  axis.  If  the 
solenoid  on  the  second  arm  of  the  spectrometer  provides  the  same 
field  integral  for  the  scattered  neutron  and  if  the  scattering  is  elastic, 
the  neutron  will  precess  through  the  same  phase  angle  as  along  the 
first  arm  and  will  end  up  with  the  original  spin  orientation  at  the 
final  tt/2  flipper  (3).  This  tt/2  flipper  rotates  the  spin  back  into  the 
horizontal  plane.  If  the  neutron  is  scattered  quasielastically,  it  will 
precess  a  slightly  different  number  of  times  in  the  second  arm  of  the 
spectrometer  and  end  up  at  the  second  tt/2  flipper  rotated  by  some 
angle  that  is  proportional  to  the  wavelength  shift.  The  analyzer  (7) 
projects  the  component  of  the  neutron  spin  that  is  parallel  to  the 
field  direction  onto  the  detector. 

Figure  2  shows  the  detector  count  rate  as  a  function  of  current 
in  the  phase  coil  (5),  which  changes  the  field  integral  on  one  arm 
of  the  spectrometer.  The  field  integrals  for  the  first  and  second  arms 
of  the  spectrometer  are  equal  at  the  echo  point  where  the  amplitude 
is  at  a  maximum,  since  for  an  elastic  scatterer  the  neutrons  all 
arrive  at  the  analyzer  in  phase.  As  the  phase  current  is  changed, 
the  neutron  spins  rotate  away  from  the  polarization  axis  of  the 
analyzer  (7)  and  the  count  rate  is  changed.  The  period  of  this  oscil- 
lation is  proportional  to  the  mean  wavelength.  The  envelope  of  the 
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FIGURE  1.  (Top)  Schematic  plan  view  of  the  NSE  spectrometer.  The  numbers 
refer  to  spectrometer  elements  as  described  in  the  text.  (Bottom)  Changes 
of  neutron  spin  orientation  passing  through  the  spectrometer  elements  for 
elastically  (black)  and  inelastically  (blue)  scattered  neutrons. 
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FIGURE  2.  Phase  scan  of  an  echo  for  an  elastic  scatterer  (Grafoil).  A  is  the 
amplitude  at  the  echo  point;  N0H  and  N0FF  are  counting  rates  with  the  it  flipper 
on  and  off,  respectively,  and  the  rr/2  flippers  off. 
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FIGURE  3.  Normalized  and  resolution-corrected  measurements  of  the 
intermediate  scattering  function  at  the  indicated  Q-values  for  10  %  SDS  in  D20 
at  room  temperature.  The  resolution  curve  is  a  normalized  measurement  of  an 
elastic  scatterer  (Grafoil)  at  10  A. 
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oscillations  gives  the  wavelength  distribution.  The  amplitude  A  is 
proportional  to  I(Q,t),  which  can  be  normalized  by  measuring  the 
difference  in  intensity  with  the  tt  flipper  on  and  off,  NQN-  N0FF  For 
the  normalizing  measurement,  which  is  proportional  to  I(Q,t  =  0), 
the  tt/2  flippers  are  off  since  there  is  no  precession  of  the  neutrons 
in  this  configuration.  Instrumental  resolution  effects  are  removed 
by  dividing  the  normalized  sample  signal  by  the  normalized  values 
from  an  elastic  scatterer. 

We  have  also  verified  the  operation  of  the  correction  elements 
(6)  which  allow  non-axial  and  divergent  neutrons  to  satisfy  the  echo 
condition.  These  are  essential  for  operation  at  high  fields.  The  value 
of  the  resolution  amplitude  at  25  ns  in  Fig.  3  would  be  at  least 
a  factor  of  ten  smaller  without  the  contribution  of  the  correction 
elements. 

As  an  example  of  the  science  that  is  available  through  NSE, 
we  have  reproduced  the  measurements  of  Hayter  and  Penfold  [6]  on 
a  micellar  solution  of  10  %  sodium  dodecylsulfate  (SDS)  in  D,0. 
Fig.  3  shows  the  normalized  intermediate  scattering  function  versus 
time  at  several  Q  values  at,  well  below,  and  well  above  the  interac- 
tion peak  in  the  structure  factor.  The  effective  diffusion  constant 
Deff(Q),  obtained  from  I(Q,t) ;  exp(-Dcfl(Q)Q2t),  varies  inversely 
with  the  structure  factor,  reaching  a  minimum  at  the  structure  factor 
peak,  an  effect  analogous  to  DeGennes  narrowing.  As  a  result,  in 
a  dilute  solution  of  interacting  micelles,  one  may  unambiguously 
separate  the  structure  factor  from  the  form  factor,  which  may  not  be 
done  with  only  a  SANS  measurement. 
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SERVING  THE  SCIENTIFIC  AND 
TECHNOLOGICAL  COMMUNITIES 


Providing  neutron  beam  methods  to  the  U.S.  research  com- 
munity is  a  central  purpose  of  the  NCNR.  Intense  neutron 
sources  are  scarce,  as  are  the  advanced  instruments  needed  to  fully 
exploit  them.  Expert  scientists  will  always  recognize  opportunities 
to  perform  research  with  the  most  powerful  tools  available,  so 
that  the  existence  of  the  NCNR  as  the  Nation's  premier  reactor- 
based  neutron  facility  stimulates  considerable  demand  from  the 
community.  However,  in  order  to  foster  the  best  possible  science, 
many  of  our  procedures  and  activities  are  designed  to  encourage 
researchers  to  learn  more  about  our  capabilities,  and  eventually  to 
obtain  instrument  time  through  proposals  and  collaborations.  The 
first  instruments  became  operational  in  the  NCNR  guide  hall  in 
1991,  and  a  formal  proposal  system  was  started  to  accommodate 
users  from  universities,  industry,  and  government  laboratories. 
Since  then,  the  quantity  and  the  diversity  of  research  carried  out 
at  the  NCNR  has  grown  steadily,  and  the  number  of  research 
participants  (Fig.  1)  is  now  several  times  larger  than  it  was  just  a 
few  years  ago. 
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In  FY99  Research  Participants  were  from: 
23  NIST  Divisions  and  Offices 
34  U.S.  Government  Laboratories 

105  U.S.  Universities 
55  U.S.  Industrial  Laboratories 


500 


86   87   88   89   90   91    92   93   94   95   96   97   98  99 


Fiscal  Year 


FIGURE  1.  Research  participants  at  the  NCNR. 


THE  USER  PROGRAM 

NIST  has  always  relied  on  advice  from  accomplished  scientists 
to  assist  in  formulating  policy.  The  Program  Advisory  Committee 
(PAC)  is  the  body  primarily  responsible  for  proposal  review  and 
recommending  user  policies  for  the  NCNR,  working  closely  with 
the  Center's  Director  and  staff.  Its  current  membership  includes 
Sanat  Kumar  (Penn  State  University,  chair),  Robert  M.  Briber 
(University  of  Maryland),  Michael  K.  Crawford  (DuPont),  Dieter 
K.  Schneider  (Brookhaven  National  Laboratory),  Thomas  P.  Russell 
(University  of  Massachusetts),  Sunil  K.  Sinha  (Argonne  National 
Laboratory),  Emile  A.  Schweikert  (Texas  A&M  University), 
Laurence  Passell  (Brookhaven  National  Laboratory),  and  Gabrielle 
G.  Long  (NIST). 

At  their  meeting  in  May  1999,  the  PAC  made  several  recom- 
mendations for  improving  the  current  system.  They  felt  that  pro- 
gram proposals  for  longer-term  projects  had  not  had  their  intended 
effect  of  reducing  the  burden  for  external  reviewers,  and  should 
be  discontinued  for  the  time  being.  Noting  a  potential  for  large 
increases  in  instrument  time  for  biology-related  proposals  in  the 
next  few  years,  the  addition  of  another  PAC  member  with  expertise 
in  that  area  was  recommended.  The  PAC  discussed  the  results  of 
a  recent  independent  survey  of  user  satisfaction,  coordinated  by 
Anne  Mayes  of  MIT.  Ideas  for  obtaining  more  feedback  from  users 
were  also  considered.  Incidentally,  Professor  Mayes,  the  head  of  the 
NCNR  User's  Group,  recently  won  the  American  Physical  Society's 
Dillon  Medal,  partly  for  work  carried  out  using  NCNR  instruments. 

Most  of  the  PAC  meeting  was  devoted  to  proposal  review. 
With  the  help  of  several  written  reviews  for  each  proposal  from 
selected  experts,  the  PAC  considered  151  proposals  requesting  1138 
instrument-days,  and  allocated  514  instrument-days  to  95  proposals 
for  SANS,  reflectometry,  SPINS,  and  time-of-flight  instruments. 
The  PAC  usually  meets  twice  a  year,  but  because  of  a  planned 
shutdown  for  cold  source  and  reactor  upgrades,  only  one  meeting 
was  held  in  1999. 
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FIGURE  2.  Nick  Rosov  explains  the  operation  of  the  Neutron  Spin  Echo 
Spectrometer  to  summer  school  participants. 


SUMMER  SCHOOL  ON  METHODS  AND 
APPLICATIONS  OF  NEUTRON  SPECTROSCOPY 

Each  summer  for  the  past  five  years,  the  NCNR  has  held  a  one- 
week  school  to  introduce  researchers  to  the  methods  of  neutron 
scattering.  Organized  by  senior  staff  members  John  Copley  and 
Peter  Gehring,  the  1999  summer  school  was  devoted  to  spectrosco- 
py with  neutrons,  with  an  emphasis  on  the  opportunities  afforded  by 
new  instruments  in  the  guide  hall.  Most  of  the  32  participants,  cho- 
sen from  almost  60  who  applied,  were  graduate  students  and  post- 
doctoral fellows.  The  curriculum  consisted  of  lectures  by  NCNR 
staff  and  experts  from  other  institutions,  tours  of  the  facility,  and 
four  three-hour,  hands-on,  sessions  at  instruments  in  the  guide  hall. 
Informal  discussions  among  students  and  the  resident  staff  were 
also  an  important  part  of  the  school's  activities. 

Beginning  with  lectures  on  the  fundamental  concepts  of 
nuclear  and  magnetic  scattering,  the  agenda  turned  to  talks  on  appli- 
cations of  neutron  spectroscopy  to  a  wide  range  of  research  topics. 
Subsequent  presentations  covered  the  reactor  and  cold  source, 
specific  scattering  techniques,  computer  modeling  and  ab  initio 
calculations,  and  the  operating  principles  of  various  spectrometers 
used  in  the  experimental  sessions.  In  the  latter,  students  and  staff 
used  time-of-flight,  backscattering,  triple-axis,  and  neutron  spin- 
echo  spectrometers  to  measure,  respectively,  quasielastic  scattering 
from  water,  rotational  tunneling  in  solid  methyl  iodide,  magnetic 
excitations  in  a  geometrically  frustrated  antiferromagnet,  and 


coherent  scattering  from  the  time  domain  in  a  system  of  spherical 
micelles.  (The  last  example  is  illustrated  in  the  article  on  the 
Neutron  Spin  Echo  Spectrometer.)  The  lecture  materials  will  be 
placed  on  the  NCNR  website,  in  order  to  reach  as  wide  an  audience 
as  possible. 

THE  CENTER  FOR  HIGH-RESOLUTION  NEUTRON 
SCATTERING  (CHRNS) 

Supported  by  the  National  Science  Foundation  (NSF),  CHRNS  is  a 
very  important  component  of  the  user  program.  It  operates  a  suite  of 
three  instruments,  including  a  30  m  SANS  machine,  the  SPINS  tri- 
ple-axis spectrometer,  and  a  double-crystal,  high-resolution  SANS. 
The  last  of  these  is  under  construction.  Approximately  40  %  of  the 
instrument  time  allocated  by  the  PAC  goes  to  experiments  carried 
out  on  CHRNS  instruments.  The  NSF  is  currently  reviewing  the 
CHRNS  renewal  proposal,  which  requests  support  for  upgrading 
existing  instruments  and  including  an  improved  9  m  SANS  instru- 
ment within  CHRNS.  Approval  would  mean  an  increase  of  more 
than  50  %  in  CHRNS'  capacity  to  serve  users. 

COLLABORATIONS 

Direct  collaborations  on  specific  experiments  remain  a  common 
way  for  users  to  pursue  their  ideas  using  NCNR  facilities,  account- 
ing for  approximately  60  %  of  the  number  of  instrument-days.  The 
thermal-neutron  triple-axis  spectrometers  are  mainly  scheduled  in 
this  way.  Most  of  the  time  reserved  for  NIST  on  these  and  all 
other  NCNR  instruments  is  also  devoted  to  experiments  that  are 
collaborations  with  non-NIST  users. 

Another  mode  of  access  to  the  NCNR  is  through  Participating 
Research  Teams  (PRTs).  In  this  case,  groups  of  researchers  from 
various  institutions  join  forces  to  build  and  operate  an  instrument. 
Typically,  50  %  to  75  %  of  the  time  on  the  instrument  is  then 
reserved  for  the  PRT,  and  the  remaining  time  is  allocated  to  general 
user  proposals.  For  example,  a  PRT  involving  ExxonMobil,  the 
University  of  Minnesota,  Texaco  R&D  and  NIST  cooperates  on 
the  NG-7  30  m  SANS  instrument.  Similar  arrangements  involving 
other  PRTs  apply  to  the  horizontal-sample  refiectometer,  the  high- 
resolution  powder  diffractometer,  and  the  neutron  spin-echo  spec- 
trometer. 
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INDEPENDENT  PROGRAMS 

There  are  a  number  of  programs  of  long  standing  located  at  the 
NCNR  which  involve  other  parts  of  NIST,  universities,  industrial 
laboratories,  or  other  government  agencies. 

The  Polymers  Division  of  the  Materials  Science  and 
Engineering  Laboratory  has  two  major  program  elements  at  the 
NCNR.  In  the  first,  the  purpose  is  to  help  the  U.S.  microelectronics 
industry  in  addressing  their  most  pressing  materials  measurement 
and  standards  issues.  In  today's  ICs  and  packages,  the  feature  size 
on  a  chip  is  ever  shrinking,  approaching  250  nm,  while  the  size 
of  a  polymer  molecule  is  typically  5  nm  to  10  nm.  As  feature  size 
shrinks,  the  structure  and  properties  of  interfaces  play  an  increas- 
ingly important  role  in  controlling  the  properties  of  the  polymer 
layers  used  in  interconnects  and  packages.  NIST  scientists  use  both 
neutron  reflectivity  and  other  neutron  scattering  methods  to  charac- 
terize polymer/metal  interfaces  with  regard  to  local  chain  mobility, 
moisture  absorption,  glass  transition  temperature,  and  crystalline 
structure.  In  the  second  program  element,  the  objective  is  to  under- 
stand underlying  principles  of  phase  behavior  and  phase  separation 
kinetics  of  polymer  blends,  both  in  the  bulk  and  on  surfaces,  in 
order  to  help  control  morphology  and  structure  during  processing. 
SANS  and  reflectivity  measurements  in  equilibrium,  in  transient 
conditions,  and  under  external  fields  provide  essential  information 
for  general  understanding  as  well  as  for  specific  application  of 
polymer  blend/alloy  systems.  Customers  include  material  producers 
and  users,  ranging  from  chemical,  rubber,  the,  and  automotive  com- 
panies, to  small  molding  and  compounding  companies.  The  focus 
of  research  on  polymeric  materials  includes  commodity,  engineering 
and  specialty  plastic  resins,  elastomers,  coatings,  adhesives,  films, 
foams,  and  fibers. 

The  ExxonMobil  Research  and  Engineering  Company 
is  a  member  of  the  Participating  Research  Team  (PRT)  that  operates, 
maintains,  and  conducts  research  at  the  NG-7  30  m  SANS  instru- 
ment and  the  recently  commissioned  NG-5  Neutron  Spin  Echo 
Spectrometer.  The  mission  is  to  use  those  instruments,  as  well  as 


other  neutron  scattering  techniques,  in  activities  that  complement 
research  at  ExxonMobil's  main  laboratories  as  well  as  at  its  affili- 
ates' laboratories  around  the  world.  The  aim  of  these  activities  is 
to  deepen  understanding  of  the  nature  of  ExxonMobil's  products 
and  processes,  so  as  to  improve  customer  service  and  to  improve 
the  return  on  shareholders'  investment.  Accordingly,  and  taking  full 
advantage  of  the  unique  properties  of  neutrons,  most  of  the  experi- 
ments use  SANS  or  other  neutron  techniques  to  study  the  structure 
and  dynamics  of  hydrocarbon  materials,  especially  in  the  fields 
of  polymers,  complex  fluids,  and  petroleum  mixtures.  ExxonMobil 
regards  its  participation  in  the  NCNR  and  collaborations  with  NIST 
and  other  PRT  members  not  only  as  an  excellent  investment  for  the 
company,  but  also  as  a  good  way  to  contribute  to  the  scientific  health 
of  the  Nation. 

The  Nuclear  Methods  Group  (Analytical  Chemistry 
Division,  Chemical  Science  and  Technology  Laboratory)  has  as  its 
principal  goals  the  development  and  application  of  nuclear  analytical 
techniques  for  the  determination  of  elemental  compositions  with 
greater  accuracy,  higher  sensitivity,  and  better  selectivity.  A  high 
level  of  competence  has  been  developed  in  both  instrumental  and 
radiochemical  neutron  activation  analysis  (INAA  and  RNAA).  In 
addition,  the  group  has  pioneered  the  use  of  cold  neutron  beams 
as  analytical  probes  with  both  prompt  gamma  activation  analysis 
(PGAA)  and  neutron  depth  profiling  (NDP).  PGAA  measures  the 
total  amount  of  a  particular  analyte  present  throughout  a  sample 
by  the  analysis  of  the  prompt  gamma-rays  emitted  during  neutron 
capture.  NDP  on  the  other  hand,  determines  concentrations  of 
several  important  elements  (isotopes)  as  a  function  of  depth  within 
the  first  few  micrometers  of  a  surface  by  energy  analysis  of 
the  prompt  charged-particles  emitted  during  neutron  bombardment. 
These  techniques  (INAA,  RNAA,  PGAA,  and  NDP)  provide  a 
powerful  combination  of  complementary  tools  to  address  a  wide 
variety  of  analytical  problems  of  great  importance  in  science  and 
technology,  and  are  used  to  help  certify  a  large  number  of  NIST 
Standard  Reference  Materials.  During  the  past  several  years,  a  large 
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part  of  the  Group's  efforts  has  been  directed  towards  the  exploitation 
of  the  analytical  applications  of  the  guided  cold-neutron  beams 
available  at  the  NIST  Center  for  Neutron  Research.  The  Group's 
involvement  has  been  to  design  and  construct  state-of-the-art  cold 
neutron  instruments  for  both  PGAA  and  NDP  and  provide  facilities 
and  measurements  for  outside  users,  while  retaining  and  utilizing 
our  existing  expertise  in  INAA  and  RNAA. 

The  Center  for  Food  Safety  and  Applied  Nutrition  of  the 
U.  S.  Food  and  Drug  Administration  (FDA)  directs  and  maintains 
a  neutron  activation  analysis  (NAA)  facility  at  the  NCNR.  This 
facility  provides  agency-wide  analytical  support  for  special  investi- 
gations and  applications  research,  complementing  other  analytical 
techniques  used  at  FDA  with  instrumental  (INAA),  neutron-capture 
prompt-gamma  (PGAA),  and  Radiochemical  Neutron  Activation 
Analysis  (RNAA)  procedures,  radioisotope  X-ray  fluorescence  spec- 
trometry (RXRFS),  and  low-level  gamma-ray  detection.  This  combi- 
nation of  analytical  techniques  enables  diverse  multi-element  and 
radiological  information  to  be  obtained  for  foods  and  related  materi- 
als. The  NAA  facility  supports  agency  quality  assurance  programs 
by  developing  in-house  reference  materials,  by  characterizing  food- 
related  reference  materials  with  NIST  and  other  agencies,  and  by 
verifying  analyses  for  FDA's  Total  Diet  Study  Program  annually. 
Other  studies  include  the  development  of  RXRFS  methods  for 
screening  foodware  for  the  presence  of  Pb,  Cd,  and  other  potentially 
toxic  elements,  use  of  INAA  to  investigate  bromate  residues  in 
bread  products,  and  use  of  PGAA  to  investigate  boron  nutrition 
and  its  relation  to  bone  strength.  The  FDA's  NAA  laboratory  person- 
nel frequently  provide  intra-agency  technical  assistance,  the  most 
recent  example  being  participation  in  the  production  of  the  docu- 
ment "Accidental  Radioactive  Contamination  of  Human  Food  and 
Animal  Feeds:  Recommendations  for  State  and  Local  Agencies"  by 
the  Center  for  Devices  and  Radiological  Health. 

The  Neutron  Interactions  and  Dosimetry  Group 
(Physics  Laboratory)  provides  measurement  services,  standards,  and 
fundamental  research  in  support  of  NIST's  mission  as  it  relates 


to  neutron  technology  and  neutron  physics.  The  national  and  indus- 
trial interests  served  include  scientific  instrument  calibration,  elec- 
tric power  production,  radiation  protection,  defense  nuclear  energy 
systems,  radiation  therapy,  neutron  radiography,  and  magnetic  reso- 
nance imaging.  The  Group's  activities  may  be  represented  as  three 
major  activities.  The  first  is  Fundamental  Neutron  Physics  -  includ- 
ing operation  of  a  neutron  interferometry  and  optics  facility,  devel- 
opment of  neutron  spin  filters  based  on  laser  polarization  of  -He, 
measurement  of  the  beta  decay  lifetime  of  the  neutron,  and  investi- 
gations of  other  coupling  constants  and  symmetries  of  the  weak 
interaction.  This  project  involves  a  large  number  of  collaborators 
from  universities  and  national  laboratories.  The  second  is  Standard 
Neutron  Fields  and  Applications  -  utilizing  both  thermal  and  fast 
neutron  fields  for  materials  dosimetry  in  nuclear  reactor  applications 
and  for  personnel  dosimetry  in  radiation  protection.  These  neutron 
fields  include  thermal  neutron  beams,  "white"  and  monochromatic 
cold  neutron  beams,  a  thermal-neutron-induced  235U  fission  neutron 
field,  and  252Cf  fission  neutron  fields,  both  moderated  and  unmoder- 
ated.  The  third  is  Neutron  Cross  Section  Standards  -  including 
experimental  advancement  of  the  accuracy  of  neutron  cross  section 
standards,  as  well  as  evaluation,  compilation  and  dissemination  of 
these  standards. 

Several  universities  have  also  established  long  term  programs 
at  the  NCNR.  The  University  of  Maryland  is  heavily  involved 
in  the  use  of  the  NCNR,  and  maintains  several  researchers  at 
the  facility.  Johns  Hopkins  University  participates  in  research 
programs  in  solid  state  physics  and  in  instrument  development  at  the 
NCNR.  The  University  of  Pennsylvania  is  working  to  help  develop 
biological  applications  of  neutron  scattering.  It  is  also  participating 
in  the  construction  of  a  new  filter  analyzer  neutron  spectrometer, 
along  with  the  University  of  California  at  Santa  Barbara,  DuPont, 
Hughes,  and  Allied  Signal.  The  University  of  Minnesota  partici- 
pates in  two  PRTs,  the  NG-7  30  m  SANS  and  the  NG-7  refiectom- 
eter.  The  University  of  Massachusetts  also  participates  in  the  latter 
PRT 
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REACTOR  OPERATIONS 
AND  ENGINEERING 


!he  reactor  operated  for  250  full  power  (20  MW)  days  during  the 
past  year,  or  68  %  of  real  time.  This  meets  our  goal,  and  was 


achieved  in  spite  of  a  shutdown  for  maintenance  early  in  the  fiscal 
year.  The  reactor  operated  with  better  than  90  %  predictability— i.e. , 
days  operated  on  the  day  originally  scheduled.  This  is  substantially 
better  operation  than  is  specified  in  national  goal  of  85  %  reliability 
for  user  facilities.  Just  after  the  end  of  the  fiscal  year,  another  fuel 
shipment  was  made,  thus  clearing  additional  space  in  the  spent  fuel 
storage  pool.  This  completes  our  current  round  of  scheduled  spent 
fuel  shipments.  There  is  now  ample  storage  space  for  at  least  5  years 
of  continuous  20  MW  operation. 


TRTR  Chairman-Elect,  Pedro  Perez,  presents  NIST  Director,  Raymond  G.  Kammer, 
with  a  plaque  stating: 

"The  National  Organization  of  Test,  Research  and  Training  Reactors  recognizes 
the  National  Institute  of  Standards  and  Technology  on  the  occasion  of  its 
centennial  as  the  foremost  measurements  and  standards  laboratory  in  the 
world,  and  honors  NIST  for  100  years  of  distinguished  service  to  the  nation  and 
to  all  the  segments  of  science  and  technology,  and  particularly  for  its  pioneer- 
ing efforts  in  nuclear  research  and  radiation  standards  and  in  appreciation  of 
its  leadership  in  the  formation  of  TRTR." 


As  part  of  preparing  for  another  20  years  of  operation  past 
2004,  we  have  been  assessing  the  present  cooling  tower.  During  the 
past  year,  a  contract  for  construction  of  a  hybrid  wet/dry  tower  was 
signed.  This  new  tower  will  not  only  meet  our  cooling  needs;  it 
will  also  reduce  the  visible  plume  given  off  by  the  tower  during  the 
winter  months,  a  substantial  additional  benefit.  The  new  tower  will 
be  erected  adjacent  to  the  current  one,  and  will  not  require  an  outage 
until  the  switch  is  made.  This  outage  will  be  scheduled  when  the  new 
tower  is  fully  ready  for  use,  so  as  to  minimize  down  time.  Extensive 
planning  for  the  piping  changes  that  will  be  needed  is  now  underway. 
It  is  expected  that  the  new  second-generation  cold  source  will  be 
installed  at  the  same  time,  so  as  to  maximize  operating  time. 

During  the  coming  year,  we  will  take  advantage  of  the  normal 
Christmas  shutdown  to  replace  the  shim  control  arms.  This  operation 
is  required  every  4  to  5  years,  since  the  cadmium  in  the  control 
arms  is  burned  up  in  service.  By  utilizing  the  time  period  of  the 
holiday  shutdown,  the  entire  operation  will  be  done  with  only  3 
weeks  of  additional  shutdown  time.  Since  this  period  will  also  be 
used  for  required  training  and  requalification  of  operations  personnel, 
there  will  be  little  lost  operating  time  compared  to  our  best  possible 
schedule. 

A  major  highlight  of  the  year  was  the  annual  conference  of 
the  National  Organization  of  Test,  Research  and  Training  Reactors 
(TRTR)  held  at  NIST,  chaired  and  organized  by  the  group.  By  all 
accounts,  it  was  highly  successful  with  the  largest  number  of  partici- 
pants ever.  Many  distinguished  speakers  addressed  the  conference 
including  the  Chairman  of  the  Nuclear  Regulatory  Commission  and 
the  Under  Secretary  of  Energy.  Also,  TRTR  became  the  first  profes- 
sional organization  to  honor  NIST  on  its  forthcoming  centennial. 

Finally,  preparations  for  a  license  renewal  application  are 
proceeding  as  planned.  We  are  fortunate  to  have  attracted  Dr. 
Seymour  H.  (Sy)  Weiss  to  our  staff,  as  Deputy  to  NCNR  Director.  He 
will  be  leading  the  relicensing  effort  and  associated  reactor  upgrades 
for  the  coming  years,  and  has  already  started  several  tasks  related 
to  this  effort. 
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INSTRUMENTATION 
DEVELOPMENTS 


IMPROVEMENTS  TO  THE 
PERFORMANCE  OF  HFBS 

The  recent  commissioning  of  the  high  flux  backscattering  spectrom- 
eter (HFBS)  has  opened  up  opportunities  to  a  variety  of  users  from 
around  the  country  to  study  the  dynamics  of  condensed  matter 
systems  with  time  scales  on  the  order  of  nanoseconds  (10  9s)  and 
length  scales  up  to  10  A.  A  careful  design  which  incorporates  state- 
of-the-art  neutron  optics  has  resulted  in  a  flux  on  sample  as  high 
as  that  of  any  other  backscattering  spectrometer  in  the  world  with 
comparable  energy  resolution.  The  first  call  for  proposals  for  the 
HFBS  was  made  in  FY  1999  and  the  response  was  outstanding. 
Several  successful  user  experiments  have  already  been  performed 
(see  Fig.  1)  with  more  scheduled  in  the  near  future.  Although  the 
initial  performance  of  the  instrument  was  impressive,  improvements 
continue  to  be  made. 

One  of  the  recent  changes  to  the  instrument  has  been  a  modi- 
fication to  the  data  acquisition  system.  Initially  the  data  were  binned 
to  the  velocity  of  the  reciprocating  silicon  monochromator.  However 
the  method  in  which  the  velocity  was  determined  overemphasized 
local  details  of  the  monochromator  motion  which  were  not  observed 


FIGURE  1.  Measurement  of  quasi-elastic  scattering  from  monolayer  coverage 
of  alkanes  on  grafoil  performed  on  the  HFBS  by  H.  Taub,  D.  Fuhrmann,  L. 
Criswell,  and  K.  Herwig.  The  low  temperature  data  (blue)  is  resolution  limited 
with  a  full-width  at  half  maximum  of  0.96  \ie\l.  The  higher  temperature  data 
(red)  clearly  displays  broadening  indicative  of  the  diffusive  motion  of  the 
alkanes. 


in  the  neutron  data.  A  detailed  analysis  of  the  performance  of  the 
monochromator  system  indicated  that  the  quality  of  the  data  could 
be  improved  by  making  a  rather  simple  change  in  the  way  that 
data  is  collected  and  stored,  namely  by  binning  data  versus  time 
during  a  period  of  the  Doppler  monochromator  motion.  The  data 
are  later  re-binned  to  velocity  using  the  average  monochromator 
motion  profile.  This  straightforward  change  to  the  way  that  the  data 
are  collected  has  improved  the  quality  of  the  data,  suppressing  the 
spurious  effects  observed  in  the  velocity-binning  mode. 

Improvements  have  also  been  made  to  the  Doppler  mono- 
chromator system  which  have  improved  its  reliability  and  perfor- 
mance. One  of  the  problems  encountered  is  associated  with  the 
fact  that  the  monochromator  is  being  driven  to  higher  energy  trans- 
fers than  any  other  Doppler  monochromator  system  in  the  world. 
The  monochromator,  which  is  vibrated  over  a  distance  of  9  cm  at 
high  frequencies,  can  experience  accelerations  in  excess  of  100  g's. 
Prolonged  vibration  of  the  monochromator  resulted  in  the  silicon 
wafers  coming  off  of  the  monochromator  surface.  Recent  changes 
in  the  gluing  method  have  dramatically  increased  the  lifetime  of  the 
silicon  wafers  on  the  monochromator. 

Efforts  to  optimize  shielding  in  various  parts  of  the  instru- 
ment have  resulted  in  improvements  in  the  signal  to  noise  ratio  in 
the  detectors.  During  normal  operation  the  instrument  vessel  (which 
is  made  up  of  the  sample,  analyzer  crystals,  and  the  detectors)  is 
evacuated  in  order  to  reduce  the  background  due  to  air  scattering. 
Additional  shielding  has  been  installed  on  the  detector  assembly 
and  on  parts  of  the  analyzer  system  in  the  vessel.  Improvements 
have  also  been  made  to  the  shielding  on  the  phase  space  transform 
chopper  to  handle  the  highly  divergent  beam  from  the  converging 
neutron  guide.  This  has  resulted  in  a  substantial  increase  of  the 
signal  to  noise  ratio  to  almost  600  to  1 . 

THE  DISK  CHOPPER  SPECTROMETER 

The  Disk  Chopper  Spectrometer  (DCS)  is  an  extremely  versatile 
time-of-flight  instrument,  designed  to  achieve  a  broad  range  of 
energy  resolution  full  widths,  from  12  jxeV  to  1  meV,  by  changing 
chopper  speed,  wavelength,  and/or  beam  width.  Once  commis- 
sioned, the  DCS  will  be  in  high  demand  for  experiments  on  a 
variety  of  systems  such  as  proteins,  molecular  crystals,  disordered 
materials,  and  metal-hydrogen  systems.  The  instrument  is  to  be 
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FIGURE  2.  Photograph  of  the  DCS  detector  bank  and  flight  chamber. 

included,  on  a  limited  basis,  in  the  first  Call  for  Proposals  in 
FY2000. 

Over  the  past  twelve  months  the  DCS  has  seen  major  prog- 
ress in  a  number  of  important  areas.  The  sample  chamber  and  flight 
chamber  were  equipped  with  beam-handling  components,  and  the 
insides  of  both  chambers  were  lined  with  cadmium.  Flow  tests  were 
performed  in  order  to  complete  the  design  of  the  gas  handling 
system  for  the  sample  and  flight  chambers;  the  system  had  been 
fabricated  and  had  largely  been  assembled  by  the  end  of  the  fiscal 
year.  An  overhead  platform,  to  be  used  for  sample  environment 
preparation  and  installation  into  the  sample  chamber,  has  also  been 
designed. 

Problems  with  some  of  the  amplifiers  for  the  DCS  detectors 
prompted  us  to  send  the  complete  inventory  to  the  manufacturer 
for  mutually  agreed  modifications.  On  their  return  discriminator 
thresholds  were  individually  set,  and  the  full  complement  of  913 
detectors  and  amplifiers  was  installed  on  the  detector  racks,  which 
were  then  installed  and  aligned  at  the  spectrometer.  The  detectors 
and  voltage  distribution  assemblies  were  cabled  and  fully  tested 
prior  to  installation  of  the  outer  shields.  Much  work  was  devoted 
to  troubleshooting  and  greatly  improving  the  reliability  and  perfor- 
mance of  the  VME  data  acquisition  electronics,  and  to  the  develop- 
ment of  data  acquisition  software. 


PERFECT  CRYSTAL  SANS  DIFFRACTOMETER 

A  perfect  crystal  diffractometer  (PCD)  used  for  very  high  resolution 
small  angle  neutron  scattering  (SANS)  is  being  developed  jointly  by 
the  NCNR  and  NSF  as  part  of  the  CHRNS  facility.  The  higher  reso- 
lution obtained  by  using  perfect  silicon  crystals  increases  the  maxi- 
mum size  of  features  that  can  be  measured  from  0. 1  \im  obtained 
using  the  current  NCNR's  two  30  m,  pinhole  geometry  SANS 
instruments,  to  10  \xm  with  the  new  instrument.  The  monochroma- 
tor  and  analyzer  use  triple  reflections  before  and  after  the  sample. 
Using  large  channel-cut  silicon  crystals  suppresses  the  "wings" 
of  the  beam  profile,  improving  the  signal-to-noise  ratio  to  values 
comparable  to  that  obtained  from  the  pinhole  instruments.  The  PCD 
will  cover  a  Q-range  from  0.0004  nnr'  to  0.1  nrrr1  with  an  expected 
beam  current  of  50,000  s '. 

Located  on  the  BT-5  beam  tube  in  the  Confinement  Building, 
the  layout  of  the  instrument  is  shown  in  Fig.  3.  A  vertically  and 
horizontally  focusing  graphite  premonochromator  directs  the  neu- 
tron beam  away  from  the  main  reactor  beam  towards  the  perfect 
crystal  monochromator.  The  monochromator  reflects  a  highly  col- 
limated  beam  to  the  sample.  After  the  sample,  a  high  precision 
rotation  stage  rotates  the  analyzer  crystal  to  scan  the  small  angle 
scattering  obtained  from  the  sample,  which  is  then  collected  by 
the  detector.  The  shields  supporting  the  premonochromator  are  set 
on  kinematic  mounts  for  accurate  repositioning.  The  monochroma- 
tor and  analyzer  are  isolated  from  room  vibrations  using  a  heavy 
pneumatic  vibration  isolation  table. 


Triple-bounce  Premonochromator 


Close-up  Sample  Collimator 
Position 

FIGURE  3.  Schematic  layout  of  the  BT-5  perfect  crystal  diffractometer. 
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In  1999,  all  detailed  design  work  was  finished  and  all  pur- 
chased parts  received.  The  beam  shutter  and  premonochromator 
shielding  have  been  installed.  Major  components  yet  to  be  installed 
are  the  detector  housing,  sample  position  table,  and  the  small  mono- 
chromator  shield.  Preliminary  characterization  measurements  lead- 
ing to  instrument  commissioning  are  planned  for  early  2000. 

THE  FILTER  ANALYZER  NEUTRON 
SPECTROMETER 

The  new  Filter  Analyzer  Neutron  Spectrometer  (FANS)  is  designed 
to  give  U.S.  researchers  access  to  unprecedented  sensitivity  for 
measuring  the  vibrational  spectra  of  a  wide  variety  of  materials. 
This  spectrometer,  which  will  replace  the  current  BT-4  instrument, 
uses  cryogenically  cooled  polycrystals  as  low  energy  band-pass 
filters  for  neutrons  scattered  from  the  sample.  A  dramatic  gain  in 
signal  over  the  existing  filter  analyzer  is  achieved  primarily  through 
a  large  increase  in  the  solid  angle  covered  by  the  secondary  spec- 
trometer. The  FANS  instrument  is  being  developed  in  two  phases. 
Phase  I  includes  the  first  of  two  new  filter  assemblies,  and  Phase  II 
includes  the  second  filter  analyzer  and  a  new  monochromator  and 
monochromator  drum  system. 


FIGURE  4.  The  FANS  filter  assembly.  The  first  Be  filter  is  in 
the  foreground  with  the  PG  filter  second  and  the  final  Be  filter  in  the 
background  just  before  the  detector  bank.  A  radial  collimator  is 
visible  in  the  middle  of  the  assembly. 


During  the  past  year,  the  detailed  design  of  the  first  filter 
wedge  was  completed,  the  vacuum  chamber  which  encloses  the 
Be  and  graphite  filters  was  delivered,  and  the  Be  and  graphite 
filter  wedges  were  assembled  (see  Fig.  4).  Recent  tests  demon- 
strated that  these  filters  can  be  cooled  to  cryogenic  temperatures 
which  is  necessary  to  maximize  the  performance  of  the  instrument. 
Furthermore  all  of  the  parts  for  the  undercarriage  for  Phase  I  have 
been  received,  assembled,  and  operated  successfully  under  load. 
Most  of  the  shielding  and  data  acquisition  system  have  also  been 
received.  Major  assembly  and  installation  of  Phase  I  will  begin  early 
in  2000. 

DEVELOPMENT  OF  AN  ACTIVE  DOUBLE 
FOCUSSING  MONOCHROMATOR  SYSTEM 

The  use  of  vertical  or  horizontal  focussing  monochromators  to 
increase  the  intensity  of  selected  neutrons  incident  on  the  sample 
is  well  known  and  widespread.  Less  common  are  systems  combin- 
ing these  features.  Often  the  machinery  employed  to  do  this  is 
cumbersome  and  results  in  extraneous  material  in  the  beam,  increas- 
ing the  background  from  unwanted  neutrons.  For  horizontal  focus- 
sing it  is  desirable  to  have  maximum  flexibility  by  adjusting  many 


FIGURE  5.  A  drawing  of  the  active  double  focussing  monochromator  system. 
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monochromator  elements  while  keeping  the  corresponding  mechani- 
cal components  away  from  the  beam  area.  By  contrast,  vertical 
focussing  can  be  confined  to  the  problem  of  adjusting  only  one 
parameter,  the  radius  of  curvature  of  the  array,  and  admits  a  cor- 
respondingly simpler  mechanical  solution. 

We  are  currently  developing  a  system  that  combines  com- 
pletely flexible  horizontal  focussing  with  vertical  focussing  that  is 
achieved  by  buckling  the  entire  system  of  separately  controllable 
vertical  crystal  arrays.  Cylindrical  curvature  is  to  be  accomplished 
by  compressing  the  system  of  specially  shaped  strips  upon  which 
the  individual  monochromating  crystals  are  mounted.  This  design 
eliminates  extraneous  material  in  the  beam,  which  should  result  in 
greatly  reduced  background. 

AN  ADVANCED  LIQUID  HYDROGEN 
COLD  SOURCE 

The  NIST  liquid  hydrogen  cold  source  has  completed  over  four 
years  of  service.  It  was  installed  with  three  goals:  at  least  double 
the  cold  neutron  intensity  with  respect  to  its  predecessor  (D,0 
ice);  operate  simply  and  reliably;  and  pose  no  safety  threat  to  the 
reactor  or  personnel.  It  has  successfully  met  or  exceeded  all  these 
goals.  The  cold  neutron  flux  increased  by  a  factor  of  4  to  6,  for 
wavelengths  in  the  range  of  0.2  nm  to  2  nm.  The  availability  of 
the  source  has  been  nearly  99  %  of  the  time  that  the  reactor  was 
available  (the  reactor  is  shut  down  if  the  source  is  inoperable).  And 
there  have  been  no  hydrogen  leaks,  nor  have  any  of  the  insulating 
vacuums  or  helium  containments  been  compromised. 

Even  as  Unit  1  was  being  installed  in  1995,  however, 
improvements  in  the  MCNP  model  (used  for  Monte  Carlo  simula- 
tions) of  the  NIST  reactor  were  pointing  toward  a  new,  but  more 
complicated  cryostat  assembly,  with  a  possible  additional  gain  of  a 
factor  of  two.  Better  coupling  between  the  reactor  fuel  and  the  cold 
source  can  be  achieved  by  expanding  the  cooling  jacket  into  the 
volume  now  occupied  by  the  insulating  vacuum,  so  that  it  partially 
surrounds  the  moderator  vessel.  In  this  way,  the  D20  coolant  also 
serves  as  an  extension  of  the  reactor  reflector.  Additional  changes 
in  the  new  source  are  based  on  our  operating  experience  with  the 
existing  LH2  source  and  extensive  MCNP  calculations.  Unit  2  will 
be  an  ellipsoidal  shell  rather  than  spherical,  it  will  have  an  evacuated 
inner  vessel,  rather  than  vapor-filled,  and  the  LH2  layer  will  be  up 


FIGURE  6.  NCNR  Mechanical  Engineering  Technician  Scott  Slifer  welds  the 
aluminum  moderator  vessel  for  the  advanced  hydrogen  cold  source. 


to  30  mm  thick,  rather  than  20  mm,  without  increasing  the  LH2 
volume. 

The  advanced  liquid  hydrogen  cold  source  is  currently  being 
fabricated  and  will  be  installed  in  the  NIST  reactor  next  year. 
Enhanced  mechanical  design  and  manufacturing  tools  are  exploited 
in  the  fabrication  of  the  advanced  source.  Components  of  the  hydro- 
gen, insulating  vacuum,  helium  containment,  and  D20  vessels  are 
cut  from  solid  blocks  of  Al  6061  on  a  computer-controlled,  high- 
speed mill  at  the  NIST  Instrumentation  shop,  and  are  then  welded 
and  thoroughly  tested  by  NCNR  staff  (see  Fig.  6).  It  is  expected  that 
the  flux  of  cold  neutrons  will  increase  by  a  factor  of  1.8. 

IMPROVEMENTS  TO  NCNR  SAMPLE 
ENVIRONMENT  EQUIPMENT 

The  sample  environment  equipment  has  seen  a  number  of 
changes  in  FY  1999  at  the  NCNR.  One  of  the  most  visible  accom- 
plishments has  been  the  development  of  informative  webpages 
that  detail  sample  environment  resources,  specifications,  and  even 
the  current  operating  condition  and  location  for  specific  devices. 
This  information  is  easily  accessible  through  the  NCNR  homepage 
(http://www.ncnr.nist.gov/),  allowing  guest  researchers  to  plan  their 
experiments  and  design  appropriate  sample  holders. 
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INSTRUMENTATION  DEVELOPMENTS 


Commissioning  has  begun  on  a  new  high  magnetic  field/ 
low  temperature  superconducting  magnet  system  which  is  financed 
through  a  joint  collaboration  of  Johns  Hopkins  University,  NCNR, 
NEC  Research  Institute,  Princeton,  the  University  of  Maryland,  and 
a  National  Science  Foundation  IMR  grant.  The  0.022  K  dilution 
refrigerator  of  this  new  system  has  already  been  tested  and  used  on 
triple  axis  spectrometers.  A  7  T  magnet  has  been  temporarily  outfit- 
ted in  the  system  at  this  time,  with  replacement  by  a  much  stronger 
superconducting  magnet  of  10  T  to  12  T  in  the  coming  months. 
The  magnet  is  capable  of  operating  with  either  a  homogeneous 
magnetic  field  at  the  sample  position  or  with  a  field  gradient,  a 
great  aid  for  polarized  neutron  beam  experiments.  This  system  is 


Table  1  Recent  Sample  Environment  Acquisitions 


Helium  flow  cryostat 

1.5  K  to  300  K,  dedicated  to 

backscattering  spectrometer 

Helium  flow  cryostat 

1 .5  K  to  300  K,  dedicated  to  disk  chopper 

time-of-flight  spectrometer 

Helium  flow  cryostat 

1.5  Kto  300  K 

Pumped  helium-3  cryostat  insert 

0.30  K  to  300  K,  for  use  with  7  T 

vertical  field  magnet 

Closed-cycle  helium  refrigerator 

7  K  to  320  K 

Closed-cycle  helium  refrigerator 

10  Kto  320  K,  modified  for 

backscattering  spectrometer 

Poiseuilte  flow  shear  cell 

shear  rate  up  to  130,000  s1  near  surface 

Light  scattering  particle  sizer 

2  nm  to  100  nm  diameter 

Rheometer 

1.7 x  103  Pasto2.7x108  Pas, 

-60  °C  to  600  °C,  for  in-situ  measurement 

during  SANS  experiments 

Lyophilizer 

freeze  dryer  for  sample  preparation 

Kare  Fisher  titrator 

quantitative  analysis  of  water  in  a  sample 

Freezer 

-20  °C,  for  longer  term  storage  of 

biological  samples 

Superconducting  magnet 

OA  to  120  A,  bi-directional 

power  supply 

Dual  channel  lock-in  amplifier 

1  mHz  to  1 02  kHz  frequency  range 

Four  temperature  controllers 

programmable,  remote  operation 

Dual  channel  lock-in  amplifier 

OA  to  120  A,  bi-directional 

Four  temperature  controllers 

programmable,  remote  operation 

Single-crystal  windowed  tail 

for  SANS  low-temperature 

electromagnet  experiments 

Sensitive  equipment 

custom-designed  to  protect 

transport  carts 

valuable  equipment 

top  loading,  which  provides  an  important  capability  for  experiments 
requiring  multiple  sample  changes  or  for  quick  changeover  between 
consecutive  experiments. 

From  new  capabilities  to  improved  sample  preparation  and 
screening,  the  sample  environment  resources  have  been  significantly 
expanded  during  the  past  year.  Listed  in  Table  1  are  further  acquisi- 
tions, notable  modifications,  and  improvements  that  combine  to 
offer  the  researcher  more  control  over  the  sample  conditions.  Fig. 
7  displays  a  typical  cryorefrigerator  sample  environment.  More 
details  are  available  at  the  web  site  mentioned  above. 


NCNR  Sample  Environment  Geometry 
Cryorefrigerators 


First  Stage 


Aluminum 


FIGURE  7.  Schematic  diagram  of  the  sample  environment  in  a 
cryorefrigerator  in  the  neutron  beam. 
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INSTRUMENTS  AND 
CONTACTS 


High  resolution  powder  diffractometer  (BT-1) 

J.  K.  Stalick,  (301)  975-6223,  judy.stalick@nist.gov 

B.  H.  Toby,  (301)  975-4297,  brian.toby@nist.gov 

Residual  stress  and  texture  diffractometer  (BT-8) 

H.  J.  Prask,  (301)  975-6226,  hank@nist.gov 
P.  C.  Brand,  (301)  975-5072,  pbrand@nist.gov 

30  m  SANS  instrument  (NG-7) 

C.  J.  Glinka,  (301)  975-6242,  cglinka@nist.gov 
J.  B.  Barker,  (301)  975-6732,  barker@nist.gov 
S.-M.  Choi,  (301)  975-4863,  sungmin@nist.gov 

30  m  SANS  instrument  (NG-3)  (NIST/NSF-CHRNS) 

B.  Hammouda,  (301)  975-3961,  hammouda@nist.gov 
S.  R.  Kline,  (301)  975-6243,  steven.kline@nist.gov 
S.  Krueger,  (301)  975-6734,  krueger@nist.gov 

8  m  SANS  instrument  (NG-1) 

A.  Karim,  (301)  975-6588,  karim@nist.gov 

Perfect  crystal  SANS  (BT-5)  (NIST/NSF-CHRNS) 

J.  B.  Barker,  (301)  975-6732,  barker@nist.gov 

Cold  neutron  reflectometer  -  vertical  sample-polarized  beam  option  (NG- 

C.  F.  Majkrzak,  (301)  975-5251,  cmajkrzak@nist.gov 
J.  A.  Dura,  (301)  975-6251,  jdura@nist.gov 

Cold  neutron  reflectometer  -  horizontal  sample  (NG-7) 

S.  K.  Satija,  (301)  975-5250,  satija@nist.gov 
R.  Ivkov,  (301)  975-6737,  rivkov@nist.gov 

Triple-axis  polarized-beam  spectrometer  (BT-2) 

J.  W.  Lynn,  (301)  975-6246,  jeff.lynn@nist.gov 

Triple-axis  spectrometer  (BT-9) 

R.  W.  Erwin,  (301)  975-6245,  rerwin@nist.gov 
P.  M.  Gehring,  (301)  975-3946,  pgehring@nist.gov 

Spin-polarized  triple-axis  spectrometer  (NG-5)  (NIST/NSF-CHRNS) 

P.  M.  Gehring,  (301)  975-3946,  pgehring@nist.gov 
S.  -H.  Lee,  (301)  975-4257,  seung-hun.lee@nist.gov 

Filter-analyzer  neutron  spectrometer  (BT-4) 

T  J.  Udovic,  (301)  975-6241,  udovic@nist.gov 

P.  Papanek,  (301)  975-5049,  pete@jazz.ncnr.nist.gov 

D.  A.  Neumann,  (301)  975-5252,  dan@nist.gov 


Fermi-chopper  time-of-flight  spectrometer  (NG-6) 

C.  M.  Brown,  (301)  975-5134,  craig.brown@nist.gov 
T.  J.  Udovic,  (301)  975-6241,  udovic@nist.gov 

Disk-chopper  time-of-flight  spectrometer  (NG-4) 

J.  R.  D.  Copley,  (301)  975-5133,  jcopley@nist.gov 
J.  C.  Cook,  (301)  975-6403,  jeremy.cook@nist.gov 

High-flux  backscattering  spectrometer  (NG-2) 

R.  M.  Dimeo,  (301)  975-8135,  robert.dimeo@nist.gov 

D.  A.  Neumann,  (301)  975-5252,  dan@nist.gov 

Neutron  spin  echo  spectrometer  (NG-5) 

N.  S.  Rosov,  (301)  975-5254,  nrosov@nist.gov 

Prompt-gamma  neutron  activation  analysis  (NG-7) 

R.  M.  Lindstrom,  (301)  975-6281,  dick.lindstrom@nist.gov 
R.  L.  Paul,  (301)  975-6287,  rpaul@nist.gov 

Other  activation  analysis  facilities 

R.  R.  Greenberg,  (301)  975-6285,  rgreenberg@nist.gov 

Cold  neutron  depth  profiling  (NG-O) 

G.  Lamaze,  (301)  975-6202,  lamaze@nist.gov 

Instrument  development  station  (NG-O) 

D.  F.  R.  Mildner,  (301)  975-6366,  mildner@nist.gov 

H.  H.  Chen,  (301)  975-3782,  chenmayer@nist.gov 

Neutron  interferometer  (NG-7) 

M.  Arif,  (301)  975-6303,  muhammad.arif@nist.gov 

Fundamental  neutron  physics  station  (NG-6) 

M.  S.  Dewey,  (301)  975-4843,  mdewey@nist.gov 

Theory  and  modelling 

N.  F.  Berk,  (301)  975-6224,  nfb@nist.gov 
I  Yildirim,  (301)  975-6228,  taner@nist.gov 

Sample  environment 

D.  C.  Dender,  (301)  975-6225,  dender@nist.gov 

L.  J.  Santodonato,  (301)  975-6657,  lsantodonato@nist.gov 
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For  additional  information  on  the  facility, 
contact: 

J.  Michael  Rowe 

(301)  975-6210 
mike.rowe@nist.gov 

John  J.  Rush 

(301)  975-6231 
john.rush@nist.gov 

To  obtain  guidelines  for  preparing  proposals 
to  conduct  research  at  the  facility,  contact: 

William  Kamitakahara 

(301)  975-6878 
william.kamitakahara@nist.gov 

Location  of  all  contacts: 

National  Institute  of  Standards  and  Technology 
100  Bureau  Drive,  MS8562 
Gaithersburg,  MD  20899-8562 

Details  are  also  available  on  the 

IMCNR  Web  site:  http://www.ncnr.nist.gov/ 


